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ABSTRACT 

The  fatigue  failure  of  a  centre-box  lug.  found  in  the  Seahawk’s  stabilator, 
highlighted  deficiencies  in  the  fatigue  life  assessment  of  helicopter  airframe 
structure.  Thus  a  methodology  for  fatigue  life  assessment  was  developed  using 
this  cracked  lug  as  a  demonstrator.  Indirect  measurements  from  the  Flight 
Loads  Survey  on  the  Black  Hawk,  although  not  ideal,  allowed  the  estimation 
of  the  amplitude  and  frequency  of  the  lug’s  loading.  Of  the  sixteen  assumptions 
made  in  this  fatigue  analysis,  the  most  restrictive  was  that  a  high  amplitude 
and  high  frequency  loading  acts  for  the  entire  flying  time.  Three  different 
fatigue  lives  were  obtained  based  on  different  levels  of  conservatism  in  the 
loading  estimates.  These  three  loadings  were:  the  worst-case  scenario,  a  high 
loading  scenario,  and  a  best-case  scenario.  The  worst-case  and  high-loading 
scenarios  resulted  in  low  fatigue  lives,  while  the  best-case  scenario  resulted 
in  an  unlimited  life  for  the  cracked  lug.  It  was  surprising  to  find  such  low 
fatigue  lives  for  two  of  these  scenarios,  but  these  low  lives  may  be  due  to  the 
conservative  assumptions  used  in  the  analysis. 
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Fatigue  Life  Estimate  of  Centre-Box  Lug  in  Seahawk’s 

Stabilator 


EXECUTIVE  SUMMARY 

During  the  routine  folding  of  a  Seahawk  stabilator  in  June  2003,  a  crack  was  found 
in  the  centre-box  lug,  which  had  accrued  only  402  flight  hours.  Had  this  cracked  lug 
gone  unnoticed,  the  eventual  loss  of  the  stabilator  mid-flight  could  have  resulted  in  loss 
of  control — and  the  potential  loss  of  the  crew  and  aircraft.  The  simplified  fatigue  analysis 
undertaken  in  this  report  investigates  the  possibility  of  this  fatigue  cracking  occurring  in 
other  Seahawk  helicopters  operated  by  the  Australian  Defence  Force  (ADF).  However, 
the  primary  intention  was  not  to  define  an  alternative  fatigue  life,  but  to  (i)  explore  a 
methodology  for  and  (ii)  develop  DSTO’s  capability  in  assessing  the  fatigue  life  of  helicopter 
airframe  structure.  In  other  words,  the  fatigue  analysis  undertaken  on  the  lug  was  used 
as  a  demonstrator  of  this  developed  methodology. 

A  significant  limitation  was  that  only  indirect  loading  information  was  available  for 
the  Seahawk’s  stabilator.  We  used  the  measurements  taken  in  2000  during  the  Flight 
Loads  Survey  for  the  Black  Hawk  as  the  indirect  information  to  estimate  the  lug  loading. 
The  indirect  nature  of  the  loading  measurements  meant  that  the  fatigue  analysis  required 
sixteen  assumptions.  The  two  coarsest  assumptions  were:  (1)  The  estimated  high  loading 
acted  for  the  entire  flight  time.  (2)  The  cycle  counting  (for  fatigue  purposes)  could  be 
estimated  from  the  frequency  decomposition  of  the  loading’s  time-history. 

The  fatigue  analysis  relied  upon  bending  stresses  resulting  from  stabilator  buffeting 
during  flight.  These  bending  stresses  were  measured  by  a  strain  gauge  bridge  that  was 
located  9.2  inches  from  the  cracked  lug  and  directly  above  the  forward  spar  on  the  stabi- 
lator’s  skin.  Thus  the  loading  on  the  cracked  lug  itself  had  to  be  calculated  indirectly  by 
assuming  a  particular  lifting  distribution  on  the  stabilator. 

The  lug’s  fatigue  life  was  calculated  for  three  cases,  which  included  a  case  termed  the 
worst-case  scenario.  This  worst-case  considered  fatigue  under  a  “harsh”  loading  environ¬ 
ment  and  resulted  in  a  short  fatigue  life.  This  short  life  was  of  concern  because  the  cracked 
lug  was  life  unlimited.  We  would  have  expected  that  such  a  conservatively  designed  lug 
(one  with  an  unlimited  life)  would  yield  a  long  fatigue  life  even  under  the  worst-case  sce¬ 
nario.  In  fact,  the  opposite  was  found— namely,  a  short  component  life — but  this  result 
is  possibly  due  to  the  two  coarse  assumptions  listed  above.  The  analysis  in  this  report 
suggests  that  fatigue  might  be  a  problem  in  other  Seahawks  in  the  ADF’s  fleet.  However, 
a  full  fatigue  analysis  (using  the  suggested  refinements  made  in  this  report)  would  be 
required  to  further  quantify  the  fatigue  susceptibility  of  the  centre-box  lug  in  question. 

The  load  path  from  the  stabilator  to  the  centre-box  is  redundant,  and  so  a  crack 
detection  (instead  of  crack  prevention)  program  would  be  sufficient  for  safety.  Within  the 
fidelity  of  results  obtained  so  far,  the  work  carried  out  in  this  report  does  not  conflict  with 
Sikorsky’s  assessment  of  an  unlimited  fatigue  life. 
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Notation 


Roman  Symbols 

a  crack  length 

A*i  cross-sectional  area  of  zth  rectangular  region  (where  the  subscript  V’ 
denotes  either  “a”  or  “f  ’  respectively  for  aft  or  forward  spar) 

Ac  cross-sectional  area  of  cracked  lug 

C  coefficient  in  law  for  simple  crack  growth  rate 
Ft  u  tensile  ultimate  strength 

Fiy  tensile  yield  strength 

h  distance  from  neutral  axis  to  top  of  beam 
height  of  half  the  aft  spar 
hf  height  of  half  the  forward  spar 

H  horizontal  load 

H  frequency  response  function 

Hai  horizontal  load  on  zth  rectangular  region  of  aft  spar 

Ha  horizontal  load  on  zth  rectangular  region  of  forward  spar 

Hu  ultimate  horizontal  load  of  lug 

i  counting  variable  for  crack  striations 
i  counting  variable  for  enumeration  of  run  number 

ft  imaginary  unit  vector  (that  is,  yf—1) 

I  second  moment  of  area 

j  counting  variable  used  in  Fourier  transform  (time  domain) 

k  counting  variable  used  in  Fourier  transform  (frequency  domain) 

Kt  stress  concentration  factor 

AK  range  of  stress  intensity  factor 

L  length  of  stabilator  spar 

m  characteristic  slope  for  curve  of  crack  growth  rate 
M  bending  moment  for  free-body  diagram 

Mg  bending  moment  at  gauge 

Mr  bending  moment  at  root 

M*  bending  moment  at  root,  not  including  lift  between  gauge  and  root 

n  number  of  points  in  time-history  vector  for  Fourier  transform 
N  loading  cycles 

jVf  cycles  to  fatigue  failure 

p  lifting  distribution  acting  on  stabilator’s  spar 
P  maximum  of  elliptic  lifting  distribution 
R  stress  ratio  (for  use  in  cycles  to  failure  equation) 
s  direct  stress  measurement 

s0%  minimum  (Oth  percentile)  of  directly  measured  stress 

ss%  5th  percentile  of  directly  measured  stress 

s95%  95th  percentile  of  directly  measured  stress 

S\oo%  maximum  (100th  percentile)  of  directly  measured  stress 

S\  turning  point  at  start  of  loading  cycle 

S‘2  turning  point  at  end  of  loading  cycle 

Seq  equivalent  stress  (for  use  in  cycles  to  failure  equation) 

notation  continued  on  next  page  . . . 
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. . .  notation  continued  from  previous  page 

Sniax  maximum  stress  (for  use  in  cycles  to  failure  equation) 

Sm in  minimum  stress  (for  use  in  cycles  to  failure  equation) 

u  vector  containing  time-history  data  for  Fourier  transform 
v  vector  containing  frequency  data  for  Fourier  transform 
V  vertical  force  for  free-body  diagram 
Vmax  maximum  vertical  force  on  stabilator  wing  panel 
x  coordinate  of  crack  growth  striation 
x  spanwise  distance  measured  from  stabilator’s  root 
xg  spanwise  location  of  strain  gauge 

X  frequency  input  function 

y  coordinate  of  crack  growth  striation 

y  vertical  distance  of  a  fibre  from  neutral  axis 
y?  distance  between  top  and  bottom  lugs  in  spar 
y  frequency  output  function 

z  coordinate  of  crack  growth  striation 

Greek  Symbols 

a  non-dimensional  length  coefficient  (used  in  range  of  stress  intensity  factor) 

7  correction  factor  (compensates  for  lift  ignored  between  gauge  and  root) 

?/  fillet  height  on  forward  spar’s  lug 

Aa  centroidal  distance  of  aft  spar’s  lug 

Af  centroidal  distance  of  fillet  for  forward  spar’s  lug 

d  angle  from  horizontal  made  by  lug  of  aft  spar 

v  vibratory  component  of  stress 

z/95^  95 %  percentile  vibratory  component  of  stress 

zzaVg  average  vibratory  component  of  stress 

z/max  maximum  vibratory  component  of  stress 

£  distance  from  stabilator  wing  panel  root 

a  steady  component  of  stress 

<ra  bending  stress  on  aft  spar  of  stabilator 

<7aVg  average  steady  component  of  stress 

err  bending  stress  on  forward  spar  of  stabilator 

<jg  bending  stress  at  gauge 

doc  principal  stress 

<rmax  maximum  steady  component  of  stress 
ermiU  minimum  steady  component  of  stress 
<j a  stress  at  lug 

i  frequency  of  worst-case  loading 

cj)me(\  frequency  of  typically  high  loading 

X  constant:  first  component  of  bending  moment  at  gauge 

■ij)  constant:  second  component  of  bending  moment  at  gauge 

u)  frequency 

ul}1v  main  rotor  frequency 

u>t  tail  rotor  frequency 

notation  continued  on  next  page  . . . 
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. . .  notation  continued  from  previous  page 

Miscellaneous  Symbols  (in  the  following  notation  a  is  a  dummy  variable) 
a  vector 

aav g  subscript  for  average 

aa  subscript  for  aft  spar 

af  subscript  for  forward  spar 

ag  subscript  for  bending  bridge  gauge 

at  subscript  for  cracked  lug 

ar  subscript  for  spar’s  root 

awp  subscript  denoting  lugs  from  stabilator’s  wing  panel 
R  real  field 
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1  Introduction 

Following  flying  operations  on  the  29th  of  June  2003,  a  potentially  significant  crack 
was  found  on  a  Seahawk  helicopter.  This  crack  occurred  on  the  stabilator’s  centre-box  (see 
Figure  1.1)  and  was  found  during  a  routine  folding  of  the  stabilator.  Had  this  crack  gone 
unnoticed  further  fractures  would  have  developed.  Due  to  these  additional  fractures,  the 
loss  of  the  stabilator  mid-flight  could  have  meant  a  loss  of  control  for  the  helicopter— and 
the  potential  loss  of  the  crew  and  the  aircraft. 

The  Australian  Defence  Force’s  (ADF)  response  to  this  lug  failure  was  guided  by 
advice  [17]  from  the  Rotary  Wing  Section  of  the  Director  General  Technical  Airworthiness. 
This  advice  concluded  that  the  damage  was  of  a  fail  safe  nature  as  intended  by  the 
original  equipment  manufacturer  (which  was  the  Sikorsky  Aircraft  Corporation).  Hence 
the  recommendation  [17]  was  for  a  minor  amendment  of  the  inspection  program  to  detect 
the  unlikely  fracture  of  this  lug. 

The  work  carried  out  in  this  report  investigates  the  fatigue  loading  the  lug  would 
probably  experience,  and  hence  the  possibility  of  this  type  of  fatigue  cracking  occurring  in 
other  Seahawks.  However,  the  primary  intention  is  not  to  define  an  alternative  fatigue  life, 
but  to  explore  a  methodology  for  assessing  the  fatigue  life  of  helicopter  airframe  structure, 
for  which  there  are  currently  no  robust  methodologies.  Furthermore,  the  development  of 


Figure  1.1:  Line  drawing  of  the  starboard  side  stabilator  and  centre-box 
from  the  Seahawkfs  tail  An  arrow  points  to  the  cracked  lug,  which  is  the 
aft  portion  of  the  clevis  lug.  This  lug  is  located  on  the  top  starboard  side 
of  the  forward  spar .  The  two  remaining  arrows  point  to  the  forward  and 
aft  spars.  (Drawings  adapted  from  Black  Hawk  manuals  [24,  25].) 
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this  methodology  will  improve  DSTO’s  capability  in  the  fatigue  life  assessment  of  helicopter 
airframe  structure.  As  a  demonstrator,  we  develop  this  methodology  using  the  lug  failure 
in  the  Seahawk  stabilator  together  with  the  extensive  data  that  DSTO  acquired  from 
Flight  Loads  Survey  [3]  of  the  Black  Hawk. 

The  fatigue  characteristics  of  a  Seahawk  lug  are  investigated  using  simplified  proce¬ 
dures.  First,  a  brief  description  of  this  cracked  lug  is  given,  the  philosophy  behind  the 
fatigue  analysis  explained,  and  the  assumptions  made  throughout  this  report  summarised. 
In  order  to  calculate  the  stress  on  this  lug,  the  geometry  of  the  loading  spars  and  gauge 
locations  are  then  given.  The  horizontal  and  vertical  loading  on  these  spars  are  then  ap¬ 
proximated  using  skin  bending  stresses.  These  bending  stresses  were  obtained  from  the 
Flight  Loads  Survey  carried  out  on  a  Black  Hawk  helicopter.  The  loading  frequency  and 
the  stress  on  the  lug  are  then  estimated,  and  the  fatigue  properties  of  the  cracked  lug- 
approximated.  The  final  section  summarises  the  finding  made  within  this  report. 


1.1  Description  of  Cracked  Lug  and  Design  Review 

The  crack  occurred  in  the  receiving  clevis  lug  (only  the  aft  portion)  on  the  starboard 
side  of  the  stabilator  centre-box.  (For  brevity,  this  lug  will  be  called  the  cracked  lug  in  this 
report.)  From  new,  this  lug  component  had  experienced  a  mere  402  flight  hours,  which  is 
below  the  more  than  3000  flight  hours  the  fleet  leader  has  accumulated.  Figure  1.1  shows 
the  crack  in  relation  to  the  tail  of  the  aircraft,  progressively  zooming-in  on  the  cracked 
lug.  The  centre-box  and  starboard  stabilator  are  shown  in  this  figure,  and  an  arrow  marks 
the  cracked  lug. 

As  can  be  seen  from  Figure  1.1,  the  cracked  lug  is  part  of  a  clevis,  which  is  one  of  a  set 
of  four  clevis  joints  that  attach  the  outboard  portion  of  the  stabilators  to  the  centre  box. 
This  design  has  a  level  of  inherent  damage  tolerance  because  the  failure  of  any  single  lug 
will  not  result  in  the  loss  of  the  stabilator  wing.  The  clevis  design  maintains  independent 
and  redundant  load  paths  [1,  p.  H-13]  by  preventing  cracks  in  one  lug  from  propagating 
to  the  other  lugs.  The  design  of  the  centre-box  lugs  considered  only  static  loads,  that  is, 
fatigue  loading  was  not  considered. 

Figure  1.2  shows  a  top  view  of  the  centre-box  (bottom  left)  attached  to  the  starboard 
stabilator  (right)  and  tail  pylon  (top  left).  The  white  arrow  marks  the  cracked  lug  and 
the  hollow  arrow  to  the  left  of  this  photograph  defines  the  forward  direction.  Both  sides 
of  the  lug’s  aft  portion  were  completely  cracked  and  only  the  lug’s  bolt  was  holding  the 
outer  tip  of  the  lug  in  place. 

A  DSTO  minute  by  Byrnes  [5]  and  a  DSTO  investigation  report,  also  by  Byrnes  [6], 
describe  this  cracking  problem  further.  Byrnes’  forensic  investigation  report  also  details 
the  folding  property  of  the  Seahawk’s  stabilator  and  provides  more  detailed  photographs 
of  the  lug  fracture. 

Figure  1.3  shows  the  centre-box  detached  from  the  stabilator,  and  a  red  circle  highlights 
the  cracked  lug.  A  close-up  of  the  lug’s  broken  portion  is  shown  along  with  a  ruler  on  the 
bottom  photograph  of  this  figure.  This  bottom  photograph  shows  striation  marks,  which 
are  typical  of  fatigue  cracking  (see  Byrnes  [6]). 
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Figure  1.2:  The  aft  portion  of  the  clevis  lug  (centre  of  photograph)  has 
completely  cracked,  only  the  bolt  is  holding  the  outer  tip  of  the  lug  in 
place.  Portions  of  the  stabilator  (right),  tail  pylon  (top  left),  and  centre- 
box  (bottom  left)  are  shown  in  this  photo.  (Photograph  courtesy  of  Jamie 
Edwards). 


The  stabilator  assembly,  composed  of  the  stabilator  wing  panels  and  centre-box,  was 
fatigue  tested  by  Sikorsky.1  According  to  Sikorsky  [1],  steel  lugs  attach  the  wing  panels 
to  the  centre-box;  it  is  unclear  whether  this  “steer  reference  is  simply  a  typing  error.  If 
the  testing  was  indeed  conducted  on  a  stabilator  with  steel  lugs,  then  the  aluminium2 
centre-box  lugs  found  in  the  ADF’s  Seahawks  may  not  have  been  tested  by  Sikorsky  for 
fatigue. 

The  crack  initiation  phase  of  Sikorsky’s  fatigue  testing  was  conducted  at  a  constant 
vibratory  stress  determined  from  high-speed  level  flight.  When  a  crack  was  detected  in 
the  centre-box,  the  load  was  reduced  to  a  conservative  loading  simulation.  One  hour  of 
this  loading  simulation  consisted  of  62  000  cycles  of  flight  loading  (namely,  high-speed  level 
flight)  followed  by  4  ground-air-ground  (GAG)  cycles.  Sikorsky  define  a  GAG  cycle  as  the 
maximum  single  stress  cycle  (steady  plus  vibratory)  occurring  in  a  complete  flight.  To 
assess  damage  tolerance,  the  testing  continued  beyond  crack  detection  until  the  stabilator 
was  no  longer  able  to  sustain  the  testing  load. 

For  helicopters,  a  working  fatigue  curve  is  obtained  by  reducing  the  mean  fatigue  curve, 
which  is  achieved  by  shifting  the  curve  down  and  to  the  left  on  the  S-N  plane.  The  down¬ 
shift  is  obtained  by  multiplying  the  stresses,  on  the  mean  fatigue  curve,  by  a  reducing 

1  See  in  particular  Section  H.  pages  8-13,  of  this  Sikorsky  report  [1]. 

2For  a  discussion  of  the  material  properties  of  the  cracked  lug  see  Section  5.1  on  page  28. 
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Figure'  1.3:  Stahilator  centre -box  (top  photograph)  and  close-up  of  lug's 
broken  portion  (bottom  photograph).  In  the  s t a b ila to r  p ho tograpl i,  the  red 
circle  highlights  the  cracked  lug.  In  the  close-up  photograph,  the  ruler's 
units  arc  centinietrcs.  (Photographs  courtesy  of  Rohan  Byrnes.) 
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factor.  Similarly,  the  left-shift  is  obtained  by  multiplying  the  cycles  to  failure,  on  the 
mean  fatigue  curve,  by  a  different  reducing  factor.  The  results  described  in  the  following 
two  paragraphs  use  multiplying  factors  of: 

•  0.61  and  0.80  for  the  down-shift  and 

•  1/5  and  1/3  for  the  left-shift ,  respectively. 

For  further  details  about  how  to  develop  a  reduced  (or  working)  fatigue  curve  from  the 
mean  fatigue  curve  see,  for  example,  Lombardo  [18,  p.  32]. 

In  the  stabilator’s  wing  panels,  it  was  the  forward  spar  that  failed  in  this  fatigue  test. 
This  forward  spar  (on  the  right-hand  panel)  failed  in  cracking  without  chafing.  The  sta- 
bilator  wing  panels  were  given  component  retirement  lives  of  1600  hours  and  18000  hours 
for  respectively  a  0.61  and  0.80  reduction  of  the  mean  fatigue  curve.  The  damage  tolerant 
testing  showed  that  after  this  spar  cracking  was  detected,  the  stabilator  could  endure  a 
further  250  hours  of  equivalent  mission  flight  before  failure. 

In  the  centre-box,  it  was  a  rivet  hole  on  the  top  flange  of  the  aft  fitting  that  cracked. 
This  rivet  hole  was  located  on  the  spar  that  connects  the  left  and  right  aft  lugs  of  the  centre- 
box.  The  centre-box  was  given  component  retirement  lives  of  660  hours  and  6100  hours 
for  respectively  a  0.61  and  0.80  reduction  of  the  mean  fatigue  curve.  The  damage  tolerant 
testing  showed  that  after  this  rivet  crack  was  detected,  the  stabilator  could  endure  a 
further  160  hours  of  equivalent  mission  flight  before  failure. 

Sikorsky  concluded  that  the  fail-safety  of  the  stabilator  had  been  demonstrated,  and 
recommended  the  stabilator  assembly  be  replaced  on-condition  [1]. 


1.2  Philosophy  behind  Fatigue  Analysis 

The  analysis  carried  out  in  this  report  involved  the  use  of  several  quick- and- dirty 
methods,  which  required  several  assumptions.  The  nature  of  this  analysis  was  made 
necessary  by  either:  (i)  a  lack  of  information  or  (ii)  a  requirement  for  a  timely  solution. 
The  lack  of  information  meant  that  there  was  often  no  benefit  in  using  a  more  sophisticated 
or  complicated  analysis. 

As  a  first  order  approximation,  we  were  trying  to  assess  the  worst-case  and  best-case 
scenarios  (which  are  defined  below)  for  the  fatigue  cracking  of  the  lug. 

Worst- Case  Scenario:  assumes  the  helicopter  is  continuously  operated  in  a  “harsh” 
manner.  If  the  estimated  stresses  experienced  by  the  lug  under  this  harsh  scenario 
were  below  the  lug’s  run-out  stress,  then  we  could  effectively  rule  out  fatigue  cracking 
as  a  problem  in  the  centre-box  lugs  of  ordinary  Seahawk  stabilators.  This  scenario 
is  essentially  a  bound  from  above  on  stress . 

Best-Case  Scenario:  assumes  the  helicopter  is  continuously  operated  in  a  “gentle”  man¬ 
ner.  If  the  estimated  stresses  experienced  by  the  lug  were  significant  in  terms  of  fa¬ 
tigue,  then  we  could  conjecture  that  centre-box  lugs  in  ordinary  Seahawk  stabilators 
might  be  at  risk  from  fatigue.  This  scenario  is  essentially  a  bound  from  below  on 
stress . 


5 


DSTO-TR-1590 


It  is  important  to  note  at  this  stage  that  no  conclusion  (or  at  least  only  a  tentative 
conclusion)  could  be  drawn  if  we  obtained  either  of  the  two  opposing  scenarios;  namely, 
(i)  high  stresses  under  the  harsh  operations  assumption  or  (ii)  low  stresses  under  the  gentle 
operations  assumption.  These  two  scenarios  will  be  termed  the  null  conclusions . 

To  explain  these  null  conclusions  more  concretely,  let  us  consider  the  first-  case  of  high 
stresses  under  harsh  operations.  If  we  obtained  this  null  conclusion,  then  we  could  not  be 
sure  whether  these  high  stresses  were  due  to  (i)  a  possible  operating  environment  or  (ii)  an 
over-pessimistic  view  of  the  operating  environment.  In  contrast,  we  could  confidently  say 
that  fatigue  was  not  a  problem  if  we  obtained  low  stresses  under  the  worst-case  scenario 
(namely,  harsh  operations).  An  analogous  argument  can  be  made  for  the  best-case  scenario 
and  its  associated  null  conclusion. 

To  summarise,  this  quick  analysis  was  really  a  process  of  elimination  rather  than  a 
complete  solution.  Throughout  this  report  suggested  improvements  to  the  analysis  are 
made;  along  with  some  comments  on  the  added  complexity  required  and  possible  accuracy 
gained  in  using  these  suggested  improvements. 

In  addition  to  these  extreme  cases,  an  intermediate  loading  case  is  also  analysed.  This 
intermediate-case  estimates  the  lug  loading  by  using  high  loads  from  a  level  flight  manoeu¬ 
vre,  in  other  words,  a  typical  high-loading  case.  Unlike  the  worst-case,  this  intermediate- 
case  is  not  necessarily  conservative,  and  hence  it  is  used  more  as  a  check  on  the  level  of 
conservatism  inherent  in  the  worst-case  scenario  than  as  a  solution  in  its  own  right. 


1.3  Synopsis  of  Assumptions 

In  this  section  we  summarise  the  assumptions  made  during  the  analysis  of  the  lug 
cracking  problem.  The  assumptions  in  this  report  were  made  for  one  of  two  reasons: 


•  Lack  of  information:  There  was  insufficient  information  to  carry  out  a  full  analysis 
or  the  available  information  was  ambiguous. 

•  Expediency  of  the  solution:  In  order  to  obtain  a  quick- and- dirty  solution  a  full  anal¬ 
ysis  was  often  omitted,  and  instead  a  simplified  analysis  was  undertaken. 


In  this  report,  the  symbol  X  is  printed  in  the  left  margin  to  emphasise  the  introduction 
of  an  assumption. 

The  calculations  in  this  report  use  dimensions  of  the  load  bearing  components  obtained 
from  Sikorsky  maintenance  drawings  (see  Appendix  A  starting  on  page  51).  In  using  these 
Sikorsky  drawings  an  implicit  assumption  was  made: 


Assumption  1  The  horizontal  stabilators  of  the  Seahawk  are  identical  to  those  of 
the  Black  Hawk.  Or  at  least  from  a  loading  perspective,  the  structural  information  ob¬ 
tained  from  the  Sikorsky  drawings  are  a  good  representation  of  the  stabilator  structure 
found  on  the  Seahawk. 
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This  assumption  was  made  because  only  Black  Hawk  (and  not  Seahawk)  component 
drawings  were  readily  available.  The  Sikorsky  maintenance  drawings  of  the  Black  Hawk 
available  at  DSTO  appear  be  a  mixed  bag  of  drawings.  For  example,  there  are  two 
drawings  of  the  aft  stabilator  spar,  compare  Figure  A6  with  Figure  A7  (on  pages  57  and  58, 
respectively).  These  two  drawings,  of  what  should  be  the  same  component,  are  completely 
different:  the  angled  aft  spar  in  Figure  A7  probably  belongs  to  the  UH-60  as  compared 
to  the  straight  aft  spar  of  Figure  A7.  Similarly,  compare  the  top  stabilator  views  in 
Figures  A2  and  A3  (on  pages  53  and  54,  respectively). 

The  specific  assumptions  made  throughout  this  report  are  summarised  below.  More 
detailed  information  on  each  of  these  assumptions  is  given  on  the  page  where  the  as¬ 
sumption  is  first  introduced  (references  to  these  pages  are  given  next  to  the  assumption 
number). 

The  following  assumptions  were  made  during  the  calculation  of  lug  stress: 

•  Assumption  2  (p.  8):  how  the  bending  bridges  on  the  stabilator  are  mounted. 

•  Assumption  3  (p.  10):  the  dimensions  of  the  spar  caps  in  the  stabilator  are  constant. 

•  Assumption  4  (p-  12):  the  dimensions  of  spars  in  the  stabilator  are  estimated. 

•  Assumption  5  (p.  14):  torsional  and  chordwise  loads  in  the  stabilator  are  ignored. 

•  Assumption  6  (p.  16):  the  lifting  load  on  the  stabilator  is  elliptic. 

•  Assumption  7  (p.  23):  manoeuvres  within  a  group  are  similar. 

•  Assumption  8  (p.  28):  the  fatigue  properties  of  7075  aluminium  alloy  with  T6  and 
T7  temperings  are  similar. 

•  Assumption  9  (p.  29):  the  cycles  to  failure  equation  given  by  the  military  handbook 
for  metallic  materials  can  be  extrapolated. 

•  Assumption  10  (p.  31):  the  stabilator’s  percentile  stresses  characterise  fatigue  load¬ 
ing  in  the  cracked  lug. 

•  Assumption  11  (p.  33):  level  flight  stresses  are  typical  all  the  time. 

•  Assumption  12  (p.  34):  all  manoeuvres  consume  the  same  fraction  of  usage  time 
and  experience  the  same  stress. 

•  Assumption  13  (p.  35):  the  steady  and  vibratory  percentile  stresses  yield  conserva¬ 
tive  estimates  of  fatigue. 

•  Assumption  14  (p-  36):  the  Flight  Loads  Survey  is  representative  of  the  in-service 
flight  loads. 

•  Assumption  15  (p.  37):  the  typical  loading  frequency  dominates  fatigue  usage. 

•  Assumption  16  (p.  37):  all  manoeuvres  exhibit  the  typical  loading  frequency. 

Not  all  of  these  assumptions  were  used  for  the  three  loading  cases  (worst-,  best-, 
and  intermediate-cases).  For  example,  Assumption  11  is  not  required  for  the  worst-case 
scenario. 
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2  Stabilator’s  Geometry  and  Gauge  Locations 


In  order  to  determine  the  vertical  and  horizontal  loading  on  the  cracked  lug  we  must 
know  the  loading  on  the  stabilator.  The  only  information  readily  available  were  strain 
gauge  and  accelerometer  measurements  from  the  Flight  Loads  Survey  on  the  Black  Hawk. 
The  strain  gauge  information  came  in  the  form  of  bending  stress  on  the  upper  and  lower 
skins  of  the  stabilator  (directly  above  and  below  the  forward  and  aft  spars).  The  spars’ 
geometries  were  needed  to  convert  this  bending  stress  to  a  lug  loading.  Hence  this  section 
provides  information  on  both  the  gauge  locations  and  the  spars’  geometries. 


2.1  Stabilator  Strain  Gauge  and  Accelerometer  Locations 
during  Flight  Loads  Survey 

During  the  Flight  Loads  Survey  [3]  carried  out  on  a  Black  Hawk,  five  bending  stresses 
were  measured  on  the  both  the  port  and  starboard  stabilators  [8,  pp.  76-8].  For  each  of 
these  ten  stabilator  locations  (five  on  the  port  side  and  five  on  the  starboard  side),  four 
strain  gauges  were  configured  as  a  full-bending  bridge.  In  total,  forty  strain  gauges  were 
mounted  on  the  stabilators  (4  gauges  x  5  locations  x  2  stabilators  =  40  gauges).  Figure  2.1 
shows  where  these  five  bending  bridges  were  mounted  for  the  port  side  stabilator. 

Friend  [8]  does  not  state  how  these  bending  bridges  were  mounted,  and  so  we  make 
the  following  assumption: 

Assumption  2  Each  bending  bridge  mounted  on  the  stabilator  consisted  oj  four 
strain  gauges ,  two  gauges  on  the  forward  spar  and  two  gauges  on  the  aft  spar.  For 
both  the  forward  and  aft  spars,  a  strain  gauge  was  mounted  on  both  the  top  and  bottom 
skins  directly  above  and  beloiv  each  spar.  The  bending  stresses  on  the  front  and  aft 
spars  (obtained,  from  these  strain  gauges)  were  then  combined  to  obtain  an  average 
stabilator  bending  at  a  particular  butt  line.  A  positive  bending  measurement  implies  a 
compressive  stress  on  the  top  skin. 

Figure  2.1  shows  the  five  butt  line  locations  (note  that  distances  in  this  figure  are  shown 
from  the  inboard  edge  of  the  stabilator,  which  is  9.00"  from  the  helicopter’s  centreline).  We 
have  assumed  the  convention  that  a  positive  bending  measurement  implies  a  compressive 
stress  on  the  stabilator’s  top  skin;  and  hence  a  compressive  load  on  the  cracked  lug. 

In  Friend’s  report  [8],  these  bending  bridges  on  the  left-  and  right-hand  sides  were 
respectively  given  the  mnemonics  STBNBMiL  and  STBNBM?'R,  where  i  =  0, 1,2, 3,4  as  shown 
in  Figure  2.1. 

In  addition  to  these  strain  gauges,  there  were  also  four  accelerometers  mounted  on  the 
stabilator  and  three  strain  gauges  mounted  on  the  stabilator’s  centre-box. 

The  four  accelerometers  (two  on  each  side  of  the  stabilator)  were  mounted  on  the 
forward  spar  of  the  stabilator  approximately  37"  and  74"  from  the  inboard  edge  [8,  pp.  90- 
1],  see  Figure  2.1.  In  Friend’s  report,  these  accelerometers  on  the  left-  and  right-hand  sides 
were  respectively  given  the  mnemonics  VSTABxL  and  VSTAB.xR,  where  x  =  T,  M  as  shown  in 
Figure  2.1. 
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Figure  2.1:  Top  view  of  the  port  side  stabilator  showing  the  location  of 
the  strain  gauges  H  (for  i  =  0,1, 2, 3,4)  and  accelerometers  ©  (for 
x  =  T,  mounted  during  the  Flight  Loads  Survey.  Distances  shown 
are  from  the  inboard  edge  of  the  stabilator,  which  is  9.00"  from  the 
helicopter’s  centreline.  The  schematic  of  the  zoomed  section  represents 
a  perpendicular  (that  is,  chordwise)  cross-section  of  the  stabilator.  This 
cross-section  shows  the  assumed  location  of  the  strain  gauges  within  a 
bending  bridge.  (Stabilator  drawing  adapted  from  a  GTRI  report  [8, 
p.  78].) 
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The  three  centre-box  strain  gauges  were  mounted  in  a  45°  rosette  pattern  on  the 
horizontal  stabilator.  More  specifically  these  strain  gauges  were  mounted  on  the  forward 
spar’s  forward  face  of  the  centre-box.  In  Friend’s  report  [8,  pp.  76-7],  these  strain  gauges 
were  given  the  mnemonics  STBCBRly,  where  y  =  A,  B,  and  C  respectively  represent  the 
gauges  aligned  0°,  45°,  and  90°  from  the  vertical.  These  gauges  are  not  shown  in  Figure  2.1 
because  the  centre-box  is  not  included  in  this  illustration. 

2.2  Geometry  of  Forward  and  Aft  Stabilator  Spars 

The  approximate  dimensions  of  the  forward  and  aft  spars  were  obtained  from  Sikorsky 
maintenance  drawings.  Figure  2.2  shows  drawings  of  the  forward  and  aft  spars: 

•  Box  A  shows  the  middle  rib  (along  with  the  forward  and  aft  spars)  of  the  horizontal 
stabilator  at  BL  28.  This  box  is  a  zoomed  section  from  Sikorsky  drawing  70201-07052 
grid  reference  D8,  see  Figure  A9  (on  page  60)  for  the  complete  drawing. 

•  Box  B  shows  the  dimensions  of  the  top  cap  from  the  forward  spar  at  BL  9.  This 
box  is  a  zoomed  section  from  Sikorsky  drawings  70202-07051  grid  reference  K21,  see 
Figure  A5  (on  page  56)  for  the  complete  drawing. 

•  Box  C  shows  the  bottom  cap  of  the  aft  spar.  This  box  is  a  zoomed  section  from 
Sikorsky  drawings  70202-27001  grid  reference  B5,  see  Figure  A8  (on  page  59)  for  the 
complete  drawing. 

The  drawings  shown  in  Figures  A4  and  A10  (on  pages  55  and  61,  respectively)  may  also 
be  of  interest. 

It  was  unclear  from  these  drawings  how  the  spar  cap  dimensions  changed  with  spanwise 
location,  and  so  we  make  the  following  assumption: 

Assumption  3  In  both  the  forward  and  aft  spars  of  the  horizontal  stabilator,  the 
cross-sectional  dimensions  of  the  spar  caps  remained  constant  with  spanwise  location. 

This  assumption  was  made  for  expediency,  because  from  the  forward  and  aft  spar  drawings 
these  caps  do  taper.  In  particular,  the  caps  appear  to  taper  in  at  least  the  following 
ways:  (i)  the  width  of  the  spar  flanges  taper,  (ii)  the  thickness  of  the  spar  flanges  taper, 
(iii)  the  height  of  the  cap’s  web  tapers,  and  (iv)  the  boron  composite  (bonded  to  the  caps) 
tapers.  For  additional  details  on  cap  tapering  see  Figures  A5  and  A6  (on  pages  56  and 
57,  respectively). 

What  are  the  implications  of  Assumption  3?  Given  that  we  are  using  the  cap  dimen¬ 
sions  obtained  from  BL  9.00  (see  Assumption  4  below),  if  the  caps  were  tapered  then  the 
true  stresses  on  the  lug  would  be  lower  than  the  stress  estimations  obtained  in  this  report. 
The  reason  the  stress  would  be  lower  is  that  the  cross-sectional  area  of  the  caps  would  be 
smaller,  which  would  reduce  the  estimated  tensile  stress  (see  Section  3.3  for  more  details). 
Given  that  we  are  primarily  chasing  a  worst-case  scenario,  this  over-estimation  of  stress 
is  fortuitous  (but  would  probably  not  be  significant  anyway). 

Using  the  information  in  Figure  2.2,  the  caps  were  assumed  to  have  the  following 
dimensions  (obtained  from  the  top  cap  of  the  forward  spar  at  BL  9.00): 


Figure  2.2:  Caps  of  both  the  forward  and  aft  spars  on  the  horizontal 
stabilator.  Box  A  shows  the  middle  rib  (along  with  the  forward  and  aft 
spars)  at  BL  28.  Box  B  shows  the  top  cap  of  the  forward  spar  at  BL  9. 
Box  C  shows  the  bottom  cap  of  the  aft  spar  at  BL  1 7. 
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Forward  Spar  Aft  Spar 

Figure  2.3:  Estimated  dimensions  of  the  forward  and  aft  spars .  The 
heights  of  the  two  spars  are  given  as  a  range,  the  maximum  and  minimum 
represents  the  heights  at  BL  9  and  BL  66,  respectively.  The  labelled 
components  of  the  forward  spar  are:  (a)  the  cap’s  flange,  (b)  the  cap’s 
web,  (c)  the  spar’s  web,  (d)  the  boron  composite  stiffener,  and  (e)  the 
adhesive  bond. 


Assumption  4  The  dimensions  of  the  forward  and  aft  spars  are  as  shown  in  Figures  2.3. 
-  Both  the  forward  and  aft  spars  have  vertical  symmetry  about  the  shown  centreline.  The 

thickness  of  the  adhesive  bond  (which  is  0.010")  was  absorbed  into  the  thickness  of  the 
boron  composite  stiffener.  The  fillets  between  the  cap  ’ s  web  and  flange  were  ignored. 
The  top  of  the  spar’s  web  is  assumed  to  be  flush  with  the  bottom  of  the  cap’s  flange. 
The  web  of  the  aft  spar  is  the  same  thickness  as  the  web  of  the  forward  spar.  The  skins 
of  the  stabilators  were  ignored.  The  heights  of  both  spars  vary  linearly  with  spanwise 
location. 


The  dimensions  shown  in  Figure  2.3  were  obtained  from  the  spar  sections  shown  in  Figure  2.2. 

Due  to  a  lack  of  information,  the  aft  spar’s  web  was  assumed  to  be  the  same  thickness 
as  the  forward  spar’s  web.  Ignoring  the  fillets  and  skins  reduced  the  cross-sectional  area, 
and  hence  we  obtain  an  over-estimation  of  the  true  stress.  Conversely,  increasing  the 
height  of  the  spar’s  web  (so  that  it’s  flush  with  the  flange)  increases  the  cross-sectional 
area  marginally.  However,  the  overall  area  should  still  decrease  because  the  skin  and  fillets 
were  ignored,  and  hence  we  end  up  with  a  conservative  estimate  of  stress. 

It  was  unclear  from  the  available  drawings  whether  the  aft  spar  also  had  a  boron 
composite  stiffener  attached  to  it.  Given  that  the  aft  spar  was  considerably  thicker  than 
the  forward  spar,  we  assumed  there  was  no  additional  stiffener  on  top  of  the  aft  spar. 
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The  cap  height  of  the  aft  spar  shown  in  Figures  2.2  (which  is  1.121")  is  slightly  different 
to  the  corresponding  length  shown  in  Figure  2.3  (which  is  1.1 25").  This  difference  is  due 
to  the  5°  declination  of  the  cap.  which  increases  this  length  by  0.38%. 

The  heights  of  the  spars  were  measured  from  the  forward  and  aft  spar  drawings  (see 
Figures  A5  and  A6  on  pages  56  and  57,  respectively).  Only  two  spar  heights  were  taken 
on  each  of  these  drawings:  one  at  BL  9  and  the  other  at  BL  66.  Both  these  spar  drawings 
appear  to  suggest  that  the  spar’s  height  tapers  linearly  from  the  stabilators  inboard  edge 
(at  BL  9)  to  its  outboard  edge  (at  BL  82).  Thus  using  the  linearly  tapering  conjecture 
(made  in  Assumption  4),  we  can  calculate  the  spar  heights  at  the  five  bending  bridge 
locations  (see  Table  2.1). 


Table  2.1:  Stabilator  spar  heights  (both  forward  and  aft  spars)  at  the  five 
bending  bridge  locations. _ 


Bending  Bridge  Number 

- 

0 

1 

2 

3 

4 

- 

Butt  Line  (inch) 

9.00 

22.4 

18.2 

25.0 

37.5 

64.4 

66.0 

2h{  = 

Forward  spar  height  (inch) 

6.1 

5.8 

5.9 

5.7 

5.4 

4.7 

4.7 

to 

p 

1! 

Aft  spar  height  (inch) 

5.2 

4.8 

4.9 

4.8 

4.4 

3.6 

3.6 
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3  Stabilator  Loading 


Ill  this  section  the  vertical  and  spanwise  loading  on  the  spars  is  approximated  using 
the  bending  stress  information  available  from  the  Flight  Loads  Survey.  First  we  make 
some  assumptions  about  the  external  loading  on  the  stabilator.  The  vertical  force  on 
the  stabilator  is  then  estimated  assuming  a  particular  lifting  distribution,  which  is  scaled 
according  to  the  stabilator’s  bending  stress.  Finally,  the  horizontal  loading  (in  the  spanwise 
direction)  is  estimated  using  both  the  stabilator’s  bending  stress  and  the  vertical  loading. 


3.1  External  Loads  on  Stabilator 

The  external  loading  on  the  stabilator  spars  is  due  solely  to  aerodynamic  forces  acting 
on  the  stabilator’s  skin.  Hence,  the  external  forces  on  the  stabilator  act  in  both  the  vertical 
(due  to  lift)  and  chordwise  (due  to  drag)  directions. 

Due  to  the  limited  amount  of  stress  information,  we  make  several  assumptions  about 
how  the  load  is  transmitted  through  the  stabilator. 

Assumption  5  Both  the  forward  and  aft  spars  experience  only  vertical  loading  through 
their  respective  shear  centres  (that  is,  there  is  no  drag  load).  This  vertical  loading  is 
evenly  distributed  between  the  two  spars.  There  is  no  secondary  structure  (such  as 
ribs)  between  the  two  spars.  The  skins  do  not  take  any  loading  (not  even  in  shear). 

Normally  the  stabilator  would  act  as  a  shear  box,  however,  due  to  the  limited  amount  of 
information  available,  there  is  no  way  to  resolve  the  different  types  of  loading.  Drag  (that 
is,  chordwise  loading)  on  the  stabilator  was  also  ignored.  For  shapes  with  aerofoil  cross- 
sections,  and  provided  that  the  angle  of  attack  is  not  too  far  from  zero,  the  drag  is  typically 
an  order  of  magnitude  lower  than  lift  (see,  for  example,  Gerhart  and  Gross  [10,  p.  548]). 
The  design  loading  used  by  Sikorsky  [9]  supports  this  order  of  magnitude  assumption, 
so  ignoring  drag  is  a  reasonable  assumption.  In  essence,  Assumption  5  says  that  any 
loading  experienced  by  the  bending  bridges  mounted  on  the  stabilator  are  solely  due  to 
the  independent  vertical  loading  experienced  by  each  spar. 

Unfortunately,  using  only  the  Flight  Loads  Survey  measurements  very  little  (if  any) 
accuracy  can  be  gained  by  refining  Assumption  5.  We  saw  in  Assumption  2  (on  page  8) 
that  the  bending  stress  from  the  Flight  Loads  Survey  is  an  average  of  the  bending  stresses 
at  the  forward  and  aft  spars.  There  is  no  way  we  can  deconvolute  this  average  stress  into 
the  individual  bending  stresses  on  the  forward  and  aft  spars,  and  hence  no  overall  accuracy 
is  gained  from  the  refinement  of  Assumption  5. 

Consider  the  spars  as  determinant  structures,  which  is  itself  a  mild  assumption  because 
the  upper  and  lower  receiving  lugs  would  make  the  spar  statically  indeterminant.  In 
particular,  we  choose  to  model  them  as  simply  supported  beams.  Drawing  a  free-body 
diagram  through  one  of  the  gauges  results  in  the  diagram  shown  in  Figure  3.1.  The  cross- 
section  of  the  beam,  highlighted  by  the  dashed  circle,  shows  the  profile  of  the  tensile  stress 
resulting  from  the  bending  moment.  The  three-dimensional  cross-section  (on  the  right- 
hand  side  of  Figure  3.1)  shows  this  stress  profile  in  perspective.  Note  that  the  vertical  force 
V  does  not  affect  the  bending  bridge,  but  we  will  need  it  when  calculating  the  vertical 
loading  on  the  lugs. 
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Figure  3.1:  Using  a  free-body  diagram  to  calculate  the  loads  on  the  spar’s 
cross-section .  The  free-body  slice  is  taken  atx  =  xg ,  which  is  the  location 
of  the  strain  gauge.  The  distributed  load  p(x)  is  unknown,  and  is  the  only 
external  load  the  spar  experiences. 


3.2  Vertical  Load  on  Stabilator  Spars 

Ill  this  section  we  estimate  the  vertical  loading  using  the  stress  on  a  bending  bridge 
and  an  assumed  lifting  distribution  for  the  stabilator. 

Using  a  bending  moment  to  determine  the  vertical  load  results  in  a  non-unique  problem. 
In  other  words,  an  infinite  number  of  different  loading  distributions  could  produce  the  same 
bending  moment  at  a  particular  location. 

Consider,  for  example,  a  simple  horizontal  cantilever  beam  with  two  vertical  loads. 
Let  us  say  that  these  two  load.  f\  and  fo  (located  LI  m  and  2.5  m,  respectively,  from  the 
fixed  end),  produce  a  7.8  Nm  bending  moment  at  the  root.  In  mathematical  notation,  we 
need  to  solve  for  the  loads  f\  and  f<i  in  the  equation  1.1  f\  +  2.5 f 2  =  7.8.  In  other  words, 
we  are  solving  for  two  unknowns  with  one  equation,  which  results  in  an  infinite  number 
of  solutions.  Even  worse,  changing  the  location  of  these  vertical  loads  would  still  yield 
sensible  solutions;  for  example,  there  exist  pairs  (/i,/2)  such  that  0.3/i  +  I.5/2  —  7.8. 
Distributed  loadings  are  even  more  indeterminate3  than  this  simple  point  load  example. 

Choosing  a  particular  loading  distribution  shape  results  in  a  unique  vertical  load — but 
there  is  little  evidence  to  validate  this  choice,  especially  in  an  airflow  subject  to  disturbance 
from  main  rotor  downwash.  Elliptic  loading  distributions  are  often  used  in  simple  aerofoil 
theory  (see,  for  example,  Houghton  and  Carruthers  [13]),  which  is  why  we  make  the 
following  assumption: 

3 This  greater  indeterminacy  is  more  easily  seen  if  a  distributed  load  is  modelled  as  a  large  number  of 
point  loads.  This  model  results  in  the  equation  Xifi  —  M,  where  n  is  the  number  of  point  loads  fi  at 
locations  Xi ,  and  M  is  the  resulting  bending  moment.  The  degree  of  indeterminacy  is  given  by  the  number 
of  point  loads  n.  Letting  n  tend  to  infinity  and  the  max  \xi  —  Xj\  — ►  0  for  any  1  <  i,j  <  n,  produces  a 
distributed  load — together  with  an  infinite  degree  of  indeterminacy. 
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Assumption  6  The  vertical  loading  acting  on  each  spar  is  elliptic ,  and  is  given  by 
p(x)  =  Py/l  —  (x/L)2,  where  x  is  the  spanwise  location  along  the  spar,  P  is  a  load 
scaling  constant,  and  L  is  the  spar's  length . 


Sikorsky  appears  to  also  use  an  elliptic  loading  distribution  in  their  design  analysis  [9, 
Fig.  7].' 

Both  the  forward  and  aft  spars  begin  at  BL  9  and  end  at  approximately  BL  83.  Thus 
the  spar  length  is  L  =  74  inches  and  x  is  measured  from  the  inboard  edge  of  the  stabilator 
at  BL  9  (that  is,  x  =  0  at  BL  9  and  x  =  74  at  BL  83).  The  unknown  constant  P  not  only 
scales  the  elliptic  loading,  but  is  also  this  loading’s  maximum,  which  occurs  at  the  spar’s 
root. 

One  way  to  improve  Assumption  6  is  by  using  the  bending  stresses  at  all  five  stabilator 
locations.  The  elliptic  lifting  distribution  in  Assumption  6  (or  a  more  complex  distribu¬ 
tion)  could  then  be  scaled  based  on  some  weighted  combination  of  these  five  bending 
stresses.  For  expediency  we  avoided  this  weighted  combination  approach. 

Referring  to  Figure  3.1,  the  bending  moment  at  the  gauge  location  is  given  by 

Mg  =  f  (x  —  xg)p(x)  dx ,  (3.1) 

dx  g 

where  xg  is  the  distance  of  the  bending  bridge  gauge  from  the  spar’s  root.  Substitute  the 
elliptic  loading  distribution  given  by  Assumption  6  into  Equation  (3.1);  then  solving  the 
integral  gives 

A/g  =  P(X  -  U),  (3.2) 

where  the  constants  x  and  W  are  given  by 


2  .  r  2 
'  XG 


)  \[r2  Tg1 


and 


tp  = 
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7f  —  2  tan 


-i 


Xcy 
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As  a  simple  check  of  this  equation,  consider  the  case  where  the  gauge  is  located  at 
the  root  of  the  spar,  that  is,  xg  =  0.  The  vertical  loading  (that  is,  area  of  the  quarter- 
ellipse)  is  it  PL/4  and  the  horizontal  centroid  (of  the  quarter-ellipse)  occurs  at  a  distance 
4L/(37t)  from  the  root.  Multiplying  the  loading  by  the  centroidal  distance  gives  the 
bending  moment  at  the  root  Mt  —  PL2/3,  which  is  the  same  solution  that  Equation  (3.2) 
together  with  Equation  (3.3)  gives. 

Using  simple  beam  bending  theory  (see,  for  example,  Hall  [11]  or  Roark  [28]),  the  fibre 
stress  at  a  point  that  is  a  height  y  above  the  beam’s  neutral  axis  is  given  by 


a  — 


My 
I  ’ 


(3-4) 


where  M  is  the  bending  moment  on  and  I  is  the  second  moment-of-area  of  the  beam’s 
cross-sectional  face  passing  through  the  point.  This  linear  (with  height)  stress  profile  is 
shown  in  Figure  3.1,  along  with  the  resulting  stress  profile  (shown  in  perspective)  on  the 
beam’s  cross-section. 
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Figure  3.2:  Horizontal  (spanwise)  loading,  a  product  of  stress  and  area. 


Use  Equation  (3.4)  to  substitute  the  stress  measured  at  the  bending  bridge  <7g  for  the 
bending  moment  Mg.  Then  for  a  given  bending  moment  Mg  we  can  solve  Equation  (3.2) 
for  the  maximum  load  on  the  spar 


crg  I 

h(x  -  ip)  ’ 


(3.5) 


where  I  is  the  spar’s  second  moment  of  area  and  h  is  the  vertical  distance  from  the  spar’s 
neutral  axis  to  the  skin  (which  is  where  the  bending  bridge  is  located). 

The  forward  spar  was  partitioned  into  three  rectangular  regions  (shown  in  Figure  3.2 
by  the  curly  braces)  in  order  to  calculate  the  second  moment  of  area  of  (and  later  the 
horizontal  load  on)  a  particular  cross-section.  In  the  forward  spar: 


•  the  upper  rectangular  region  contains  the  boron  composite,  the  adhesive  bond,  and 
the  cap’s  flange; 

•  the  middle  rectangular  region  contains  the  cap’s  web  and  part  of  the  spar’s  web;  and 

•  the  lower  rectangular  region  contains  only  the  remainder  of  the  spar’s  web. 

A  similar  partitioning  of  the  aft  spar  was  needed  to  calculate  the  second  moment  of 
area  on  it.  In  the  aft  spar: 

•  the  upper  rectangular  region  contains  only  the  cap’s  flange; 

•  the  middle  rectangular  region  contains  the  cap’s  web  and  part  of  the  spar’s  web;  and 

•  the  lower  rectangular  region  contains  only  the  remainder  of  the  spar’s  web. 


Due  to  the  5°  declination  of  the  cap’s  flange,  the  upper  and  middle  rectangular  regions 
overlap.  Again  due  to  the  declination,  the  maximum  stress  that  the  aft  spar  experiences 
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will  be  greater  than  the  stress  at  the  bending  bridge  gauge.  In  Figure  3.2,  the  gauge 
measuring  stress  is  depicted  as  occurring  0.087"  from  the  top  of  the  spar,  this  length  was 
obtained  by  simple  geometry  2"  x  sin(5°)  -r  2  =  0.087". 

Using  Roark  [28]  the  second  moment  of  area  of  a  rectangular  region  about  its  centroid 
is  given  by  7  =  bd3/ 12,  where  b  and  d  are  the  rectangle’s  width  and  height,  respectively. 
The  second  moment  of  area  for  the  three  rectangular  regions  (defined  in  Figure  3.2)  about 
their  centroids  are  as  follows. 

For  the  forward  spar 

if,  =  2.48 x  10~3  in4, 

I{2  =  6.46x1 0“3  in4,  and 
/f3  =  0.032 (/if  -  1.132)3/12  in4, 

where  /if  and  /ia  respectively  denote  half  the  height  of  the  forward  and  aft  spars,  see 
Figure  2.3  on  page  12.  (Thus  the  values  shown  in  Table  2.1  on  page  13  represent  2/if  and 
2ha-) 

For  the  aft  spar 

7al  =  95Gx  10~6  in4, 

7a2  =  13.1  xlO-3  in4,  and 
7a3  -  0.032(/i.a  -  0.125)3/12  in4, 

where  we  have  used  the  formula  Iq  —  bd(d, 2  cos2  9  +  b2  sin2  9)/ 12  (see  Roark  [28])  to  derive 
the  second  moment  of  area  for  rectangular  Region  1  of  the  aft  spar.  This  formula  was 
required  because  Region  1  was  aligned  at  a  6  =  5°  declination  to  the  horizontal  axis. 

The  formula  Ix  =  7/  +  Ay2  (see  Roark  [28])  relates  the  second  moment  of  area  about 
two  axes  separated  by  a  perpendicular  distance  y.  Using  this  formula,  summing  the  three 
different  regions,  and  doubling  the  result  (to  account  for  the  lower  half  of  the  spars)  gives 
the  second  moment  of  area  for  the  forward  spar 

7f  =  2^  [2.48 x  10"3  +  0.508(/tf  -  0.121)2] 

+  [6.46  x  10-3  +  0.098(/)f  -  0.687)2] 

+  jo.032(/if  -  1.132)3/12  +  0.032(/if 
=  0.0213/if3  +  1.14/if2  -  0.433/i.f  +  0.0942 
and  for  the  aft  spar 

7a  =  2^  [956 x  10-6  +  0.250 {ha  -  0.063)2] 

+  [13.1  x  10~3  +  0.156(/ia  -  0.625)2] 

+  |o.032(/i.a  -  1.125)3/12  +  0.032 (^ 

=  0.0213ha3  +  0.742/i.a2  -  0.374/ia  +  0.122. 


-  1.132)  [(/if  —  1.132)/2] 2  j  ^ 

(3.6) 


-  1.125)  [(/ia-1.125)/2]2}) 

(3.7) 
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The  vertical  load  at  the  spar's  root  is  the  integral  of  the  elliptic  loading,  that  is, 
V  =  7rPL/4.  (This  loading  is  the  same  as  the  area  of  a  quarter-ellipse  with  semi-axes 
lengths  given  by  the  spar’s  maximum  loading  and  length.)  From  Equations  (3.3)-(3.7)  we 
have  that  the  vertical  load  at  the  forward  and  aft  spars’  roots  are  respectively 


Vf=  -- 


(0.0213/if3  +  1.14/if2  -  0.433/if  +  0.0942) (37r<rgL2//if) 


2  (2L2  +  xg2)  y/L2  -  xg2  -  3 xgL2 


2  tan  1  (xg/v^^V)] 


(3.8) 


and 


Fa  = 


(0.0213/ia3  +  0.742 ha2  -  0.374/?,a  +  0.122) (37ragL2/ha 


2  (2 L2  +  xg2)  yj L2  —  xg2  —  3xgL2  tx  —  2  tan  1  (^xg/y/L2  —  xg2^j 


(3.9) 


A  negative  stabilator  bending-stress  implies  tension  in  the  cracked  lug  (see  Assumption  2), 
and  hence  Equations  (3.8)  and  (3.9)  have  negative  signs  in  front  of  them. 


3.3  Horizontal  (Spanwise)  Load  on  Stabilator  Spars 

In  this  section,  we  calculate  the  horizontal  (spanwise)  stress  on  the  spar’s  cross-section 
using  the  bending  stress  measured  on  the  skin. 

The  horizontal  loading  at  a  particular  cross-section  is  calculated  by  multiplying  the 
horizontal  stress  by  the  cross-sectional  area.  Figure  3.2  (on  page  17)  shows  how  the  hori¬ 
zontal  load  was  calculated  for  the  forward  and  aft  spars. 

Some  simple  but  time  consuming  calculations  give  the  horizontal  loads  for  the  three 
rectangular  regions  of  both  the  forward  and  aft  spars.  First,  we  need  to  determine  the 
cross-sectional  area  of  each  rectangular  region.  Second,  we  need  to  determine  the  stress 
at  the  highest  and  lowest  parts  of  these  rectangles,  which  will  involve  the  variable  height 
of  the  spar’s  web. 

The  areas  and  maximum  stresses  for  the  rectangular  Regions  1,  2,  and  3  of  the’  forward 
spar  are  respectively: 


An  =  0.508, 

max(<7fi)  =  ex  f, 

(3.10) 

An  =  0.098, 

max(<7f2)  =  crf(l 

-  0.242/ftf), 

(3.11) 

Aa  =  0.032 (%  -  1.132), 

and  max(<Xf3)  =  <7f(l 

-  1.132/%), 

(3.12) 

where  An  and  max(<jfi)  respectively  denote  the  area  and  maximum  stress  of  the  ith  rect¬ 
angular  region  of  the  forward  spar.  The  symbols  <jf  and  cra  denote  the  stresses  measured 
by  the  bending  bridges  on  the  forward  and  aft  spars,  respectively  (see  the  tops  of  the  stress 
profiles  in  Figure  3.2). 

The  areas  and  maximum  stresses  for  the  rectangular  Regions  1,2,  and  3  of  the  aft  spar 
are  respectively: 


An  =  0.250, 

max(orai)  =  cra(l  +  0.087/%), 

(3.13) 

Aa2  =  0.157, 

max(era2)  —  <ra(l  —  0.125/%), 

(3.14) 

Aa3  =  0.032(%  -  1.125), 

and  max(o-a3)  =  cra(l  -  1.125/%), 

(3.15) 
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where  Aixi  and  max(aa#)  respectively  denote  the  area  and  maximum  stress  of  the  ith 
rectangular  region  of  the  aft  spar. 

The  minimum  stress  in  a  rectangular  region  is  the  same  as  the  maximum  stress  in  the 
region  immediately  below  it.  Mathematically,  the  minimum  stress  of  the  three  rectangular 
regions  belonging  to  the  aft  spar  are 


min(<raj) 


max  [<ra(i+1)]  if  i  =  1, 2, 
0  if  i  =  3. 


(3.16) 


An  analogous  equation  arises  for  the  minimum  stresses  of  the  forward  spar.  The  minimum 
stress  of  rectangular  Region  3  (in  both  spars)  is  zero  because  this  point  is  the  centre  line 
of  symmetry,  and  therefore  the  neutral  axis. 

Within  any  of  these  rectangular  regions,  the  stress  varies  linearly  between  the  region’s 
maximum  and  minimum  stress.  All  of  the  rectangular  regions  in  the  forward  and  aft 
spars  are  either  symmetric  or  antisymmetric  about  the  horizontal  line  passing  through 
each  rectangle’s  vertical  midpoint.  In  other  words,  the  centroid  coincides  with  the  vertical 
midpoint.  Thus  the  equation  for  the  horizontal  load  in  the  ith  region  is  simply  given  by 
the  stress  in  each  region’s  midpoint, 


Hi  =  Ai 


max((7j)  +  min(<7j) 
2 


(3.17) 


where  A{  is  the  area  of  the  ?th  region,  and  max  (o',)  and  min(<7j)  are  the  maximum  and 
minimum  stress  on  the  ith  region. 

Using  Equations  (3.10)-(3.12)  and  Equation  (3.17)  we  obtain  the  horizontal  loads  on 
the  three  rectangular  regions  of  the  forward  spar. 

Hn  =  —0.508<7f(l  -  0.121/7if),  (3.18) 

Ha  =  — 0.098<7f(l  -  0.687//if),  and  (3.19) 

Ha  =  — 0.032<7f(/if  -  1.132)(0.5  -  0.566/%).  (3.20) 


Analogously,  using  Equations  (3. 13)— (3. 1 5)  and  Equation  (3.17)  we  obtain  the  horizontal 
loads  on  the  three  rectangular  regions  of  the  aft.  spar. 

Ha\  =  -0.250era(l  -  0.019//ta),  (3-21) 

Hh2  =  — 0.157(Ta(l  -  0.625/%),  and  (3.22) 

H&  =  — 0.032<ra(%  -  1.125)(0.5  -  0.563/%).  (3.23) 

A  negative  stabilator  bending-stress  implies  tension  in  the  cracked  lug  (see  Assumption  2), 
and  hence  Equations  (3.18) -(3.23)  have  negative  signs  in  front  of  them. 

The  horizontal  loads  given  by  Equations  (3.18)-(3.23)  under-estimate  the  stress  on  the 
cracked  lug  for  the  following  reason:  the  portion  of  distributed  loading  between  the  spar’s 
root  and  the  bending  bridge  (used  to  estimate  the  horizontal  loading)  were  excluded.  This 
loading  abutting  the  root  is  excluded  by  the  fact  that  the  bending  stress  is  measured  at 
some  point  away  from  the  stabilator’s  root. 
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In  order  to  correct  this  under-estimation  the  loading  profile  given  by  Assumption  6  is 
needed.  We  have  already  seen  that  the  bending  moment  at  the  root,  due  to  an  elliptic 
loading  profile,  is 


Mr  —  P  [  xy/l  -  ( x/L)2dx 

Jo 


PL 2 

~y~ 


If  the  vertical  loading  between  the  gauge  and  the  root  is  ignored  the  resulting  bending 
moment  at  the  root  is 


fL 

M*  =  P  x\J\  —  (x/L)'2dx 

Jx g 


PL2 

3 


remember  that  xg  is  the  distance  from  the  spar’s  root  to  the  gauge.  Taking  the  ratio  of 
these  two  bending  moments  gives  the  correction  factor  for  the  horizontal  loading  as 


7  -  mv/m;  = 


The  total  horizontal  (spanwise)  load  on  the  top  half  of  each  spar  is  given  by  the  sum 
of  the  three  forces  on  each  spar.  Thus  approximations  to  the  total  horizontal  load  on  each 
spar’s  top  half  are  given  by 

Hf  =  -O.OI6O7  af(fif-0.189)(fif +  35.8)/fif  (3.24) 

and 

Ha  =  —0.0160 7  <7a  (ha  -  0.221)(fia  4-  23A)/ha  (3.25) 

respectively  for  forward  and  aft  spars. 

As  is  shown  next,  the  horizontal  loading  can  also  be  calculated  using  information  from 
the  vertical  loading.  When  we  later  calculate  stress  on  the  cracked  lug  (in  Section  6.2), 
we  will  choose  the  higher  of  these  two  values  for  horizontal  loading. 


3.4  Using  Vertical  Load  to  Calculate  Horizontal  Load  on 
Stabilator  Spars 

We  can  use  the  estimate  of  vertical  loading  given  by  Equations  (3.8)  and  (3.9)  to 
calculate  the  horizontal  (spanwise)  load  on  the  spar’s  root.  Below  are  the  calculations 
used  to  obtain  this  horizontal  load  estimate. 

It  is  possible  to  determine  the  bending  moment  (and  hence  horizontal  load)  at  the 
spar’s  root  using  the  estimate  of  vertical  load.  Multiplying  the  vertical  load  V  and  the 
centroidal  distance  of  the  loading  from  the  root  (namely,  4L/(37r))  gives  the  bending 
moment  at  the  root 

4  L 
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where  L  is  the  length  of  the  spar.  In  the  above  bending  moment  equation  we  have  implicitly 
used  the  elliptic  loading  profile  given  by  Assumption  6. 

It  is  clear-  from  Figure  3.1  that  any  bending  moment  at  the  spar’s  root  is  reacted  out 
as  a  couple  by  the  spar’s  supporting  lugs.  Let  }j(  be  the  distance  between  the  upper  and 
lower  supports,  that  is,  upper  and  lower  clevis  lugs  (see  Figure  3.1).  From  Figures  A4  and 
A12  (on  pages  55  and  63,  respectively)  the  distance  between  the  upper  lug  and  lower  lug 
is  ye  =  5".  Summing  moments  about  the  lower  support,  we  obtain  the  horizontal  load  on 
the  upper  support  as 


H  =  Mr/yt 

=  JL  V 

Smjc 

=  6.28F, 


(3.26) 


where  V  is  given  by  either  Equation  (3.8)  or  Equation  (3.9)  respectively  for  either  the 
forward  or  aft  spar. 
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4  Information  from  Flight  Loads  Survey 

All  the  loading  information  used  in  this  report  was  derived  from  the  measurements 
taken  during  the  Flight  Loads  Survey.  In  order  to  make  some  approximations  about  the 
“typical”  loading  experienced  by  the  cracked  lug  we  must  first  aggregate  the  Flight  Loads 
Survey  data  into  loading  groups.  Instead  of  using  the  raw  data  for  the  fatigue  analysis, 
the  Survey’s  statistically  reduced  data  set  was  used.  This  reduced  data  set  is  explained  in 
the  second  part  of  this  section. 


4.1  Partitioning  526  Manoeuvres  into  12  Groups 

All  up,  there  were  526  individually  named  manoeuvres  (and  approximately  3800  runs) 
in  the  Flight  Loads  Survey.  To  gain  some  appreciation  of  whether  particular  manoeu¬ 
vres  induced  high  stress  in  the  stabilators,  these  526  manoeuvres  were  partitioned  into 
12  groups.  And  so  we  make  the  following  assumption: 

Assumption  7  All  manoeuvres  within  a  particular  group  result  in  comparable  bend¬ 
ing  stresses  on  the  two  stabilator  spars. 

This  assumption  merely  states  that  we  can  sensibly  group  manoeuvres  together,  which  is 
reasonable  as  far  as  determining  general  trends  goes. 

The  526  manoeuvres  were  initially  partitioned  according  to  the  work  of  King  and 
Lombardo  [15,  pp.  62-4],  who  list  21  groups.  These  groups  were  then  further  aggregated 
into  12  groups.  (The  individual  manoeuvres  in  each  group  are  listed  in  Tables  B1-B12  on 
pages  69-74.)  This  modification  of  the  King  and  Lombardo  (K-L)  grouping  was  carried 
to  (i)  have  approximately  a  dozen  groups  (instead  of  the  21  groups  K-L  formulated)  and 
(ii)  obtain  sensible  stress  plots  (see  Section  5.3). 

As  can  be  seen  from  Table  4.1,  which  compares  the  groupings  of  this  report  with  K-L, 
approximately  half  of  the  groups  K-L  used  were  absorbed  into  the  miscellaneous  group 
in  this  report.  On  the  other  hand,  the  pullouts  group  of  K-L  was  partitioned  into  two 
groups  in  this  report.  The  K-L  groupings  were  modified  for  one  of  two  reasons,  either  a 
group  contained  too  few  or  too  many  manoeuvres.  The  groups  that  contained  too  few 
manoeuvres  were  absorbed  into  the  miscellaneous  group.  In  contrast,  the  pullouts  group 
from  the  K-L  report  contained  too  many  manoeuvres,  and  so  it  was  partitioned  into  two 
groups:  symmetric  pullouts  and  rolling  pullouts. 

A  few  of  the  groups  formulated  by  King  and  Lombardo  could  not  be  mapped  onto  the 
set  of  manoeuvres  from  the  Flight  Loads  Survey  for  one  of  two  reasons: 

•  The  group  contained  no  manoeuvres  from  the  Flight  Loads  Survey  (for  example,  the 
droop-stop  pounding  manoeuvre). 

•  The  group  was  not  explicitly  named  in  the  Flight  Loads  Survey  (for  example,  the 
ground- air-ground  cycles  were  decomposed  into  their  intermediate  manoeuvres). 


Manoeuvres  were  partitioned  into  the  12  groups  with  a  clear  goal  in  mind—to  estimate 
the  stress  on  the  cracked  stabilator  lug.  As  such,  caution  should  be  exercised  when  using 
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Talkie  4.1:  Comparing  the  manoeuvre  groupings  used  by  King  and  Lom¬ 
bardo  (K-L)  and  This  report.  The  symbol  &  denotes  approximately  the 
sam,e  group,  $  a  non-existent  group ,  and  0  a  group  with  no  manoeuvres. 
The  symbols  ©n  and  Qn  denote  a  group  containing  n  manoeuvres  that 
was  absorbed  into  the  miscellaneous  or  climb  group ,  respectively. 


K-L  report 

This  report 

autos 

break  turn 

0 

climb 

dive 

©2 

droop-stop  pounding 

0 

ground  run 

0 

ground-air-ground  cycles 

$ 

heavy  manoeuvres 

0 

hover 

©15 

landing 
level  flight 

©2 

$ 

miscellaneous 

nap  of  earth 

©2 

partial  power  descent 

05 

pullouts 

J  symmetric  pullouts 
\  rolling  pullouts 

reversals 

rotor  start/stop 
side/rearward  flight 

©9 

O' 

sideslip 

take  off 

ril 

®5 

taxi 

©IS 

turns 

these  manoeuvre  groups  for  any  other  purpose.  Appendix  B  lists  the  twelve  groups  into 
which  the  manoeuvres  from  the  Flight  Loads  Survey  were  partitioned. 
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4.2  Statistically  Reduced  Data  Set  of  Flight  Loads  Survey 

Instead  of  using  the  time-history  data  from  the  Flight  Loads  Survey,  a  statistically 
reduced  data  set4  of  the  Survey  was  used.  This  reduced  data  contained  ten  parameters  (see 
Table  4.2)  that  described  the  time-history  data  set.  We  denote  the  steady  and  vibratory 
components  of  the  stress  measurements  respectively  by  the  symbols  “<r”  and  CV”,  and 
the  direct  stress  measurements  by  the  symbol  “s”.  On  any  of  these  three  symbols,  the 
subscripts  “min”,  “avg”,  “max”,  “5%”,  and  “95%”  respectively  denote  the  minimum, 
average,  maximum,  and  5th  and  95th  percentile  stress. 


Table  4.2:  The  ten  statistical  parameters  that  describe  the  time-history 
data  set  for  the  Flight  Loads  Survey .  The  maximum  and.  minimum, 
stresses  are  not  shown  in  this  section's  stress  plots,  and  hence  were  not 
assigned  symbols .  The  mnemonics  shown  were  assigned,  by  GTRI. 


Symbol 

Mnemonic 

Parameter  description 

1 

0min 

MIN.STDY 

minimum  of  steady  stress  component 

steady  < 

<7aVg 

AVGJ3TDY 

average  of  steady  stress  component 

1 

<  &  max 

MAX.STDY 

maximum  of  steady  stress  component 

1 

^avg 

AVG.VIB 

average  of  vibratory  stress  component 

vibratory  < 

^max 

MAX.VIB 

maximum  of  vibratory  stress  component 

1 

PCL.VIB 

95th  percentile  of  vibratory  stress  component 

1 

ABMIN_STDY 

minimum  directly  measured"  stress 

direct  < 

- 

ABMAX_STDY 

maximum  directly  measured  stress 

S5% 

PCCJ)IR 

5th  percentile  of  directly  measured  stress 

| 

k  S  95% 

PCLJDIR 

95th  percentile  of  directly  measured  stress 

According  to  King  and  Boykett  [14],  one  way  to  calculate  the  steady  and  vibratory 
components  of  stress  from  the  time-history  data  is  as  follows: 

•  Determine  all  stress  cycles ,  which  are  defined  as  two  consecutive  turning  points  (Si 
and  52)  from  the  time-history  data. 

•  The  steady  component  of  the  stress  is  defined  as  the  mean  of  these  two  turning 
points,  that  is,  a  =  (Si  +  S2V2. 

•  The  vibratory  component  of  the  stress  is  defined  as  half  the  difference  (in  modulus) 
of  these  two  turning  points,  that  is,  v  =  \S\  —  £2 1/2. 

In  contrast,  before  performing  the  steps  described  above  Sikorsky5  pre-process  the 
data  for  each  flight  run  to  reduce  the  number  of  turning  points.  This  pre-processing 

4The  time-history  data  set  was  approximately  300  GB  in  size  and  was  stored  on  25  CDs.  In  contrast, 
the  statistically  reduced  data  set  was  approximately  0.1  GB  in  size  and  fits  on  a  single  CD. 

5Sikorsky  calculated  these  statistical  parameters  for  a  GTRI  report. 
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determined  the  peaks  and  troughs  inside  fixed  time  windows.  Sikorsky  chose  the  tail  rotor 
frequency  (which  is  19.8  Hz)  as  the  most  important  loading  frequency,  and  hence  used  time 
windows  of  0.051  seconds  (which  is  the  inverse  19.8  Hz).  Compared  to  the  raw  sampling 
rate  of  833  Hz  used  during  the  Flight  Loads  Survey,  Sikorsky’s  pre-processing  results  in  a 
reduction  in  the  number  of  turning  points.  This  turning  point  reduction  yields  larger  (and 
hence  more  conservative)  vibratory  stresses,  which  are  used  to  derive  fatigue  damage. 

Like  King  and  Boykett  [14],  we  have  not  used  Sikorsky’s  data  pre-processing  technique 
to  determine  the  GTRI  statistical  parameters. 

Table  4.3  compares  the  results  obtained  by  GTRI  with  the  results  obtained  using  the 
more  accurate  procedure  outlined  above.  If  “s”  is  a  direct,  vibratory,  or  steady  stress, 
then  the  percentage  difference  shown  in  this  table  was  calculated  as  follows: 

%-diff  =  (sG  —  sR)/  sG, 

where  sG  and  sH  are  the  stresses  calculated  respectively  by  GTRI  and  this  report 


Table  4.3:  Percentage  difference  of  statistical  parameters  obtained  by 
GTRI  and  this  report  (which  are  labelled  Our).  All  units  are  in  psi, 
except  for  the  error  which  is  percentage  error.  The  stress  on  bulkhead 
FS308  was  measured  during  a  45°  left  turn.  The  stress  on  the  stabilator 
was  measured  during  level  flight. 


Bulkhead  (psi)  Stabilator  (psi) 


Parameter 

GTRI 

Our 

Diff. 

GTRI 

Our 

Diff. 

1 

f  O'min 

-8237 

-8832 

-7.2% 

-18102 

-24016 

-33% 

steady  < 

1 

|  °kvg 

-7903 

-7844 

0.75% 

-16644 

-16520 

0.75% 

1 

l  ^max 

-7460 

-7133 

4.4% 

-15120 

-9062 

40% 

i 

j'  J'avg 

621 

265 

57%, 

6664 

1636 

75% 

vibratory  < 

1  ^max 

1104 

1064 

3.6% 

9569 

8799 

8.0% 

1 

l  ^95% 

880 

565 

36% 

8527 

4808 

44% 

'  ABMIN_STDY 

-8988 

-8988 

0% 

-26853 

-26853 

0% 

ABMAX_STDY 

-6669 

-6669 

0% 

-7042 

-7042 

0% 

direct < 

sr>% 

-7311 

-7324 

-0.18% 

-11594 

-11755 

-1.4% 

.  S95f% 

-8483 

-8474 

0.11% 

-21640 

-21564 

0.35% 

The  statistical  parameters  were  compared  with  the  GTRI  results  for  two  manoeuvres: 

•  a  45°  left  turn  from  run  46  of  flight  28  (REC0VRYuLTuTN ,  UVH ,  U45DEG)  and 

#  level  flight  from  run  82  of  flight  24  (LEV .  FLTuF0RWARD ,  u0 . 9VH) . 
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The  left  turn  manoeuvre  was  chosen  to  compare  results  with  King  and  Boykett  [14].  For 
this  left  turn  manoeuvre,  measurements  were  taken  from  a  strain  gauge  on  the  FS308 
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bulkhead6  (see  the  description  given  by  Friend  [8,  p.  44]  for  gauge  BEAMAS1).  The  level 
flight  manoeuvre  was  chosen  randomly  as  an  additional  check.  For  this  level  flight  ma¬ 
noeuvre,  measurements  were  taken  from  the  right-hand  side  bending  bridge  Number  1  (see 
Figure  2.1  and  the  description  for  gauge  STBNBM1R  in  Section  2.1). 

As  can  be  seen  from  Table  4.3,  the  two  results  are  identical  for  the  absolute  maximum 
and  minimum  of  the  direct  stress.  The  5th  and  95th  percentile  are  within  1%  error,  as  is 
the  average  steady  component  (crav g).  The  remaining  parameters  are  significantly  different. 
These  large  differences  are  explained  by  the  definition  of  a  turning  point. 

We  have  used  the  mathematical  definition  of  a  turning  point,  that  is,  either  a  local 
maximum  or  minimum.  Using  this  definition  of  a  turning  point,  the  magnitudes  of  the 
steady  and  vibratory  components  are  respectively  over-  and  under-estimated  as  compared 
to  the  GTRI  magnitudes  (see  Table  4.3).  This  over-  and  under-estimation  show  that 
GTRI  ignored  certain  intermediate  turning  points  that  they  considered  insignificant  for 
the  fatigue  analysis.  As  an  example,  in  Figure  5.2  (on  page  33)  the  red  points  (which 
denote  turning  points)  were  determined  using  the  mathematical  definition  of  a  turning 
point.  The  arrow  in  this  figure  marks  a  pair  of  turning  points  GTRI  may  be  ignoring 
when  calculating  the  statistical  parameters. 


6On  the  lower  face  of  the  upper  flange  of  the  BL34.5  longeron,  just  aft  of  FS308  on  the  left-hand  side. 
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5  Bending  Stress  on  Stabilator  Spar 

In  this  section  we  determine  the  stress  experienced  by  the  stabilator  spars  in  three 
different  ways.  The  first  method  assumes  the  worst-case  scenario  and  determines  the  spar 
bending  stress  by  considering  the  combination  of  maximum  stress  and  alternating  stress 
that  yield  the  highest  fatigue.  As  an  alternative,  the  second  method  considers  a  typical 
high-loading  case  and  determines  the  bending  stress  that  more  closely  approximates  the 
cracked  lug’s  fatigue.  Unlike  the  worst-case  scenario,  however,  this  high-loading  case  is  not 
necessarily  conservative.  The  final  method  assumes  the  best-case  scenario  and  determines 
the  spar  bending  stress  by  considering  the  median  loading  for  both  the  steady  and  vibratory 
components  of  stress. 


5.1  Fatigue  Equation  for  7075  Aluminium  Alloy 

Before  we  can  determine  the  worst-case  stress  combination,  we  need  to  know  the  fatigue 
properties  of  the  cracked  lug’s  material,  which  is  7075-T7  aluminium  alloy  [6], 

According  to  Higdon  et  al.  [12,  p.  573],  the  stress  concentration  factor  to  use  in  the 
fatigue  analysis  of  a  lug  with  the  geometry  of  the  cracked  lug  is  approximately  2  (see,  in 
particular,  Figure  ll-9(c)  in  Higdon  et  al.  [12]).  Thus  the  correct  fatigue  curve  to  use  in 
the  fatigue  analysis  is  an  S-N  curve  for  notched  specimens  with  Kt  =  2. 

Let  Smax  and  Smin  respectively  denote  the  maximum  and  minimum  stresses  of  the 
loading  cycles  (given  in  ksi)  and 

R  =  Smm/Sma*  '  '  (5.1) 

denote  the  stress  ratio.  Using  MIL-HDBK-5  [7,  p.  3-411],  which  is  the  abbreviation  for 
Volume  5  of  the  Military  Handbook  for  Metallic  Materials,  the  equation'  for  cycles  to 
failure  of  Kt  —  2  notched  7075-T6  aluminium  alloy  is  given  by 

Nt  =  31.6  x  106  (Seq  -  18.6)-2’46  ,  (5.2) 

where  Seq  is  the  equivalent  stress  defined  by 

=  Smax(l  ~  R)^ •  (5-3) 

Equations  (5.1)  and  (5.3)  were  developed  using  net  stresses.  See  the  MIL-HDBK-5  for 
restrictions  on  using  Equation  (5.2). 

Since  MIL-HDBK-5  does  not  contain  fatigue  information  on  7075-T7  aluminium  alloy 
(the  material  of  the  cracked  lug)  we  make  the  following  assumption: 

<^>  Assumption  8  The  cycles  to  failure  formula  given  by  Equation  (5.2)  for  7075-T6 

-n- -  aluminium,  alloy  is  also  valid  for  7075-T7  aluminium  alloy. 

As  a  comparison  of  these  two  alloys,  the  ultimate  and  yield  tensile  strengths  for  these  two 
temperings  are: 

7ln  MIL-HDBK-5,  Equation  (5.2)  is  shown  in  logarithmic  (base  10)  form. 
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•  jPtu  —  74  ksi  and  =  63  ksi  for  the  T6  tempering  treatment  [7,  p.  3-371]  and 

•  Ftu  —  67  ksi  and  Fty  =  56  ksi  for  the  T7  tempering  treatment  [7,  p.  3-373]. 

Thus  we  might  expect  that  Equation  (5.2)  results  in  a  slight  under-estimation  of  cycles  for 
the  T7  tempering.  This  expectation  is  supported  by  experimental  data  [19,  Fig.  3.0552] 
comparing  T6  and  T7  tempering.  Hence  Assumption  8  appears  to  be  conservative. 

The  cycles  to  failure  formula  is  based  on  experimental  data  with  bounded  stress  ratios, 
and  thus  the  following  assumption  is  required: 

Assumption  9  The  cracked  lug  is  not  pre-stressed  in  tension  and  the  cycles  to  failure 
formula  given  by  Equation  (5.2)  can  be  extrapolated  into  the  ranges  R  <  —  1  and 
0.5  <  R<  1. 


Assuming  that  the  stabilator  lugs  are  not  pre-stressed  in  tension,  then  the  purely  com¬ 
pressive  load  cycles  ( Smax  <  0  and  R  >  1)  can  be  ignored  in  the  fatigue  analysis.  Hence 
Assumption  9  is  not  affected  by  the  range  R  >  1.  We  are  considering  the  worst-case  sce¬ 
nario,  and  hence  assuming  that  the  most  fatiguing  combination  of  5max  and  R  is  acting 
during  all  manoeuvres.  Provide  we  can  extrapolate  the  cycle-to-failure  contours,  then 
from  the  zoomed  section  in  Figure  5.1  (on  page  30)  it  is  clear  that  the  highest  equivalent 
stress  occurs  outside  the  ranges  R  <  —  1  and  0.5  <  R  <  1. 


5.2  Worst-Case  Bending  Stress  from  Flight  Loads  Survey 

One  way  to  determine  the  loading  for  the  worst-case  loading  scenario  is  to  consider  the 
maximum  equivalent  stress,  which  is  given  by  Equation  (5.3).  Plotting  Smax  versus  R  and 
then  superimposing  contours  of  constant  equivalent  stress  on  this  plot  we  can  determine 
the  worst-case  loading  combination  of  Smax  and  R. 

A  log-like8  plot  of  the  maximum  stress  Smax  versus  the  stress  ratio  R  is  shown  in 
Figure  5.1,  where  we  have  taken  iSmin  =  s5%  and  5max  =  s95%.  In  other  words,  for  each  run 
from  the  Flight  Loads  Survey,  we  have  taken  the  minimum  and  maximum  stress  to  be  the 
5th  and  95th  percentile  of  the  directly  measured  stress  for  that  run.  These  3765  runs  were 
partitioned  into  the  twelve  groups  described  in  Section  4.1.  The  short  horizontal  lines  (on 
the  plot’s  right)  show  the  Smax  distribution  of  each  group. 

A  negative  bending  stress  in  the  stabilator  spars  generates  a  tensile  stress  in  the  cracked 
lug  (see  Assumption  2  on  page  8).  Thus  to  obtain  the  correct  combination  of  SmstX  and  R 
for  the  fatigue  analysis  of  the  cracked  lug,  the  sign  of  the  bending  stresses  were  reversed 
in  Figure  5.1. 

The  zoomed  section  in  Figure  5.1  shows  contours  of  cycles- to- failure,  which  were  calcu¬ 
lated  from  Equation  (5.2).  This  cycles-to-failure  equation  is  derived  from  coupons  tested 
at  stress  ratios  of  —1  <  R  <  0.5,  and  hence  the  contours  in  this  zoomed  section  were 
restricted  to  this  same  range  of  stress  ratios. 

8The  horizontal  axis  of  Figure  5.1  is  on  an  arcsinh  scale,  which  is  approximately  linear  about  zero  and 
logarithmic  elsewhere  (as  illustrated  by  the  tickmarks). 
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Figure  5.1:  Plot  of  maximum  stress  (at  bending  bridge  Number  1 )  versus  stress  ratio  for  all  runs  of  the  Fligh 
Loads  Survey.  The  “minimum”  and  “ maximum ”  stresses  were  taken  as  s5%  and  sy5%  (that  is,  5th  and  95tl 
percentile  of  the  directly  measured  stress).  The  short  horizontal  lines  to  the  right  of  the  plot  show  the  Smcu 
distribution  for  each  group.  Contours  of  equivalent  stress  Seq  are  superimposed  onto  the  zoomed  inset. 


Autos  max  (Seq) 
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For  expediency,  we  chose  s5%  and  s95%  to  calculate  fatigue,  and  so  we  require  the 
following  assumption: 

Assumption  10  The  5th  and  95th  percentile  of  stabilator  stress  measured  during  the 
Flight  Loads  Survey  characterise  the  fatigue  loading  that  the  stabilator  lug  experiences. 

For  a  quick-and-dirty  analysis,  we  require  a  single  loading  amplitude  to  calculate  fa¬ 
tigue.  Alternative  ways  to  determine  the  single  loading  amplitude  (listed  in  order  of 
increasing  fatigue  severity)  include: 

•  Using  the  average  steady  stress  combined  with  the  average  vibratory  stress;  that  is, 

Smm  =  max(<7avg.)  -  max(i'avgj)  and  Sniax  =  max(cravgi)  +  max(i/avgA 

l  l  l  % 

where  i  is  the  run  number  and  the  maximum  is  taken  over  all  runs. 

•  Using  the  maximum  and  minimum,  steady  stress  combined  with  the  maximum  vibra¬ 
tory  stress;  that  is, 

S'min  =  max(<7mini)  -  max(i/niaxj)  and  5max  =  niax(crmaxj)  +  max(Vniaxi). 

•  Using  the  maximum,  and  minimum,  direct  stress  measured;  that  is, 

S'min  =  max(s0%j)  and  Smax  —  niax(s10o%^). 

Intuitively,  we  thought  the  stress  given  by  the  first  alternative  would  yield  results  that 
were  unconservative;  while  the  stresses  given  by  the  remaining  two  alternatives  would 
yield  results  that  were  too  conservative. 

The  single  load  amplitude  resulting  from  Assumption  10  yields  a  coarse  approximation, 
which  was  required  to  produce  timely  results.  Without  this  assumption  we  would  need 
to  carefully  analyse  the  usage  spectrum  of  the  helicopter  in  order  to  say  approximately 
how  it  is  used.  Assumption  10  can  be  refined  by  assumptions  that  are  successively  better 
models  of  reality.  Two  examples  of  more  refined  assumptions  are: 

•  Determine  approximately  what  percentage  of  time  the  helicopter  spends  in  each 
manoeuvre  group  and  then  calculate  approximate  stresses  for  each  of  these  groups. 

•  Or  even  more  realistically,  determine  approximately  what  percentage  of  time  the  heli¬ 
copter  spends  in  each  individual  manoeuvre  and  then  calculate  approximate  stresses 
for  each  of  these  manoeuvres. 

Both  of  these  refined  assumptions  head  towards  the  goal  of  recalculating  the  fatigue  life 
of  the  cracked  lug  based  on  data  from  the  Flight  Loads  Survey.  It  must  be  remembered, 
however,  that  the  stress  at  the  lug  is  approximated  using  the  bending  stress  on  the  stabi- 
lator’s  spar.  Hence  it  would  be  questionable  to  aim  for  the  higher  accuracy  suggested  by 
the  recommended  refinements  listed  above  when  the  other  aspects  of  the  fatigue  analysis 
were  of  a  lower  accuracy. 
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From  Equation  (5.3)  and  Figure  5.1  the  most  fatiguing  run  was  a  symmetric  pushover 
manoeuvre,  which  had  a  maximum  stress  of  Smax  =  —30.3  ksi  and  a  stress  ratio  of  R  = 
-0.186.  Substituting  these  values  into  Equation  (5.3)  gives  the  equivalent  stress  as 

Se q  =  -33.2  ksi  (worst-case  loading),  (5.4) 

which  comes  from  run  9  of  flight  75  (with  run  mnemonic  SYMMuPUSH0VER)UVH,u.25G). 
This  maximum  equivalent  stress  is  singled-out  by  the  arrow  in  Figure  5.1.  Remember 
that  the  sign  of  the  stresses  in  Figure  5.1  were  reversed  to  obtain  the  correct  stress  in  the 
cracked  lug.  This  sign  reversal  explains  the  negative  sign  of  the  equivalent  stress  value  in 
Equation  (5.4). 


5.3  High  Bending  Stress  from  Flight  Loads  Survey 

As  will  be  seen  later,  the  fatigue  life  predicted  by  the  worst-case  scenario  is  very  short. 
To  gain  some  appreciation  for  how  over-conservative  this  worst-case  is,  a  typical  fatigue 
life  (based  on  high  loading)  is  now  calculated.  The  difference  between  this  section  and  the 
previous  section  is  that  we  are  now  using  a  high  loading  (based  on  a  typical  level  flight 
run)  and  not  the  worst- case  loading  (determined  over  all  runs).  The  strong  caveat  is  that 
stresses  obtained  in  this  section  are  not  necessarily  conservative,  and  hence  should  only 
be  used  in  making  comparisons  with  the  worst-case  scenario. 

In  this  section  we  determine  a  “typically  high’5  stress  on  the  stabilator  at  a  location 
that  is  9"  outboard  of  the  cracked  lug.  This  typically  high  stress  is  determined  using 
information  from  the  Flight  Loads  Survey  and  engineering  judgements  as  to  what  consti¬ 
tutes  typically  high.  (More  precise  definitions  of  “typically  high”  are  given  later.)  In  this 
section  these  engineering  judgements  are  typeset  as  assumptions  and  explained  when  they 
are  introduced. 

Figure  5.2  shows  a  legend  for  the  stress  plot  in  this  section  (and  the  stress  plots  in 
Appendix  D).  The  time-history  data  in  this  plot  are  taken  from  a  level  flight  manoeu¬ 
vre  (run  82  of  flight  24,  LEV .  FLTuF0RWARD ,  u0 . 9VH) .  Only  the  first  0.481  seconds  of  this 
15.5  second  manoeuvre  are  shown,  that  is,  approximately  3%  of  the  time-history  data  are 
shown.  The  stress  for  the  remainder  of  this  manoeuvre  has  a  similar  structure  to  the  stress 
encountered  in  the  first  3%  of  the  manoeuvre.  The  complete  15.5  second  time-history  data 
(not  just  the  first  3%)  was  used  to  calculate  the  statistical  parameters  in  the  colour  coded 
legend  of  this  example. 

The  legend  in  Figure  5.2  summarises  the  most  relevant  statistical  parameters  calculated 
for  the  data  from  the  Flight  Loads  Survey.  The  definitions  of  the  different  symbols  found 
in  this  legend  are  given  in  Table  4.2  (on  page  25).  The  upper  and  lower  halves  of  this 
legend  are  not  exactly  symmetric.  The  uppermost  and  lowermost  points  of  the  red  region 
are  over-conservative  bounds  on  the  true  maximum  and  minimum,  respectively.  Although 
possibly  unconservative,  better  approximations  to  the  maximum  and  minimum  stresses 
are  given  by  the  extreme  points  of  the  green  section:  namely,  a niax  +  I'avg  and  am\n  ~  ^avg? 
respectively. 

The  stress  for  the  level  flight  group  is  shown  in  Figure  5.3  (for  further  details,  and 
the  plots  for  other  groups,  see  Appendix  D  on  page  83).  This  figure  shows  the  stress 
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for  all  manoeuvres  that  were  categorised  as  level  flight.  The  stress  is  sorted  by  average 
steady  stress  for  each  run,  which  explains  why  the  average  steady  stress  (depicted  by  the 
white  line)  is  monotonically  increasing.  The  vertical  axis  measures  stress  (in  ksi)  and 
the  horizontal  axis  enumerates  the  sorted  runs  (for  the  level  flight  group  there  were  338 
runs).  These  stress  measurements  were  for  bending  bridge  Number  1  on  the  right-hand 
stabilator,  which  is  9.2"  from  the  stabilator  inboard  edge  and  has  mnemonic  STBNBM1R 
(see  Figure  2.1  on  page  9). 

As  can  be  seen  from  Figure  5.3,  the  average  of  the  steady  component  <javg  varies  from 
approximately  -25  ksi  to  +5  ksi,  while  90%  of  the  directly  measured  stresses  (that  is,  s5% 
to  <s95%)  ranged  from  approximately  -30  ksi  to  +8  ksi. 

In  order  to  estimate  the  stress  on  the  cracked  stabilator  lug  we  make  two  assumptions. 
The  first  assumption  essentially  defines  a  “typical”  stress. 

Assumption  11  The  stresses  experienced  by  the  helicopter  (at  least  in  the  stabilator) 
during  level  flight  are  typical  of  normal  flying. 

This  assumption  of  a  single  load  amplitude  yields  a  coarse  approximation,  which  could  be 
refined  analogously  to  the  suggested  refinements  of  Assumption  10. 

The  second  assumption  goes  on  to  say  that  the  manoeuvres  within  this  typical  group 
are  homogeneous  in  stress: 


Time  (sec) 

Figure  5.2:  Legend  for  the  statistical  parameter  plots  of  bending  stress. 
Time-history  data  for  first  O.4.8I  seconds  of  the  15.5  seconds  it  took  to 
complete  this  level  flight  manoeuvre.  The  symbols  “o  ”  and  uv  ”  respec¬ 
tively  denote  the  steady  and  vibratory  component  of  the  stress,  while  the 
symbol  (ls  *  denotes  the  direct  measurement  of  stress. 
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Enumeration  of  runs 

Figure  5.3:  Plot  of  stabilator  spar  bending  stress  for  the  level  flight  ma¬ 
noeuvre  group.  These  stresses  were  sorted  by  the  average  steady  stress 
for  each  run.  These  stress  measurement  were  taken  at  bending  bridge 
Number  1 ,  which  is  located  9.2 "  from  the  stabilator  root. 


Assumption  12  The  helicopter  spends  the  same  amount  of  time  in  each  manoeuvre 
within  the  level  flight  group  and  each  of  these  manoeuvres  experiences  the  same  stress . 

Again,  this  coarse  assumption  was  required  to  produce  timely  results.  The  suggested 
refinements  given  for  Assumption  10  are  also  valid  for  Assumption  12. 

A  fatigue  strength  analysis  of  the  whole  stabilator  assembly9  by  Sikorsky  [1]  makes 
simplifications  similar  to  Assumptions  11  and  12.  Sikorsky  states:  (i)  that  transient  ma¬ 
noeuvres  did  not  produce  stresses  that  were  significantly  different  to  high  speed  level  flight; 
and  (ii)  that  the  percentage  of  time  spent  in  transient  manoeuvres  was  small  as  compared 
to  level  flight.  Using  these  two  observations  Sikorsky  [1,  p.  H-10]  concludes  that  fatigue 
testing  could  be  conducted  at  the  lug  stresses  observed  during  high  speed  level  flight. 

To  be  conservative  in  the  fatigue  analysis  that  follows,  we  want  to  obtain  a  typically 
high  fatigue  stress.  One  way  to  obtain  this  typically-high  steady-stress  is  to  consider  the 
95th  percentile  steady  stress  shown  in  Figure  5.3.  Negative  bending  stresses  imply  tension 
in  the  lug,  and  hence  we  want  the  5th  percentile  of  (jmin  to  use  in  fatigue  calculations. 

9The  stabilator  assembly  that  was  fatigue  tested  by  Sikorsky  consisted  of  the  left  and  right  stabilators 
and  the  centre-box. 
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From  338  runs,  the  5th  percentile  occurs  approximately  at  run  17  (see  Figure  5.3),  which 
has  a  minimum  steady  stress  of  <rmjn  ~  —23  ksi.  Using  a  similar  procedure  to  determine 
the  95th  percentile  run  and  its  corresponding  95th  percentile  vibratory  stress10  yields 
~  9  ksi,  which  is  the  typically-high  vibratory- stress.  In  summary,  the  typically  high 
stress  on  bending  bridge  Number  1  has  a  steady  component  of  —23  ksi  and  a  vibratory 
component  of  9  ksi.  Substituting  the  resulting  maximum  stress  Smax  =  —32  ksi  and  stress 
ratio  R  =  0.44  into  Equation  (5.3)  gives  the  equivalent  stress  as 

Seq  =  —23.4  ksi  (typically  high  loading).  (5.5) 

Again,  remember  that  a  negative  bending  stress  on  the  stabilator  implies  a  tensile  stress 
on  the  cracked  lug.  This  explains  why  the  maximum  stress  5max  is  more  negative  than 
the  steady  stress.  In  other  words,  the  operator  in  the  calculation  of  maximum  stress  is 
subtraction  not  addition,  that  is, 

Smax  =  (^miii)95%  —  (^9  5%)95%*  (5*6) 


The  95th  percentile  stresses  were  chosen  based  on  the  engineering  judgement  that 
these  percentiles  would  best  predict  fatigue.  In  contrast,  choosing  average  stresses  or  the 
maximum  stresses  would  result  in  fatigue  predictions  that  were  respectively  either  uncon¬ 
servative  or  too  conservative.  This  engineering  judgement  (to  choose  the  95th  percentile) 
requires  the  following  assumption: 


Assumption  13  The  95th  percentile  stress  (for  both  the  steady  and  vibratory  com¬ 
ponents)  yields  conservative  approximations  of  fatigue. 


This  assumption  was  made  to  obtain  a  less  conservative  fatigue  estimate  than  the  worst- 
case  scenario.  The  problem  with  Assumption  13  is  that  we  would  expect  any  single¬ 
amplitude  estimate  of  fatigue  loading  to  be  close  to  the  run-out  stress  of  the  component; 
that  is,  we  expect  the  single-amplitude  load  to  be  at  the  flat-end  of  the  fatigue  curve.  In 
this  run-out  region,  small  variations  in  loading  can  lead  to  large  variations  in  cycles  to 
failure. 

The  typical  high  stress  was  chosen  based  on  percentiles,  and  hence  there  is  an  implicit 
assumption  that  the  Flight  Loads  Survey  [3]  was  representative  of  in-service  usage.  More 
specifically,  that  the  number  of  runs  for  each  manoeuvre  carried  out  during  the  Flight 
Loads  Survey  is  representative  of  the  in-service  flight-time  fraction  of  those  manoeuvres. 
This  implicit  assumption  is  probably  reasonable,  but  not  for  the  reasons  that  may  first 
spring  to  mind.  The  Flight  Loads  Survey  was  carried  out  to  obtain  load  estimates  for 
critical  components — and  not  to  determine  the  Black  Hawk’s  usage  spectrum.  Thus  we 
might  expect  that  this  implicit  assumption  may  not  be  accurate.  However,  the  stress 
range  obtained  from  the  Survey  is  probably  representative  of  the  stress  range  components 
would  experience  in-service.  So  although  the  fraction  of  time  spent  in  different  manoeuvres 
during  the  Flight  Loads  Survey  is  not  representative  of  in-service  usage,  the  stress  range 
probably  is  representative  of  in-service  stress! 

For  expediency,  we  choose  a  single  flight  group  to  be  representative  of  typical  loads: 
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Assumption  14  The  “typically  high ”  stress  experienced  by  the  stabilator  during  the 
level  flight  group  of  manoeuvres  is  typical  of  in-service  flight  for  fatigue  purposes. 

Level  flight  consumes  most  of  the  in-service  flight  time,  and  also  consumes  a  fair  proportion 
of  fatigue  life  in  dynamically  loaded  components  from  the  rotor  system  (see  Krake  [16]  and 
Polanco  [22]).  We  would  then  expect,  from  a  fatigue  perspective,  that  using  level  flight  to 
approximate  typical  flying  is  a  reasonable  assumption. 


5.4  Best-Case  Bending  Stress  from  Flight  Loads  Survey 

In  this  section,  the  best-case  bending  stress  on  the  stabilator  is  determined  using 
measurements  from  the  Flight  Loads  Survey.  Remember  that  if  the  cracked  lug  is  still 
prone  to  fatigue  under  this  “gentle”  loading,  then  we  would  be  concerned  that  other 
Seahawks  may  also  be  susceptible  to  the  same  centre-box  lug  cracking. 

Using  engineering  judgement,  the  median  loading  should  represent  a  gentle  loading 
scenario.  In  other  words,  the  true  fatigue  life  of  the  component  should  be  shorter  than 
the  fatigue  life  predicted  using  the  median  loading.  Hence  the  median  loading  is  used  to 
determine  this  “best-case”  scenario. 

Since  by  definition  the  median  is  the  50th  percentile,  the  analysis  in  this  section  is 
analogous  to  the  analysis  in  the  previous  section.  The  only  difference  between  these  two 
sections  is  that  we  are  now  considering  median  steady  and  vibratory  loads  (as  opposed  to 
95th  percentile  steady  and  vibratory  loads  used  in  the  previous  section).  The  assumptions 
made  in  Section  5.3  are  also  required  in  this  section’s  analysis.  For  brevity,  however,  these 
assumptions  will  not  be  explicitly  repeated  in  this  section. 

From  the  bending  stress  plot  shown  in  Figure  5.3,  the  median  steady  load  occurs 
approximately  at  run  169  (which  is  half  of  the  338  runs  shown  in  this  plot).  Analogously 
to  the  previous  section,  see  Equation  (5.6),  the  maximum  stress  is  calculated  as 

^max  “  (<Javg)50%  “  (zyavg)50<v  , 

where  (x)r)07  denotes  the  50th  percentile  (that  is,  the  median)  of  the  argument  x.  The 
median  steady  and  vibratory  components  are  —5.5  ksi  and  3.5  ksi,  respectively,  which  give 
a  maximum  stress  of  Snmx  —  —  9.0  ksi  and  a  stress  ratio  of  R  =  0.22.  Substituting  these 
values  for  maximum  stress  and  stress  ratio  into  Equation  (5.3)  gives  the  equivalent  stress 

Seq  =  —7.9  ksi  (best-case  loading).  (5.7) 


10Unlike  the  high  steady  stress,  the  high  vibratory  stress  cannot  be  easily  determined  from  Figure  5.3. 
However,  the  procedure  for  the  determination  of  the  high  vibratory  stress  is  analogous  to  the  calculation 
of  the  high  steady  stress. 
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6  Fatigue  of  Cracked  Lug 

In  any  fatigue  analysis,  the  two  most  important  pieces  of  information  are  the  loading 
(i)  frequency  and  (ii)  amplitude.  In  this  section  we  estimate  both  these  quantities  from 
measurements  taken  during  the  Flight  Loads  Survey.  The  loading  frequency  is  estimated 
using  a  frequency  decomposition  of  a  level  flight  manoeuvre.  While  the  loading  amplitude 
is  estimated  from  the  geometry  of  the  cracked  lug  and  the  bending  stress  on  the  stabilator. 


6.1  Stabilator  Frequency  Response  during  Level  Flight 

In  this  section  the  loading  frequency  is  estimated  in  order  to  determine  the  cracked 
lug’s  susceptibility  to  fatigue  under  typical  loading  conditions. 

To  determine  the  stabilator’s  frequency  response  at  a  location,  both  the  input  loading 
and  the  resulting  frequency  response  are  required.  From  vibration  analysis  (see,  for  ex¬ 
ample,  Thomson  [26,  p.  380])  the  frequency  response  of  a  system  is  given  as  the  ratio  of 
output  to  input 

W(  \  _  _  FT  of  output 

10  X(uj)  FT  of  input 

where  a;  is  a  particular  frequency  and  H ,  X,  and  y  are  respectively  the  Fourier  transforms 
(FTs)  of  the  response,  input,  and  output  functions. 

For  the  Flight  Loads  Survey,  we  do  not  have  the  input  loading;  and  hence  cannot 
calculate  the  response  function.  Thus  we  require  the  following  assumption: 

Assumption  15  The  typical  loading  frequency  is  calculated  from  the  stabilator 's  fre¬ 
quency  response  during  a  level  flight  manoeuvre  from  the  Flight  Loads  Survey.  The 
typical  loading  frequency  is  defined  to  be  the  highest  frequency  out  of  all  the  large 
amplitude  loading  frequencies  (excluding  the  zero  frequency).  In  term,s  of  fatigue,  this 
typical  loading  is  assumed,  to  be  the  dominant  (that  is,  most  important)  loading. 

This  assumption  was  made  for  expediency.  A  more  accurate  fatigue  analysis  would  require 
cycle  counting  of  the  stress  history,  which  is  different  to  the  simple  frequency  decomposition 
undertaken  in  this  section. 

Under  Assumption  15  the  typical  loading  frequency  is  chosen  based  on  two  require¬ 
ments:  (i)  a  large  loading  amplitude  and  (ii)  the  highest  frequency  from  this  large  ampli¬ 
tude  group.  The  highest  frequency  aspect  of  this  assumption  should  make  this  assumption 
conservative—in  fact,  probably  too  conservative!  For  expediency  we  leave  a  more  refined 
analysis  for  future  work,  and  merely  acknowledge  the  over-conservatism  in  this  report. 

We  need  one  further  assumption,  this  time  about  how  the  typical  loading  frequency 
varies  among  different  manoeuvres: 


Assumption  16  The  typical  loading  frequency  is  the  same  for  all  manoeuvres,  and  is 
determined  from  the  stabilator  bending  stress  (STBNBM1R)  of  a  level  flight  manoeuvre 

(LEV . FLTuFORWARD , u0 . 9VH ). 
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This  level  flight  manoeuvre  was  deemed  to  be  typical  because  the  level  flight  group  of 
manoeuvres  consumes  the  largest  portion  of  flying  time.  (See  Krake  [16,  p.  25]  for  the 
Seahawk  usage  spectrum,  which  show  level  flight  consuming  approximately  70%  of  the 
flight-time.) 

The  frequency  responses  from  a  level  flight  manoeuvre  (see  Figure  6.1)  can  be  compared 
with  that  of  a  rolling  pullout  and  a  rough  approach  (see  Figures  E1-E2  on  pages  110- 
111).  This  comparison  of  the  highest  significant  frequencies  adds  weight  to  Assumption  16, 
namely,  that  the  different  manoeuvres  are  frequency  homogeneous  (at  least  in  terms  of 
significant  fatigue  frequency). 

The  frequency  response  of  the  stabilator,  at  bending  bridge  Number  1  (STBNBM1R),  is 
shown  in  Figure  6.1.  This  response  plot  is  for  the  level  flight  manoeuvre  from  run  82  of 
flight  24  (LEV .  FLTuF0RWARD ,  u0 . 9VH)  of  the  Flight  Loads  Survey.  During  the  Flight  Loads 
Survey  measurements  from  bending  bridge  Number  1  were  recorded  at  833  Hz,  which  gives 
the  Nyquist  frequency  of  416  Hz  shown  in  Figure  6.1. 

The  frequency  response  was  calculated  using  Mathematical  [27] 

Fourier  [w,FourierParaineters->{-l ,  1}] 


command,  which  returns  the  discrete  Fourier  transform  of  the  time-history  data  u .  The 
Fourier  transform  of  the  vector  u  is  given  by 


1  n 

Vk  -  -  Uj  exp  [2nd(k  -  1  )(j  -  l)/n]  . 


j= 1 


for  k  =  1,2  ,...,n, 


where  i  —  \f—\  and  uj  and  are  the  j th  and  kth  component  of  the  n-dimensional 
complex  vectors  u  and  u,  respectively. 

Because  the  time-history  data  are  real  (that  is,  u  €  Rn),  we  require  only  half  of  the 
complex  vector  v .  (The  second  half  of  this  vector  is  the  complex  conjugate  of  the  first 
half;  see,  for  example,  Press  et  al.  [23,  p.  497].)  Except  for  the  first  element  (which 
represents  the  mean  of  the  time-history),  the  first  half  of  the  discrete  Fourier  transform 
v  was  multiplied  by  2  in  order  to  compensate  for  the  discarded  second  half.  In  summary, 
the  frequency  response  figure  shown  in  this  section  plots  stress  amplitude  (|iq|  and  2\vk\ 
for  k  =  2, 3, . . . ,  n/2)  versus  frequency. 

The  thick  black  line  in  Figure  6.1  represents  a  polynomial  fit  of  the  frequency  data.  For 
numerical  accuracy,  this  polynomial  of  degree  twenty  was  fitted  as  a  Chebyshev  polynomial 
of  the  first  kind11  and  was  scaled  to  the  range  of  the  horizontal  axis. 

The  spikes  in  the  top  plot  of  Figure  6.1  occur  near  multiples  of  the  tail  rotor  frequency 
uju  which  is  19.8  Hz.  Multiples  of  the  tail  rotor  frequency  are  shown  by  the  red  vertical 
lines.  The  boxed  labels  attached  to  these  lines  denote  which  multiple  of  the  tail  rotor 
frequency  the  attached  red  line  represents. 

A  close-up  of  the  first  80  Hz  (shown  in  the  bottom  plot  of  Figure  6.1)  also  shows  the 
effect  of  the  main  rotor  frequency  ujm.  which  is  4.3  Hz.  The  thick  black  line  in  this  figure 
represents  a  first-kind  Chebyshev  polynomial  fit  to  only  the  first  80  Hz  of  the  frequency 
data. 

11  For  a  definition  of  Chebyshev  polynomials  see,  for  example,  Abramowitz  and  Stegun  [2,  p.  778] 
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Figure  6.1:  Frequency  response  of  stabilator  during  level  flight .  The 
bottom  plot  is  a  zoomed  section  of  the  dashed  box  in  the  top  plot.  The 
red  vertical  lines  labelled  with  a  box  denote  multiples  of  either  the  main 
or  tail  rotor  frequencies,  which  are  u)m  =  1^.3  Hz  and  u)t  —  19.8  Hz, 
respectively. 
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Surprisingly,  the  first  multiple  of  the  tail  rotor  frequency  is  almost  a  vibration  node; 
that  is,  the  vibrating  amplitude  at  19.8  Hz  is  lower  than  the  amplitude  at  the  surrounding 
frequencies.  However,  the  stabilator  does  strongly  resonate  at  the  second  and  fourth 
multiple  of  the  tail  rotor  frequency  (that  is,  at  2 uJt  =  39.6  Hz  and  4 ujt  —  79.2  Hz). 

The  resonating  frequencies  that  occur  at  approximately  62.3  Hz  and  367  Hz  do  not 
appear  to  coincide  with  multiples  of  either  the  main  or  tail  rotor  frequencies.  We  suspect 
that  these  resonating  frequencies  may  be  natural  frequencies  of  the  stabilator.  Alter¬ 
natively,  these  frequencies  may  be  due  to  aerodynamic  excitations  from  the  main  rotor 
blade — perhaps  buffeting? 

To  check  that  the  frequency  response  shown  in  Figures  6.1  was  typical,  two  additional 
manoeuvres  were  investigated.  The  results  for  these  manoeuvres  (a  rolling  pullout  and  a 
rough  approach)  are  shown  in  Appendix  E  on  pages  110  and  111. 

The  highest  frequency  within  the  high  loading  group  of  frequencies  in  Figure  6.1  occurs 
at  four  times  the  tail  rotor  frequency 

<fr hi  =  4 u>t  —  79.2  Hz 

sa  290  x  103  cycles/hour,  (6.1) 

which  is  the  frequency  of  the  worst-case  scenario.  As  has  already  been  mentioned,  this 
frequency  value  is  probably  too  conservative,  and  so  we  also  determine  a  lower  frequency 
value  to  reduce  this  conservatism. 

Figure  5.2  (on  page  33)  shows  a  time  history  plot  of  the  stabilator’s  bending  stress  for 
a  level  flight  manoeuvre  (see  Section  5.3  for  further  details  of  this  plot).  In  this  plot,  the 
peak-to-peak  timing  of  “significant”  stresses  is  approximately  0.1  seconds.  (A  significant 
stress  cycle  is  loosely  defined  to  be  a  cycle  that  exceeds  the  5th  and  95th  percentile  of 
directly  measured  stress.)  Inverting  the  period  of  this  cycle  gives  what  we  term  the  medium 
frequency  of  a  typically  high  loading 

Allied  =  10  Hz 

«  36xl03  cycles/hour.  (6.2) 

To  provide  some  perspective  on  the  two  loading  frequency  values  given  by  Equations  (6.1) 
and  (6.2),  the  fatigue  striation  marks  on  the  cracked  lug  were  analysed  fractographically. 

A  limited  amount  of  fractographic  data  was  generated  by  sampling  five  small  regions 
on  the  cracked  lug  (see  Appendix  G  on  page  120  for  further  details).  Extrapolating  the 
number  of  visible  striations  found  on  these  five  sampled  regions  provided  a  lower  bound 
on  the  number  of  loading  cycles  over  the  entire  crack  length,  namely,  20  000  cycles.  Thus 
a  lower  bound  for  the  loading  frequency  that  was  fatigue  damaging12  is  0.01  Hz,  which  is 
three  orders  of  magnitude  below  10  Hz  (the  lowest  of  the  two  loading  frequency  estimates 
calculated  above). 

One  explanation  for  the  discrepancy  between  these  loading  frequencies  is  that  the  dam¬ 
age  was  not  evenly  spread  over  the  service  life  of  402  hours,  but  instead  was  concentrated 
into  a  short  crack  life.  However,  the  crack’s  contamination  by  lubricant  suggests  a  crack 
life  of  weeks,  rather  than  a  crack  life  of  hours  (which  would  be  implied  by  a  10  Hz  loading 
frequency) . 

lzThis  fractographic  loading  frequency  was  calculated  as  20000/(402x60x60)  «  0.01  Hz. 
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If  we  consider  that  several  of  the  assumptions  used  in  this  analysis  had  an  undefined 
level  of  conservatism  (particularly  Assumption  11),  then  the  fractography  data  suggests 
that  these  were  extremely  conservative.  It  is  not  within  the  scope  of  this  report  to  reduce 
the  level  of  conservatism,  but  it  may  be  possible  to  target  these  reductions  in  future  efforts. 


6.2  Stress  and  Fatigue  on  Centre-Box  Lug 

In  this  section  we  determine  the  stress  on  the  cracked  lug  using  the  loading  estimates 
developed  in  earlier  sections.  First,  the  geometry  of  the  centre-box  lug  (which  attaches  to 
the  stabilator)  is  determined.  Then,  a  simple  load  analysis  determines  the  stress  on  this 
lug.  Finally,  the  fatigue  on  the  centre-box  lug  is  estimated. 

Figure  6.2  illustrates  the  two  relevant  (to  the  cracked  lug)  failure  modes  of  a  lug:  (i)  by 
tension  and  (ii)  by  shear  tear  out,  see  Bruhn  [4,  p.  D1.5]. 


Figure  6.2:  Relevant  failure  modes  of  a  lug:  (a)  tension  and  (b)  shear 
tear  out.  The  combined  horizontal  and  vertical  loading  experienced  by 
the  lug  is  shown  in  part  (c). 


As  is  shown  in  Figure  6.2(c),  the  stabilator  lug  experiences  both  horizontal  and  vertical 
loading.  From  the  theory  of  combined  stresses  (see,  for  example,  Roark  [28]),  the  principal 
stress  will  occur  either  horizontally  or  vertically  and  is  given  by 


<7oo 


max(ff ,  V) 


(6.3) 


where  H  and  V  are  the  horizontal  and  vertical  loads  experienced  by  the  lug,  and  A £  is  the 
lug’s  cross-sectional  area. 

The  areas  for  the  two  relevant  failure  modes  (tension  and  shear  tear  out)  are  calcu¬ 
lated  in  the  same  way.  Figure  6.3  shows  three  zoomed  sections  from  Sikorsky  drawing 
70209-27001  of  the  centre-box,  see  Figure  A16  on  page  67  (Figure  A17  may  also  be  of 
interest).  From  Figure  6.3,  the  relevant  cross-sectional  area  of  the  top  receiving  clevis  lug 
(shown  in  Box  C)  is 

At  =  4 x  (0.437x0.440)  =  0.769  in2.  (6.4) 
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Concentrating  on  the  forward  spar,  the  horizontal  and  vertical  loads  are  given  respec¬ 
tively  by  Equations  (3.8)  and  (3.24).  The  spar’s  length  is  L  =  74  (see  page  16),  and  for 
bending  bridge  Number  1,  xg  =  9.2  (see  Figure  2.1)  and  hf  =  5.9/2  (see  Table  2.1).  Thus 
the  equations  for  the  horizontal  and  vertical  loads  simplify  to 

Hf  =  -0.862 af  (6.5) 

and 

Vf  =  — 0.137crf.  (6.6) 

Remember  that  <jf  is  the  bending  stress  on  the  forward  spar  (in  this  case  at  bending  bridge 
Number  1).  The  cracked  lug  is  in  tension  when  the  bending  stress  on  the  spar  is  negative 
(see  Assumption  2),  which  explains  the  negative  signs  in  Equations  (6.5)  and  (6.6).  We 
have  used  Equation  (3.26)  in  preference  to  Equation  (3.24)  to  estimate  the  horizontal  load 
because  the  former  produces  the  higher  loads.13 

To  provide  some  perspective  of  our  loading  calculations  we  compare  the  vertical  loading 
given  by  Equation  (6.6)  with  the  design  loads  that  were  measured  by  Sikorsky  [9].  We 
will  compare  both  the  maximum  and  minimum  vertical  loads. 

From  Figure  5.1  we  see  that  the  maximum  bending  stress  (at  bridge  Number  1)  for  the 
Flight  Loads  Survey  occurs  in  a  manoeuvre  from  the  miscellaneous  group.  This  manoeuvre 
was  a  power  dive  from  run  22  of  flight  7  (with  mnemonic  P0WERUDIVEU1 . 2VH)  and  had  a 
maximum  bending  stress  of  —35.0  ksi.  Substituting  this  stress  into  Equation  (6.6)  gives 
the  vertical  loading  on  the  forward  spar  as  4800  lb.  In  the  same  way14  we  determine 
the  vertical  loading  on  the  aft  spar  as  2470  lb.  Adding  these  two  vertical  loads  gives  the 
maximum  vertical  loading  on  the  stabilator,  which  is  approximately  7270  lb.  According 
to  Sikorsky  [9],  however,  the  maximum  vertical  load  on  the  right  stabilator  was  1848  lb, 
which  is  approximately  one-quarter  of  the  7270  lb  load  we  would  predict.  We  will  discuss 
this  discrepancy  below. 

The  minimum  bending  stress  (at  bridge  Number  1)  for  the  Flight  Loads  Survey  occurs 
in  a  rolling  pullout  manoeuvre.  (These  minimum  stresses  are  not  shown  in  Figure  5.1.) 
This  manoeuvre  was  from  run  42  of  flight  29  (with  mnemonic  R0LLuP0uRT,u*8VH)  and 
had  a  minimum  bending  stress  of  24.0  ksi.  Using  the  same  procedure  as  in  the  previous 
paragraph,  the  minimum  vertical  load  on  the  right  stabilator  is  approximately  —4980  lb. 
The  Sikorsky  [9]  value  for  the  minimum  vertical  load  on  the  right  stabilator  was  —1552  lb, 
which  is  approximately  one-third  of  the  —4980  lb  load  we  would  predict. 

As  a  consistency  check  on  Sikorsky’s  loading,  two  different  Sikorsky  reports  are  com¬ 
pared  to  determine  the  consistency  of  the  maximum  loading  values  used  in  the  previous 
two  paragraphs.  Appendix  F  gives  the  details  of  this  consistency  check,  which  shows  that 
the  two  Sikorsky  reports  do  indeed  use  consistent  stabilator  loads. 

To  summarise,  the  different  vertical  loads  on  right  stabilator  (obtained  by  Sikorsky 
and  derived  in  this  report  using  the  Flight  Loads  Survey)  are: 

13Using  Equations  (3.24)  and  (6.6)  the  alternative  value  for  horizontal  loading  is  Hf  —  —0.594crf,  which 
is  a  lower  load  than  the  value  given  by  Equation  (6.5). 

14Using  Equation  (3.9)  and  the  value  ha  =  4.9/2  (from  Table  2.1)  the  vertical  loading  on  the  aft  spar  is 
given  by  Va  =  —  0.0707crf. 
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•  +1848  lb  (Sikorsky)  or  +7270  lb  (this  report)  for  the  maximum  load  and 

•  -1552  lb  (Sikorsky)  or  -4980  lb  (this  report)  for  the  minimum i  load. 

There  are  then  three  explanations  for  the  different  vertical  loads  shown  above:  (i)  the 
loads  estimates  we  derive  are  too  conservative,  (ii)  the  Flight  Loads  Survey  measurements 
are  erroneous,  or  (iii)  the  stabilator  loads  experienced  by  ADF  operated  Seahawks  are 
significantly  different  to  Sikorsky’s  test  helicopter.  It  is  unclear  which  of  these  options  is 
correct,  and  so  we  merely  note  that  a  difference  exists  between  Sikorsky’s  and  this  report’s 
bending  stresses. 

The  difference  in  maximum  vertical  load  between  this  report  and  Sikorsky  is  significant 
because  of  the  non-linear  relationship  between  stress  and  fatigue  life.  It  is  common  for 
the  gross  weight  of  the  Seahawk  to  be  approximately  20000  lb,  and  thus  the  magnitude 
of  the  stabilator  loads  are  a  significant  proportion  of  the  main  rotor  lift.  The  maximum 
vertical  loads  that  Sikorsky  use  would  not  fatigue  damage  the  investigated  lug.  However, 
the  intent  of  this  report  was  to  explore  a  methodology  for  airframe  fatigue,  and  not  to 
solve  a  specific  in-service  failure.  We  thus  leave  further  investigation  of  this  discrepancy 
(of  maximum  vertical  load)  for  future  work. 

The  “maximum”15  stress  experienced  by  the  cracked  lug  is  determined  from  Equations  (6.3) 
and  (6.4),  and  (because  the  horizontal  load  is  greater  than  the  vertical  load)  from  Equation  (6.5), 

<?c  =  H{/A( 

=  — L12<jf.  (6.7) 

In  Sections  5. 2-5. 4  we  obtained  values  for  the  bending  stress  on  the  stabilator ’s  forward 
spar,  see  Equations  (5.4),  (5.5)  and  (5.7).  The  stress  at  the  cracked  lug  can  be  calculated 
using  these  values  of  stress  in  Equation  (6.7).  Substituting  the  resulting  values  for  lug 
stress  into  the  cycles  to  failure  formula  given  by  Equation  (5.2)  yields: 


N{  =  24xl03 

(for  the  worst-case  loading), 

(6.8) 

N(  —  220  x  103 

(for  the  typically  high  loading), 

(6.9) 

and 

Arf  =  oo 

(for  the  best-case  loading). 

(6.10) 

Both  Equations  (6.8)  and  (6.9)  show  fewer  cycles  to  failure  than  the  run-out  cycles 

of  the 

MIL-HDBK-5  testing  [7,  p.  3-406],  and  hence  there  are  fatigue  implications.  In  contrast, 
Equation  (6.10)  shows  that  the  lug’s  fatigue  life  is  unlimited. 

Substituting  the  loading  frequencies  given  by  Equations  (6.1)  and  (6.2)  into  the  cycles 
to  failure  given  by  Equations  (6.8)  and  (6.9),  respectively,  yield  lives  of  approximately: 

•  5  minutes  for  the  worst-case  scenario, 

•  6  hours  for  the  high-loading  scenario,  and 

•  unlimited  hours  for  the  best- case  scenario. 

16 Maximum  stress  in  the  sense  of  the  95th  percentile. 
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It  must  be  remembered  that  (unlike  the  worst-case  scenario)  the  high  loading  scenario 
is  not  necessarily  a  conservative  bound  on  fatigue.  The  life  estimate  given  by  the  high 
loading  scenario  only  gives  us  a  feel  for  the  level  of  conservatism  inherent  in  the  worst-case 
estimate. 

The  best-case  scenario  was  included  in  case  it  resulted  in  a  finite  life  for  the  cracked 
lug,  which  would  have  suggested  the  lug  was  susceptible  to  fatigue.  As  explained  in  the 
introduction,  little  can  be  said  if  the  worst-case  analysis  results  in  a  finite  component  life 
or  if  the  best-case  analysis  results  in  an  unlimited  life.  Both  the  worst-  and  best-case 
scenarios  have  resulted  in  the  null  conclusion. 
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7  Summary 

The  aim  of  this  report  was  to  develop  a  methodology  for  assessing  the  fatigue  life 
of  helicopter  airframe  structure.  As  a  demonstrator,  we  estimated  the  fatigue  life  of  a 
centre-box  lug  found  on  the  Seahawk’s  stabilator. 

A  review  of  this  structure’s  design  shows  that  there  was  no  fatigue  damage  consider¬ 
ation  of  the  centre-box  lug  that  cracked.  The  design  appears  to  calculate  failure  margins 
based  on  static  loads  of  the  in-service  worst-case.  Only  the  stabilator  assembly  (composed 
of  the  wing  panels  and  centre-box)  was  fatigue  tested.  The  fatigue  sensitive  points  within 
the  stabilator  assembly  were:  (i)  the  forward  spar  from  the  right  wing  panel  and  (ii)  an 
aft  rivet  from  the  centre-box. 

The  information  required  for  this  fatigue  analysis  had  to  be  determined  indirectly,  and 
so  sixteen  assumptions  were  necessary.  The  two  coarsest  assumptions  were: 

•  that  the  worst-case  loading  applies  all  the  time  and 

•  that  cycle  counting  (for  fatigue  purposes)  can  be  approximated  using  the  frequency 
decomposition  of  the  loading. 

Although  probably  less  accurate,  the  frequency  decomposition  assumption  is  not  as  im¬ 
portant  as  the  worst-case  loading  assumption  for  the  following  reason:  The  cracked  lug 
was  not  a  lifed  component.  We  then  expected  the  loading  estimate  for  the  worst-case  to 
be  below  the  endurance  limit  for  the  cracked  lug — which  was  not  the  case!  For  loads  below 
the  cracked  lug’s  endurance  limit,  loading  cycles  would  be  irrelevant,  and  thus  the  worst- 
case  loading  assumption  is  more  important  than  the  frequency  decomposition  assumption. 
In  comparison  with  these  two  coarse  assumptions,  the  remaining  fourteen  assumptions  are 
more  accurate. 

The  estimated  fatigue  life  for  the  worst-case  scenario  was  low  for  several  reasons,  which 
are  all  related  to  the  degree  of  conservatism  inherent  in  the  assumptions.  Most  of  these 
assumptions  can  be  improved  by: 

•  Performing  a  full  fatigue  analysis  based  on  loading  amplitude  and  cycle  counting  of 
the  time-history  measurements. 

•  Or,  improving  the  measurement  information  on  which  the  fatigue  analysis  was  based. 

The  first  suggested  improvement  is  where  the  largest  gain  in  accuracy  (for  the  least  effort) 
can  be  achieved;  namely,  using  a  full  fatigue  analysis  to  refine  the  assumption  that  the 
worst-case  loading  is  acting  all  the  time  and  at  a  relatively  high  frequency. 

Low  and  unlimited  fatigue  lives  were  obtained  for  the  worst-  and  best-case  scenarios, 
respectively.  These  two  results  are  what  we  termed  the  null  conclusion,  and  so  only 
tentative  conclusions  can  be  drawn  from  these  results. 

The  typically  high  loading  (combined  with  a  loading  of  medium  frequency)  increased 
the  fatigue  life  as  compared  to  the  worst-case  scenario.  However,  even  under  this  relaxed 
scenario  the  cracked  lug  was  still  life  limited.  Only  the  best-case  scenario  resulted  in  a 
loading  that  was  below  the  endurance  limit  of  the  cracked  lug. 
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The  results  for  the  typically  high  loading  suggest  that  the  cracked  lug  may  be  suscepti¬ 
ble  to  fatigue  in  other  Seahawk  aircraft.  However,  this  fatigue  susceptible  conjecture  must 
be  tempered  by  the  results  from  the  best-case  scenario,  which  suggest  that  the  component 
may  indeed  have  an  unlimited  fatigue  life.  From  an  engineering  perspective,  a  component 
given  an  unlimited  life  would  imply  a  high  level  of  conservatism  in  design.  The  results 
we  have  obtained  in  this  report  suggest  that  this  conservatism  may  not  be  present  in  the 
cracked  lug.  Only  a  more  refined  analysis  (using  the  suggestions  made  throughout  this 
report)  could  draw  a  sharper  conclusion — either  for  or  against  the  unlimited  life  of  the 
cracked  lug.  It  must  be  stressed,  however,  that  within  the  fidelity  of  results  obtained  so 
far,  the  work  carried  out  in  this  report  does  not  conflict  with  Sikorsky’s  assessment  of  an 
unlimited  fatigue  life. 

Experience  of  fleet  service  life  provides  some  background  assessment  to  this  analysis. 
The  ADF  has  a  fleet  of  52  Hawk- variant  helicopters  that  share  the  same  folding  stabilator 
design,  with  the  fleet  leader  having  accumulated  over  3000  flight  hours  since  new.  This 
failure  mode  in  the  cracked  lug  has  not  occurred  in  any  other  of  these  52  helicopters. 
Furthermore,  using  the  authors’  limited  access  to  failure  data  no  other  instance  of  this 
failure  type  was  uncovered  in  the  remaining  worldwide  fleets.  This  worldwide  fleet  is 
composed  predominantly  of  the  US  Navy’s  fleet  of  approximately  300  helicopters,  where 
the  fleet  leader  has  reached  10000  flight  hours. 

Despite  the  evidence  of  stabilator  buffeting  recorded  during  the  Flight  Loads  Survey, 
there  is  little  available  operational  evidence  that  suggests  that  this  buffeting  commonly 
manifests  itself  in  extensive  fatigue  damage  to  this  part  of  the  stabilator.  In  contrast, 
however,  cracking  of  the  trailing  edge  skins  and  excessive  wear  in  the  stabilator  actuator 
attachment  bushes  are  both  common  problems  afflicting  the  Hawk  fleet.  These  common 
problems  were  not  evident  in  the  fatigue  test  conducted  by  Sikorsky  but  could  be  expected 
to  be  exacerbated  by  buffeting. 

We  recommend  the  refinement  of  the  fatigue  estimation  methodology  developed  in  this 
report.  In  particular,  this  refinement  should  include: 

•  a  detailed  fractographic  study  of  the  lug’s  fracture  surface  (to  identify  crack  growth 
behaviour),  and 

•  a  fatigue  life  calculation  based  on  cycle- counting  of  the  Flight  Loads  Survey  data. 

However,  since  both  these  activities  are  resource  intensive,  the  scheduling  for  such  work 
is  dependent  upon  a  review  of  current  priorities. 
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Appendix  A  Engineering  Drawings 

This  appendix  contains  engineering  drawings  of  some  stabilator  components.  Each 
caption  gives  the  drawing  number  (dn),  sheet  number  (sht:X  of  Y),  and  filename  (f)  of 
the  drawing.  For  example,  the  caption  in  Figure  A5  has  the  code  (dn:70202-07051 ,  sht:  1 
of  2,  rev:K,  f:0003658A.CAL),  which  stands  for  revision  K  of  sheet  1  (of  two  sheets)  from 
drawing  number  70202-07051 ,  with  associated  filename  0003658A. CAL. 

In  October  2003,  these  drawings  were  available  on  the  Air  Vehicle  Division’s  network 
computer,  in  the  directory  I :  \BlackuHawkuDrawings\SIKORSKY\ IMAGES. 

The  electronic  version  of  this  report  contains  the  complete  versions  of  these  drawings. 
Thus  any  section  of  these  drawings  may  be  zoomed  to  obtain  greater  detail. 
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(dn:70202-27001,  sht:l  of  2,  rev:F,  f:0003689A.CAL) 


Figure  A9:  Rib  of  stabilator  at  BL  28. 
(dn:70201-07052,  sht:l  of  1,  rev:A,  f:0003568A.CAL) 
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Figure  All:  Rib  of  stabilator  at  BL  9,  including 
(dn:70201-27002,  sht:2  of  2,  rev:F,  f:0003596A .  CH 
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Figure  A15:  Stabilator  beam  at  BL  5.40  (cent 
(dn:70201-07056,  sht:l  of  1,  rev:K,  f:0003574A 
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Appendix  B  Twelve  Manoeuvre  Groups 

This  appendix  lists  the  twelve  groups  into  which  the  manoeuvres  from  the  Flight  Loads 
Survey  were  partitioned  (see  Section  4.1  for  further  details). 


_ Table  Bl:  List  of  manoeuvres  in  the  autos  group. _ 

Autos 

auto^desc,^vbroc,M20yonr  auto^desc ,  vma,  120°/onr  auto^turnJLt ,  . 8vma,30deg 

auto^desc,^vbroc,  J)270nr  auto^desc , vma,92°/0nr  auto^turnJLt ,  .8vma,45deg 

auto^desc  ,^vh,^lOO0/onr  auto^entry ,  8vh  autojturn^lt ,  vma,30deg 

auto^desc ,  .8vh,  10070nr  auto^entry  8vma  auto^turnJLt ,  vma,45deg 

auto^desc ,  . 8vh,  116yonr  auto^entry , ^vh  auto^turn^rt ,  . 8vh, 30deg 

auto^desc,  .8vh,1200/onr  auto_entry  ,^vma  auto^turn^rt , . 8vh,45deg 

auto^desc,  .8vh,928/0nr  auto^recvry 8vh  auto^turn^rt , . 8vma,  15deg 

auto^desc,  .8vma,  1008/0nr  auto^recvry  8vma  auto^turn^rt ,  .8vma,30deg 

auto^desc ,  ,8vma,  120°/0nr  autojrecvry ,^vh  auto^turn^rt , . 8vma,45deg 

auto^desc ,  . 8vma,928/cnr  auto^recvry , ~vma  auto^turn^rt ,  vma,30deg 

auto^desc ,  vma ,  100%nr  auto^turnJLt ,  .  8vh ,  30deg  auto^turn^rt ,  vma , 45deg 

auto^desc ,  vma,  1088/onr  auto^turnJLt ,  . 8vh,45deg 

auto^desc  ,vma,  112°/enr  auto^turn  JLt ,  . 8vma,15deg 


Table  B2:  List  of  manoeuvres  in  the  climb  group. 

Climb 


climb , ..vbroc , ..irp 
climb , ..vbroc ,  ancp 
climb , ..vbroc+15 , ..irp 
climb ,  ..vbroc+15 ,  jncp 
climb ,  ..vbroc-15 , ..irp 
climb ,  _vbroc-15 ,  jncp 
cltrn..lt..l5 ,  vbroc , irp 
cltrn~lt~15, vbroc,  mcp 
cltrn..lt..l5 , vbroc+15 , irp 
cltrn..lt..l5,  vbroc+15 , mcp 
cltrnJ_t..30 , vbroc , irp 


cltrn..lt..30 ,  vbroc ,mcp 
cltrnJLt  J30,  vbroc+15  ,mcp 
cltrn^lt^45 , vbroc , irp 
cltrn..lt  ..45,  vbroc ,mcp 
cltrn^lt^45 ,  vbroc+15 ,  irp 
cltrnJLt..45,  vbroc+15,  mcp 
cltrnJLt^45 ,  vbroc-15 ,  irp 
cltrn..lt..45 ,  vbroc-15  ,mcp 
cltrn^rt ^15 , vbroc , irp 
cltrn..rt^l5 , vbroc ,mcp 
clt  rn..rt..l  5 ,  vbroc+15 , irp 


cltrn..rt..l5 ,  vbroc+15  ,mcp 
cltrn^rt  J30 ,  vbroc ,  irp 
c It  rn^r t  J30 , vbroc , mcp 
cltrn^rt  J30 ,  vbroc+15  ,mcp 
cltrn..rt..45 , vbroc , irp 
cltrn^rt^45 , vbroc , mcp 
cltrn^rt^45 ,  vbroc+15 ,  irp 
cltrn.xt.A5 ,  vbroc+15 , mcp 
cltrn.xt..45 ,  vbroc-15,  irp 
cltrn_rt^45 ,  vbroc-15,  mcp 
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Table  B3:  List  of  manoeuvres  in  the  hover  group. 


Hover 

hover  .Tge ,  ,T00T/*nr 
hoverx>ge ,  T.00  J/*nr 
hover «oge , T.oad..ext  ended 
hover^oge , ..lower  ..load 
hover^oge , ..pick^upT.oad 

hover^oge , .naiseTxad 
hover^oge , _setTLoad 
hover^oge ,  ..set -download 
hover^oge , ..sling , ..pickup 
hover^oge , ..withJLoad 

hover^turn^lef  t ,  ..15x1/ s 
hover  ^turn^left ,  J30xl/s 
hover  ..turner  ight ,  ..15..d/s 
hover ^turn^right ,  J30..d/s 
hover , ..cableT: ull-up 

Table  B4:  List  of  manoeuvres  in  the  level  flight  group. 


lev.f  ltjf  orward,  J)  .3vh 
lev .  f  It  J:  orward ,  J) .  4vh 
lev .  fit  ^forward ,  T) .  5vh 
lev.f It ^forward,  J)  .6vh 
lev . flt^f orward ,  JD .  7vh 
lev .  fit  ^forward ,  J) . 8vh 
lev . fit ^forward ,  J3 . 9vh 
lev .  f  It  J:  orward ,  ^vh 


_ Level  flight 

lev.  f  lt^fwd,^.  6vh,  ^9570nr 
lev .  f  lt~f  wd ,  6vh ,  ^97°/0nr 
lev.f  lt^fwd,^.6vh,^9970nr 
lev .  f  lt-f  wd ,  6vh ,  10l7eiir 
lev .  f  It  J^wd,  9vh,-9570nr 
lev .  f  lt^fwd,^.  9vh ,  w97#/0nr 
lev.f  lt^fwd,^.9vh,~9970nr 
lev .  f  It  wd,  9vh,  100e/0nr 


lev.flt~fwd,^.9vh,  10l70nr 
level^f  light ,  .TOO  Jtias 
level  J: light ,  J30  Jsias 
level^f  light  ,  J70~kias 
level^f light ,  J30  Jcias 
level^f light ,  J90  Jsias 


Table  B5:  List  of  manoeuvres  in  the  miscellaneous  group. 


ambient 

appr .  runon ,  40kt ,  aer  o  J>rk 
approach  Jc_runon ,  _40kt 
approach ,  ^normal 
approach,  ^operational 
approach ,  trough 
close^aptudoor 
dashJt-quickstop 
e/s/rJLt«trn, .8vh,30deg 
esss^4001^1hs  Jc^rhs^tank 
esss_4001TLhsT;ank 
esss^4001^1hs/8001^rhs 
esssT3001JLhs.J^rhsT;ank 
e  s  s  s  ^emp  t  y  T:  ue  1 
ilhsJ2001Jfc.TrhsT)l 
ilhs^4001^&^irhs^01 
ilhs^4001._&T.rhsJ2001 
ilhs^4001Jc^irhs^4001 
ilhs.T001Tc.Trhs  J5001 
ilhs.T001  JcT.rhs  J3001 
ilhs  .TOO  lTc.Trhs  ..8501 
ilhs..6001Tc.Trhs.J3501 
ilhsT3001Tc.Trhs..8501 
ilhs^8501^&^irhs^8501 
jump^take-of f 
landing , ..ship 
lhsJ31bs  Jt..rhsJ)lbs 


_ Miscellaneous _ 

lhs~2001bsTcj:hsJ)lbs 

lhs.T001bs  Jc^rhs^Olbs 

lhs.T001bs  Jc^rhs^2001bs 

lhs.T001bsTc..rhs.T001bs 

lhs.T001bsTc.Jrhs.i>001bs 

lhs.T001bs  JC^rhsTSOOlbs 

lhs^4001bs  Jc^rhs^8401bs 

lhs..6001bs  Jt^rhsTTOlbs 

lhs..800.TbsTc..rhsJ3401bs 

lhs^8601bs^&^rhs^8401bs 

lift^of  f  ^toJiover^ige 

lsf ..accelJc  ..quick-stop 

maintainer.xlimbingT.hs 

maintainer.x  limbing  .xhs 

mntnr^stnd.xlsd.JT^cowl 

mntnr^stnd.xlsd..#2.xowl 

mntnr^stnd^lhs^apu^door 

mntnr^stnd^open^#l^.cowl 

mntnr^stnd^open^#2^cowl 

mntnr^stnd^rhs^apu^door 

mntnr  ^walk.x  lsd~#  1  ~c  owl 

mntnr..walk.xlsd^#2.xowl 

mnt nr ^wal k ^on_apu^d o o r s 

mntnr^walk^open^#l^cowl 

mntnr^walk^open^#2^cowl 

no . l^engine_start 

no . 2^engine^start _ 


opening^apu^door 
power^dive^l . lvh 
power^dive^l . 2vh 
ppd , ^90kts , TOOOf pm 
ppd ,  J^Okts ,  T500f  pm 
ppd ,  J50kts , T>00f  pm 
ppd , ^vh , .T  50  Of  pm 
r  c  vr  y  Jf  r  om^par t  ^pwr  ^de  s  c 
rotor^engagmnt ,  j£  ly 
rotor^engagmnt ,  .jioTxrake 
rotor^engagmnt ,brk 
rotor^engagmnt , brk ,  Jiard 
rotor^engagmnt , brk , Tdle 
rotor  ^shutdown ,  aio^brake 
rot  or  ..shut  down , brk 
rot  or  ..shutdown ,  brk ,  407oiir 
rot  or  ..shut  down , brk , hard 
rsf..accelTc  ..quick-stop 
rwd^accelJc^quick-stop 
take-off  ..jump 
take-off  .normal 
take-off  trolling 
takeoff , .normal 
terrain^cyclic^pull 
terrain^cyclic^push 
vert icalTLanding 
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Table  B6:  List  of  manoeuvres  in  the  symmetric  pullouts  group. 

Symmetric  pullouts 


incl^sym^po , . 8vh , It , 2 . lg 
incl^sym^po ,  .  8vh , It , 2 . 4g 
incl^sym^po ,  .  8vh ,  It ,  2 . 5g 
sym^pullout , .  8vh ,  JL.  5g 
sym^pullout ,  .  8vh ,  ^1 . 6g 
sym^pullout , . 8vh ,  J. .  8g 
sym^pushover , . 8vh , ^ . 25g 
sym^pushover ,  .  8vh ,  w .  4g 
sym^pushover , . 8vh , - . 6g 
symm^pullout , ^vh , ^1 . 3g 
symm^pullout , ^vh , ^1 . 45g 
symm^pullout , ^vh , . 4g 
symm^pullout , ^vh , ^1 . 55g 
symm^poillout ,  ~vh ,  . 5g 

symm^pullout , ^vh , ^1 . 65g 
symm^pullout , ^vh , ^1 . 6g 
symm^pullout ,  ^vh ,  A . 73g 
symm^pullout , ~vh , . 75g 
symm^pullout , ^vh , ^1 . 7g 
symm^pullout , ~vh , A . 85g 
symnkpullout ,  ^vh ,  „1 . 8g 
symm^pullout ,  ^vh ,  -1 . 9g 
symm^pullout ,  „vh ,  J1 .  Og 
symm^pullout ,  ^vh ,  J2 .  lg 
symm^pullout , .  8vh ,  A . 25g 


symm^pullout , . 8vh , „i . 35g 
symm^pullout , . 8vh , . 3g 
symm^pullout ,  .  8vh ,  . 4g 

symm^pullout ,  .  8vh ,  ^1 . 5 lg 
symm^pullout , . 8vh , ^1 . 55g 
symm^pullout , . 8vh , „i . 5g 
symm^pullout , .  8vh ,  _1 . 65g 
symnupullout ,  .  8vh ,  Jl .  68g 
symm^pullout ,  .  8vh ,  ^1 . 6g 
symm^pullout , . 8vh , ^1 . 75g 
symm^pullout , . 8vh ,  JL  7g 
symm^pullout , . 8vh , A . 85g 
symm^pullout ,  .  8vh ,  A . 8g 
symm^pullout ,  .  8vh ,  J .  91g 
symm^pullout ,  .  8vh ,  A  .  95g 
symm^pullout ,  .  8vh ,  J .  9g 
symm^pullout , .  8vh ,  J2 .  Og 
symm^pullout , .  8vh ,  J2 .  lg 
symm^pullout , . 8vh ,J2.2g 
symm^pushover ,  ^vh ,  - .  lg 
symm^pushover ,  jvh , « .  21g 
symm^pushover , ^vh , „ . 25g 
symm^pushover ,  ^vh ,  - .  2g 
symnupus  hover ,  ..vh ,  w .  35g 
symm^pus  hover ,  ^.vh ,  „ .  3g 


symm^pushover , 

«vh>M 

•45g 

symm^pushover , 

^vh,^ 

•4g 

synmupus  hover , 

^vh,^ 

55g 

symm^pushover , 

-vh,^ 

57g 

symnupushover , 

“Vh,^ 

•5g 

symm^pushover , 

“Vh,^ 

•% 

symm^pushover , 

~vh,~ 

•  75g 

symm.pushover , 

-vh,- 

•7g 

symm.pushover , 

-vh,u 

•8g 

symm^pushover , 

.8vh, 

-16g 

symm^pushover , 

.8vh, 

“•lg 

symnupushover , 

.8vh, 

-21g 

symnupoishover , 

.8vh, 

“•  25g 

symm^pushover , 

.  8vh, 

-2g 

symm^pushover , 

.8vh, 

“•  35g 

symm^pushover , 

.  8vh, 

“•38g 

symm_pushover , 

.  8vh, 

“•3g 

symuLpu  shove  r , 

.  8vh, 

“•45g 

symmjpushover , 

.  8vh, 

-•47g 

symm^pushover , 

.  8vh, 

“•4g 

symm^pushover , 

.  8vh, 

-•  57g 

symm^pushover , 

.  8vh, 

„.5g 

symm^pushover , 

.  8vh, 

-6g 

symm^pushover , 

.8vh, 

“•7g 
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Table  B7:  List  of  manoeuvres  in  the  rolling  pullouts  group. 


Rolling  pullouts 


roll^po^(verify) 
roll^po_2t  ,^.8vh 
roll^po  JLt , « . 8vh , ^1 . 25g 
roll^poJLt  ,^.8vh,.T .35g 
roll^poJLt  ,^.8vh,.T  .3g 
roll^poJLt  ,^.8vh,_l  .4g 
roll^po^lt , „ . 8vh , ^1 . 55g 
roll^po^lt ,^.8vh,~l . 5g 
roll^po  JLt , 8vh , ^1 . 65g 
roll^po  JLt , ^ . 8vh , -1 . 6g 
roll^po  JLt , ~ . 8vh , . 75g 
roll^po  JLt , ^ . 8vh , . 7g 
roll^po  JLt ,  ^ .  8vh ,  . 85g 

roll^po  JLt , ^ . 8vh ,  JL .  8g 
roll^po  JLt , ~ . 8vh , ^1 . 95g 
roll^po  JLt , ~ . 8vh , -1 . 9g 
roll^pOwlt ,^.8vh,-2.0g 
roll^po  JLt ,  ~ .  8vh ,  J2 . 12g 
roll^poJLt  ,...8vh,  J2. lg 
roll^po  JLt ,  _ .  8vh ,  J2 . 25g 
roll-po  JLt ,  _ .  8vh ,  J2 . 2g 
roll^po  JLt ,  8vh ,  J2 . 3g 
roll^po^lt ,^.8vh,^2.4g 
roll^po  JLt ,  ^vh. ,  ^1 . 25g 
roll^po  JLt , ^vh , -1 . 2g 
roll^po  JLt , jv h , . 35g 


roll^po  JLt , ^vh , . 3g 
roll^po  JLt , ^vh , -1 . 45g 
roll^po  JLt , _vh , . 4g 
roll^poJLt ,  ^vh ,  ..1 . 55g 
roll^po  JLt , _vh , -1 . 5g 
roll»po*J.t , ^vh , . 65g 
roll^po  JLt , ^vh , ^1 . 6g 
roll^po  JLt , ^vh , . 75g 
roll^po  JLt , ^vh , ^1 . 7g 
roll^po  JLt , ^vh , . 85g 
roll-po  JLt ,  ^vh ,  ^1 . 8g 
roll-po  JLt , ^vh , ^1 . 95g 
roll^po  JLt , _vh , ^1 . 9g 
roll-po  JLt ,  ^vh ,  J2 .  Og 
roll^po  JLt ,  ^vh ,  J2 . 2g 
roll-poj:t 8vh,  JL .  25g 
roll^po jrt  y~.  8vh,^l  .3g 
roll^po^rt ,^.8vh,^l .45g 
roll^pojrt  8vh,^l  .4g 
roll^po jrt 8vh, ^1 . 55g 
roll^po^rt 8vh, . 5g 
roll^po jrt  ,~.  8vh,  «1 . 65g 
roll^po^rt 8vh,^l .6g 
r oll^po^rt , _ . 8vh , ^1 . 75g 
roll._po.jrt  8vh,^l .  77g 
roll._po.jrt , « .  8vh ,  . 7g 


roll-pojrt ,  _ .  8vh ,  JL .  8g 
roll-pojrt  ,_.8vh,_l  .  95g 
roll^po^rt 8vh,  JL . 9g 
roll-po ~rt ,  ^ .  8vh ,  J2 .  Og 
roll^po^rt ,  ~ .  8vh ,  J2 .  lg 
rollupojrt ,  u .  8vh ,  J2 . 2g 
r  oll^po^rt ,  ^ .  8vh ,  J2 . 3g 
r  oll^po^rt ,  _ .  8vh ,  J2 . 4g 
r  olL-po^rt , ^vh , ^1 . 25g 
r oll^po^rt , ^vh , . 2g 
r  oll_po.jrt , ^vh , ^1 . 35g 
r oll^po^rt ,  jvh , ^1 . 3g 
r oll^po^rt , ^vh ,  J. . 45g 
roll^po^rt , ^vh , -1 . 4g 
roll^po^rt , ^vh , ^1 . 55g 
roll^po^rt , ^vh , ^1 . 56g 
roll^po^rt , wVh , -1 . 5g 
roll^po^rt , ^vh , . 65g 
roll^po^rt , ^vh , ^1 . 6g 
roll^po^rt , ^vh , ^1 . 75g 
roll^po^rt , jvh , ^1 . 77g 
roll^po^rt , ^vh , -1 . 7g 
roll^po^rt ,  ^.vh ,  ^1 . 85g 
roll^po^rt ,  ^vh ,  Jl .  8g 
roll^po^rt ,  ^vh ,  J2 .  Og 
rt . ^roll^pullup^vh^l . 5g 


Table  B8:  List  of  manoeuvres  in  the  reversal  group. 


col .  rev^auto ,  .  8vh,  100°/*n 
col .  r  ev^auto ,  vma ,  10070nr 
col . rev^lev . fit , 8vh 
col . rev^lev . f It ,  ^vh 
col .rev,  Jiover 
lat .  rev^auto ,  .  8vh ,  1000/0n 
lat .  rev^auto ,  vma ,  100#/onr 


_ Reversal _ 

lat .  revJLev .  fit ,  8vh 
lat . rev^lev . fit , ^vh 
lat .  rev,  Jiover 
lng .  rev^auto ,  .  8vh ,  10070n 
lng .  rev^auto ,  vma ,  1007«nr 
long . rev^lev . fit , 8vh 
long .  revJLev .  f  It ,  ^vh 


long. rev, Jhover 
rud .  rev^anto ,  .  8vh ,  10070n 
rud .  rev^aut  o ,  vma ,  10079nr 
rudd. revJLev. f It 8vh 
rudd . rev^lev . f It , ^vh 
rudd . rev ,  Jiover 
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Table  B9:  List  of  manoeuvres  in  the  side  and  rearward  flight  group. 
Side  and  rearward  flight 


lt^side^f It^accel , ^lOkts 
lt^side^J:  lt^accel ,  ^20kts 
lt^side^f lt^accel , JSOkts 
lt^side^f lt^accel , ~35kts 
lt^side^f  Itu.accel ,  ^40kts 
lt^side^f lt^accel , ^45kts 
rear^f lt^accel , ^lOkts 
rear^jf  lt^accel ,  ^20kts 
rear^f It^accel ,  J25kts 
rear^f lt^accel , ~30kts 
rear^f lt^accel , ^35kts 
rear  Jilt^accel ,  ^40kts 
rear^f It^accel , ^45kts 
rearwards It , ^lOkts 


rearward^f  It ,  w.20kts 
rearwards It ,  J30kt  s 
rearwardT It , J35kts 
rearwards It , ^40kts 
rearward^  It ,  ^45kts 
recovry^left^side^flt 
r e  covry^r earwardT 1 ight 
recovryjright^side^f  It 
rt^side^lt^accel ,  ^lOkts 
rt^side^f lt^accel , ^20kts 
rt^side^f lt^accel , ^25kts 
rt^side^flt^accel ,  J30kts 
rt^side^f lt^accel ,  J35kts 
rt^side^f lt^accel , ^40kts 


rt^side^f lt^accel , ^45kts 
side^f  ItJLef  t , ^lOkts 
side^f It  JLef t ,  J20kts 
side^f  lt^lef  t ,  J30kts 
side^f lt^lef t ,  J35kts 
side^jf  lt^lef  t ,  ^40kts 
side^f lt^lef t , ^45kts 
side^f It  aright ,^10kts 
side^f It  aright , ^20kts 
side^f It  aright , J25kts 
side^f It  aright ,  J30kts 
sideJ  It  aright  ,^35kts 
side^J:  It  aright  ,^40kts 
side^f It ight , ^45kts 


sideslipJLt , 
sideslipJLt , 
sideslipJLt , 
sideslipJLt , 
sideslip  JLt , 
sideslip  JLt , 
sideslip  JLt , 
sideslip  JLt , 
sideslipJLt , 
sideslipJLt , 
sideslipJLt , 


Table  BIO:  List  of  manoeuvres  in  the  sideslip  group. 


^.6vh,~i0deg 
^.6vh,~15deg 
~ . 6vh ,  J20deg 
6vh,  J25deg 
_>.6vh,^5deg 
^.8vh,~10deg 
^.8vh,^15deg 
8vh,^20deg 
8vh,  J22deg 
8vh,^25deg 
^.8vh,^5deg 


Sideslip 

sideslip  JLt  9vh,  J5deg 

sideslip  JLt ,^vh,  ^lOdeg 
sideslip  JLt , ^vh , ^5deg 
sideslipjrt ,  6vh ,  ^lOdeg 
sideslip^rt , „ . 6vh , ^15deg 
sideslip^rt , ^ . 6vh , J>0deg 
sideslip^rt , ^ . 6vh , ^25deg 
sideslip  jrt , 6vh , J>deg 
sideslip  _rt , w. 8vh , _5deg 
sideslip^rt 8vh, ^lOdeg 
sideslip^rt , 8vh , ^15deg 


sideslip^rt , 8vh,  J20deg 
sideslip^t  ,~.8vh,  J25deg 
sideslip^t , 8vh ,  J5deg 
sideslip^rt  ,^vh,^10deg 
sideslip^rt , ^vh , w5deg 
sideslip^trim,^. 6vh 
sideslip^trim,^. 8vh 
sideslip^trim , . 8vh 
sideslip  J;rim,_vh 


Table  Bll:  List  of  manoeuvres  in  the  taxi  group. 


Taxi 

air^taxi^accel , ^20kts 

air^taxi^f  wd ,  -ige ,  J20kts 

taxi^s tart 

air^taxi^accel , J30kts 

air^taxi^f  wd ,  ^ige ,  J30kts 

taxi cStop 

air^taxi^accel , ~40kts 

air^taxi^f  wd ,  ^ige ,  ^40kts 

taxi_straight 

air^taxi^accel , JSOkts 

air^taxi^fwd, ^ige , ^50kts 

taxi^turn^left 

air^taxi^f wd , ^ige , ^lOkts 

air^taxi^recovry 

taxi^turn^right 
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Table  B12:  List  of  manoeuvres  in  the  turns  group. 


desc^trn^lt , 90kts , 15deg 
desc^trnJLt , 90kts , 30deg 
desc^trn  JLt ,  90kts ,  45deg 
desc^trmjrt , 90kts , lBdeg 
desc^trn^rt , 90kts , 30deg 
desc^trn^rt  ,90kts  ,45deg 
entry^lt^tn , .  8vh ,  15deg 
entryJLt^tn,  .8vh,30deg 
entryJLt^tn,  .8vh,45deg 
entryJLt^tn,  ,8vh, BOdeg 
entryJLt^tn,  .8vh, BBdeg 
entry  JLt^tn , .  8vh ,  60deg 
entryJLt^tn,vh,  lBdeg 
entry  JLt^tn ,  vh ,  30deg 
entry^lt^tn ,  vh ,  45deg 
entry  JLt^tn ,  vh ,  50deg 
entry  JLt^tn ,  vh ,  55deg 
entryJLt^tn ,  vh ,  60deg 
entry-rt^tn,  .8vh, lBdeg 
entry^rt^tn,  . 8vh,30deg 
entry^rt^tn,  . 8vh,4Bdeg 
entryjrt^tn,  .8vh,50deg 
entryjrt^tn,  .8vh, BBdeg 
entry^rt^tn,  .8vh,60deg 
entry  ^rt^tn ,  vh ,  15deg 
entry^rt^tn ,  vh ,  30deg 
entry^t^tn,  vh,4Bdeg 


_ Turns _ 

entry  ^rt^tn ,  vh ,  50deg 
entry  ^rt^tn ,  vh ,  55deg 
entry  _rt^tn ,  vh ,  60deg 
rapid^decel^turn ,  JLef t 
rapid^decelwturn ,  ..right 
recovry^lt^tn , ,8vh, lBdeg 
recovry^lt^tn, .8vh,30deg 
recovry^lt^tn, ,8vh,45deg 
recovryJLt^tn,  .8vh, BOdeg 
recovry^lt^tn,  .8vh, BBdeg 
recovry^lt^tn,  .8vh,60deg 
recovry^lt^tn, vh, 15deg 
r  ecovry  JLt>_tn , vh , 30deg 
r ecovry^lt^tn , vh , 45deg 
recovry^lt^tn , vh , BOdeg 
recovry  JLt..tn , vh , BSdeg 
recovry  Jlt^tn ,  vh ,  60deg 
recovryjrt^tn, .8vh,lBdeg 
recovry^t^tn,  .8vh,30deg 
recovry^rt^tn, .8vh,4Bdeg 
recovry^rt^tn,  .8vh,B0deg 
recovry^rt^tn, .8vh,BBdeg 
recovry^rt^tn, .8vh,60deg 
recovry^rt^tn,  vh,  lBdeg 
recovry^rt^tn ,  vh ,  30deg 
recovry^rt^tn , vh , 4Sdeg 
recovry^rt^tn , vh , BOdeg 


recovry^rt^tn ,  vh ,  BBdeg 
recovry^rt^tn , vh , 60deg 
stdyturnJLt, .8vh,lBdeg 
stdyturn  J.t , . 8vh , 30deg 
stdyturmlt,  .8vh,4Bdeg 
stdyturnJLt , . 8vh , BOdeg 
stdyturn  J.t, .8vh, BBdeg 
stdyturnJLt , . 8vh , 60deg 
stdyturnJLt , vh , lBdeg 
stdyturnJLt , vh , 30deg 
stdyturn^lt , vh , 4Bdeg 
stdyturn^lt , vh , BOdeg 
stdyturnJLt ,  vh ,  BBdeg 
stdy  tumult ,  vh ,  60deg 
stdyturn  j:t , . 8vh , lBdeg 
stdyturn^rt , . 8vh,30deg 
stdyturnjrt , . 8vh,4Bdeg 
stdytnrn^rt , . 8vh, BOdeg 
s t dy turner t , .8vh, BBdeg 
stdyturn^rt , . 8vh,60deg 
stdyturn^rt , vh , lBdeg 
stdyturn^rt , vh , 30deg 
stdyturn^rt , vh , 4Bdeg 
stdyturn^rt , vh , BOdeg 
stdyturn^rt , vh , BBdeg 
stdyturn^rt ,  vh ,  60deg 
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Appendix  C  Mathematica  Code  for  Bending 

Stress  Plots 

This  appendix  contains  the  code  used  to  generate  the  bending  stress  plots  shown  in 
Section  5.3  and  in  Appendix  D.  The  Mathematica  code  is  listed  in  Figures  C1-C7.  In 
both  the  group  plots  (shown  in  Figure  C6)  and  the  manoeuvre  plots  (shown  in  Figure  C7) 
only  the  first  plot  from  the  set  of  twelve  plots  is  shown.  Each  of  these  cells  was  collapsed 
since  these  plots  are  shown  elsewhere  in  this  report  (see  Appendix  D) . 
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stahStrainJ.nh  * 

Plots  of  Stabilator  Bending  Strain 

(using  data  from  Black  Hawk  Flight  Loads  Survey) 

Work  carried  out  for  Stabilator  lug  cracking  on  Seahawk  (Jut  —  Aug  2003) 

in  [l]  :=  <<  LinearAlgebra'MatrixManipulation'; 

<<  Graphics 'Colors"; 

<<  Graphics 'Graphics3D'; 

<<  myColors.m; 

In  [5 ]  :*  Off [General: :spelll] 

In  [6]  SetDirectoryfToFileName [Extract! 

"FileName"  /.  NotebookInformation[EvaluationNotebook[ ] ] ,  {1},  FrontEnd'FileName] ] ] 

Out  16}  e  C:\Data,  Black  Hawk  strain  survey\Stabilator\StatData 

■  Importing  Manoeuvre  Groups  (grouping  according  to  influence  on 
dynamic  components) 

In(7):*  manvDyn  =  Map[DeleteCases[#,  ""]&,  Transpose  [Import  [ "manvDyn. csv"  ,  "CSV"]]]; 

Map [Length,  manvDyn] 

OutfSJ-  (38,  33,  16,  23,  81,  75,  79,  21,  43,  32,  16,  81} 

■  Partitioning  Manoeuvres  into  Groups 

In  [9 ):*  fnPos[y_,  man_J  :=  Module  [  (pos) , 
pos  =  Position[v,  man]  ; 

Flatten [ {First [pos]  ,  Last [pos]}] 

] 

In[10]  partitionManv[mtx_,  manv_]  :  =  Module [  {manvCol}  , 

(*  Partitions  a  matrix  into  manoeuvres  *) 
manvCol  =  Transpose  [mtx]  [[1]] ; 

Table  [Take  [mtx,  fnPos  [manvCol,  manv[[i]]  ]  ]  ,  {i.  Length [manv]  }  ] 

3 

In  [11]:“  readGauge[file_J  :=  Module [ {srt} , 

(*  Reading  in  comma- separated  data  file  *) 
srt[v_]  s  =  Sort [v,  {#115]]  <#2[I5]])  &]; 
dat  =  Map  [Take  [#,  13]  &,  Import  [file,  "CSV"]]; 
colNam  =  datpj ; 

Print [TableForm[Transpose [ {Range [Length [colNam] ] ,  colNam} }  ]  ] ; 
dat  =  Rest [dat] ; 

Print  ["Number  of  manoeuvres^"  ,  Length  [manvNam  =  Union  [Transpose  [dat]  [[1]]  ]]]  ; 

Print [ " dimens ions= " ,  Dimensions [dat] ] ; 

dat = DeleteCases [dat,  {a _ ,  *'",  b _ },  co]  ; 

statDat  =  Map [srt,  partitionManv[dat ,  manvNam] ] ; 

Lengthf statDat] 

] 


17:27:53  Thursday.  August  7.  2003 


Figure  Cl:  Mathematica  code  for  stress  plots,  Page  1. 
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stabStrain2.nb 


In [12]  .-=  readGauge["stbnbmlr.csv"] 

1  MANEUVEP._DESC 

2  AVG_VIB 

3  MAX_VIB 

4  PCL_VIB 

5  AVG_STDY 

6  MIN_STDY 

7  MAX_STDY 

8  ABMIN_STDY 

9  ABMAX_STDY 

10  PCL_DIR 

11  PCC_DIR 

12  FL I GHTTE  S T_NUM 

13  RUN_NUM 

Number  of  manoeuvres- 52 6 
dimensions- {3796,  13} 

24.445  Second 
Out  [12]“  526 


■  Generic  Plot  for  Set  of  Data 

In  [13]  .■  =  ptPoly [pTop_,  pBot_,  x_,  dx_]  :=  Module  [  {npts,  ptsTop,  ptsBot}, 
npts  =  Length[pTop] ; 

ptsTop  =  Table  [  {x  +  i  dx,  pTopfi]]},  {i,  npts}]; 
ptsBot  =  Table [ {x  +  i  dx,  pBotliJ},  {i,  npts}]; 

Join[ptsTop,  Reverse [ptsBot]  ] 

]; 

ptLine[pts_,  x_,  dx]  :=  Module  [  {npts}  , 
npts  =  Length[pts] ; 

Table [ {x  +  i  dx,  ptsffil},  {i,  npts}] 

] 


maxVibCol  1 

Red  1 

percVibCol 

Yellow 

avgVlbCol 

Green 

percDirCol 

RoyalBlue 

maxStdCol 

Black 

avgStdCol  ; 

,  White 

17:27:54  Thursday,  August  7.  2003 


Figure  C2:  Mathematica  code  for  stress  plots,  Page  2. 
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stdbStrain2.nb 


In  [16]  -  picHanv  [mtx_ ,  x,  dx_]  :  = 

(*  Draws  the  statistical  data  for  a  manoeuvre  *) 

Module [ {man,  avgVib,  maxVib,  percVib,  avgStd,  minStd,  maxStd, 

absMinStd,  absMaxStd,  percDirMax,  percDirMin,  flgtNum,  runNum) , 

{man,  avgVib,  maxVib,  percVib,  avgStd,  minStd,  maxStd,  absMinStd,  absMaxStd, 
percDirMax,  percDirMin,  flgtNum,  runNum}  = Transpose [mtx] / 1000. 7 
If  [Length [mtx]  ==  1, 

{Thickness [0 .001]  ,  PointSize[0.002] , 

maxVibCol,  Line  [ptPoly  [maxStd  +  maxVib,  minStd  -  maxVib,  x,  dx]  ] , 
percVibCol,  Line  [ptPoly  [maxStd  +  percVib,  minStd  -  percVib,  x,  dx]  ] , 
avgVibCol,  Line (ptPoly [maxStd +  avgVib,  minStd  -  avgVib,  x,  dx]  ]  , 

Dashing [{ (*0.05,0.01*) >] ,  percDirCol, 

Line [ptPoly [percDirMax,  percDirMin,  x,  dx] ] , 

Dashing[{ (*0.003,0.0075*) }] ,  maxStdCol, 

Point  [ptLine [maxStd,  x,  dx]|[l]I],  Point  [ptLine  [minStd,  x,  dx]  IIJ]  , 
Dashing!  {}]  ,  avgStdCol,  Point  [ptLine  [avgStd,  x,  dx][[l]|] 

{Thickness[0 . 001] , 

maxVibCol,  Polygon  [ptPoly  [maxStd  +  maxVib,  minStd  -  maxVib,  x,  dx]  ]  , 
percVibCol,  Polygon  [ptPoly  [maxStd  +  percVib,  minStd  -  percVib,  x,  dx]  ] , 
avgVibCol,  Polygon  [ptPoly  [maxStd  +  avgVib,  minStd  -  avgVib,  x,  dx]  ] , 
Dashing [{ (*0.05, 0.01*))] ,  percDirCol, 

Polygon [ptPoly [percDirMax,  percDirMin,  x,  dx] ] , 

(★Line [ptLine [percDirMax, x,dx] ] , Line [ptLine [percDirMin, x, dx] ] , *) 
Dashing [{ (*0.003, 0.0075*) )] ,  maxStdCol, 

Line [ptLine [maxStd,  x,  dx] ] ,  Line [ptLine [minStd,  x,  dx]  ] , 

Dashing[ O ] ,  avgStdCol,  Line [ptLine [avgStd,  x,  dx] ] 

)] 


17:27:54  Thursday.  August  7,  2003 


Figure  C3:  Mathematica  code  for  stress  plots,  Page  3. 
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stabStrain2.nb  4 


■  Legend  for  Plots 

In[l7]:=  legnd =  Module [ {avgVib,  maxVib,  percVib,  avgStd,  minStd,  maxStd,  absMinStd, 
absMaxStd,  percDirMax,  percDirMin,  x=  400,  dx  =  5,  vec,  pi,  p2}, 
vec  =  {6663 . 813",  9568.719',  8526.507',  -16644.44',  -18101.65', 

-15120.44',  -26852.94',  -7042.328',  -11593.64',  -21639.94'}; 

{1800,  4000,  3000,  8000,  6000,  9400,  3400,  11400,  10000,  4700); 

{avgVib,  maxVib,  percVib,  avgStd,  minStd,  maxStd,  absMinStd, 
absMaxStd,  percDirMax,  percDirMin}  = Transpose! {vec,  vec}]; 
p2  =  Show[Graphics[{ 

Text["maxVibMax" ,  {x  +  2  dx,  (maxStd  +  maxVib)  [[1J  }  ,  {-1,  0}], 
Text["percVibMax",  {x+2dx,  (maxStd  +  percVib)  {[131 } •  {-1#  0 >  ] , 

Text  ["avgVibMax",  {x  +  2dx,  (maxStd  +  avgVib)  1 1J } ,  {-1,  0)], 

Text  [“percDirMax" ,  {x  +  2dx,  percDirMax  [I]] } ,  {-1,  0}], 

Text  ["maxStd",  {x  +  2dx,  maxStdllJ},  {-1,  0}], 

Text  [  "avgStd"  ,  {x  +  2dx,  avgStd[[l]|  }  ,  {-1,  0}], 

Text [ "minStd" ,  {x  +  2dx,  minStd[[ll},  {-1,  0}], 

Text  ["percDirMin",  {x  +  2  dx,  percDirMin  [[1J } ,  {-1,  0}], 

Text["avgVibMin" ,  {x  +  2dx,  (minStd  -  avgVib)  [[11 } ,  {-1,  0}], 

Text  [  "percVibMin" ,  {x  +  2  dx,  (minStd  -  percVib)  [[II }  ,  {-1,  0}], 

Text  [  "maxVibMin" ,  {x  +  2  dx,  (minStd  -  maxVib)  [[11 }  ,  {-1,  0}], 

Thickness [0.001] , 

maxVibCol,  Polygon [ptPoly  [maxStd  +  maxVib,  minStd  -  maxVib,  x,  dx]  ] , 
percVibCol,  Polygon  [ptPoly  [maxStd  +  percVib,  minStd  -  percVib,  x,  dx]  ] , 
avgVibCol,  Polygon  [ptPoly  [maxStd  +  avgVib,  minStd  -  avgVib,  x,  dx]  ] , 

Dashing[ { (*0 . 05, 0.01*) }] ,  percDirCol, 

Polygon [ptPoly [percDirMax,  percDirMin,  x,  dx] ] , 

Dashing[{ (*0.003,0.0075*) }] ,  maxStdCol, 

Line [ptLine [maxStd,  x,  dx] ] ,  Line [ptLine [minStd,  x,  dx] ] , 

Dashing [{}],  avgStdCol,  Line [ptLine [avgStd,  x,  dx] ] 

}] ,  Di  splay  Function-*  Identity] ; 

Print ["time  history=". 

Dimensions [datTimHist  =  ReadList [ "stbnbmlrTimeHistoryT24R082.dat"] ] ] ; 

(*  Choosing  the  first  400  pts  to  show  trend  only  *) 

pi  = 

ListPlot  [Take  [datTimHist,  400],  PlotJoined -*  True,  DisplayFunction -♦  Identity] ; 
Show[{pl,  p2}.  Frame-*  True,  FrameLabel  -*  {"time",  "BndStrnPSI" ,  ,  ""}, 

DisplayFunction-*  $DisplayFunction, 

PlotRange  -*  All,  Axes  -»  False,  FrameTicks  -►  None] 

] 

time  history- {12907} 


17:27:54  Thursday,  August  7,  2003 


Figure  C4:  Mathematica  code  for  stress  plots,  Page  4. 
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stabSfrain2.nb 


Out  [17 ]  =  -Graphics- 

■  Plots  of  Bending  Strain  for  Groups 
(sorted  by  average  bending  stress  of  run) 

In  [18]:=  manvGrp  = 

Sort  [Table  [  {Position [manvDyn,  statDat[[i,  1,  1J][[1,  1J  /  i}*  {i»  Length [ statDat] }]] ; 
manvGrp  =  Map  [Transpose  [#]  [[2]]  &,  Split  [manvGrp,  (#111]]  ==#2[[1J)  &]  ] ; 

in [20]:*  picGrp  [k_]  :=  Module  [  (grp,  mtx,  n) , 
n  =  Leng th[ manvGrp HkJ]  ; 

grp  =  Table  [j  =  manvGrp |[k,  ij  ;  statDat JjJ,  {i,  n}]; 
mtx  =  Sort  [Flatten [grp,  1],  (#1[[5]]  <  #2[[5]])  &]; 

Show  [Graphics  [picManv[mtx,  0,  1]  ] ,  Frame-*  True, 

PlotRange  ->  All,  FrameLabel -*  {manvDyn Jk,  1J  ,  "BendingStrainKSI" } ] 

] 

In  [21] :  =  pGrp  =  Table [picGrp [i] ,  {i.  Length [manvDyn] >] 


17:27:54  Thursday,  August  7,  2003 


Figure  C5:  Mathematica  code  for  stress  plots,  Page  5. 
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stabSlrain2.nb 


6 


■  Plots  of  Bending  Strain  for  Groups 
(sorted  by  average  bending  stress  of  run  within  a  manoeuvre) 

In  [22 ]:=  picMultManv[mtx_,  vec_,  lbls_]  :=  Module [  {nvec,  dx  =  0.1,  x=0,  pics,  xtks} , 
nvec  =  Length[vec]  ; 

pics  =  Table  [If  [i  >1,  x  +=  dx  Length [mtx[[j  J  ]  ,  0]  ;  j  =  vecftij  ; 

picManv[mtx[[j]l ,  x,  dx]  ,  {i,  nvec}]; 
x  =  0; 

xtks  =  Table  [j  =  vec|[i]] ;  {x  +=  dx  Length  [mtx[[;j]]  ]  ;  x  -  dx  Length  [mtx[[j]|  ]  /  2 , 
ToString [Position [manvNam,  mtx[[j  ,  1,  l]j][[l,  1]]]},  {i,  nvec}]; 
ShowfGraphics  [pics]  ,  Frame-*  True,  PlotRange  -»  All,  FrameTicks  ^ 

(xtks.  Automatic,  None,  Automatic},  AspectRatio -*■  .5,  FrameLabel -»  Ibis] 

] 

In  [23]  : -  pMnv  =  Table  [picMultManv [ statDat,  manvGrp[[i]] , 

{manvDyn[[i,  1J  ,  "BendingStrainKSI" }  ]  ,  {i.  Length [manvGrp]  }] 


1 7:27:54  Thursday.  August  7.  2003 


Figure  C6:  Matliematica  code  for  stress  plots,  Page  6. 
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sntbStniinlnh 


-  id  '• 

171  a *120  21  22  2  i .mi:  27  Z2r?n  32  Si  3-1  35  SC  37  3?  39  4 QM12  4  3  4 41 54 64 71439  5  0  51 

auto:: 


1  -  Graphics  -  ,  -  Gi apbics  - ,  -  Graphics  -  ,  -  Graphics  -  ,  -  Graphics  -  ,  -  Graphics  - 
. Graphics  - ,  -  Graphics  - ,  -  Graphics  -  ,  -  Graphics  -  ,  -  Graphics  - ,  -  Graphics  - 


■  Exporting  Plots  as  EPS  Files 

.1  jj  !  :  -  Export [ "legend. eps"  ,  legnd,  "EPS"] 

Our  [.:■;]  legend  .  eps 
Jn/.’f.;.--  Table  [ 

Export  [  "grp "<>  StringReplace[manvDyn[[i,  1J  ,  {"  "  ->  "/"  ->  ""  }  ]  <>  ".eps", 

pGrp[[i]] ,  "EPS"], 

{ i ,  Length [manvDyn] } ) 

; it; f  /.'M  ,  grpAUTOS.eps,  yrpCLIHB.eps,  grpHOVER.eps,  grpLEVEL.eps, 

qrpMlSC.eps,  grpPVHPULLOUT.eps,  grpROLLPULLOUT.eps,  grpREVERSAL. eps, 
grpS 1  PE REAR  . eps ,  grp? T DE SLIP . eps ,  grpTAXI . eps ,  grpTURNS . eps  j 

in Table  [ 

Export  ["mnv"  <>  S  tringRep  lace  [manvDyn  [i ,  lfl  ,  {  "  eps"  , 

pMnvJi]],  "EPS"], 

{i,  Length [manvDyn] } ] 

p!:f  ••  i  nmvAlJTGS .  eps  ,  nmvCLIME .  ops ,  mnvHOVER .  eps ,  mnvLEVEL .  eps , 

mnvKlSC.eps,  mnvSYKPULLOUT.eps,  mnvROLLPULLOUT.eps ,  mnvREVERSAL.eps, 
mnv?.  I  PER  EAR .  eps ,  mnvS  IPE3L1P .  eps  ,  mnvTAXl .  eps ,  mnvTURNS  .  eps  • 
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Fnniiv  C7:  Mat Roma  lira  code  for  stress  plots.  Page  7. 
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Appendix  D  Bending  Stress  from  Flight  Loads 

Survey 

This  Appendix  contains  plots  of  the  stabilator  bending  strain  at  a  location  that  is 
9.2"  from  the  root.  These  “statistical  summaries”  of  stress  were  obtained  from  the  Flight 
Loads  Survey  of  the  Black  Hawk  [8].  By  statistical  summaries  we  mean,  for  example,  the 
average  steady,  maximum  vibratory,  and  95%  percentile  stress.  For  a  legend  to  these  plots 
and  further  details  of  these  statistical  summaries  see  Section  5.3,  which  begins  on  page  32. 
In  particular,  for  an  illustration  of  the  legend  see  Figure  5.2. 

Most  of  the  526  manoeuvres  contain  multiple  runs,  but  a  few  manoeuvres  contain  only 
one  run.  There  were  some  missing  data  in  these  statistical  summaries,  in  particular,  of 
the  3796  runs  there  were  31  runs  that  were  missing  all  the  statistical  data  except  for  the 
95%  and  5%  percentile  of  the  direct  stress  measurement.  The  list  of  manoeuvres  with 
incomplete  data  is  shown  in  Table  Dl.  Note  that  the  stresses  in  this  table  are  given  in  psi 
(and  have  been  rounded  to  three  significant  figures),  whereas  the  stresses  in  all  the  plots 
within  this  appendix  are  given  in  ksi 

The  first  twelve  plots  in  this  appendix,  Figures  D1-D12,  show  the  bending  stresses 
sorted  by  the  average  steady  stress  for  each  run  of  the  Flight  Loads  Survey.  All  manoeuvres 
were  partitioned  into  one  of  twelve  groups.  The  horizontal  axis  in  these  plots  enumerates 
each  run  within  a  particular  group,  which  gives  us  an  idea  of  the  number  of  runs  within  a 
group. 

The  next  twelve  plots,  Figures  D13-D24,  show  the  bending  stresses  partitioned  into 
manoeuvres,  and  then  sorted  by  the  average  steady  stress  for  each  run  within  a  manoeuvre. 
In  these  plots,  the  tick  marks  on  the  .r-axis  denote  the  manoeuvre.  Each  of  these  tick  marks 
has  been  labelled  with  the  manoeuvre  name  in  a  very  small  font  (which  was  necessary  to 
avoid  overlap  of  the  manoeuvre  names).  These  manoeuvre  names  will  be  illegible  when 
the  plots  are  printed  on  a  standard  A4  page.  However,  these  names  are  legible  when  the 
electronic  version  of  this  report  is  viewed  in  a  zoomed  mode. 
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Table  Dl:  List  of  the  31  runs  that  were  deleted  in  the  plots  contained 
in  this  appendix.  These  runs  were  deleted  because  they  were  missing 
statistical  data  for  the  steady  and  vibratory  components  of  the  run.  Only 
the  5%  and  95%  percentile  direct  stress  measurements  (s5%  and  s,,s%9 
respectively)  were  available  for  these  manoeuvres. _ 


Flight 

Run 

Manoeuvre 

Sr.%  (PSI) 

S„r,%  (Psi) 

22 

4 

ambient 

178 

354 

24 

8 

ambient 

940 

1230 

25 

4 

ambient 

1050 

1240 

2G 

4 

ambient 

-2.44 

286 

28 

4 

ambient 

-554 

-295 

29 

4 

ambient 

568 

857 

30 

4 

ambient 

1010 

1270 

32 

4 

ambient 

277 

445 

34 

4 

ambient 

-290 

-97 

36 

4 

ambient 

-77.8 

211 

37 

4 

ambient 

-508 

-316 

39 

4 

ambient 

-304 

853 

40 

4 

ambient 

456 

840 

41 

4 

ambient 

1160 

1350 

42 

4 

ambient 

207 

378 

43 

6 

ambient 

325 

517 

44 

4 

ambient 

536 

708 

45 

4 

ambient 

24.3 

313 

46 

4 

ambient 

-453 

-164 

48 

4 

ambient 

961 

1150 

48 

16 

air  taxi  fwd,  ige, 

lOkts 

-3460 

845 

48 

18 

air  taxi  fwd,  ige, 

20kts 

-3320 

967 

48 

19 

air  taxi  accel,  30kts 

-3420 

767 

49 

4 

ambient 

766 

958 

50 

4 

ambient 

182 

471 

52 

4 

ambient 

474 

666 

52 

33 

roll  po  rt,  . 8vh, 

1  •  5g 

-2330 

7430 

52 

34 

roll  po  rt,  .8vh, 

1  •  5g 

-1410 

6870 

52 

35 

roll  po  rt,  .8vh, 

1  •  5g 

-3080 

7180 

52 

36 

roll  po  rt,  .8vh, 

1.5g 

-3770 

9080 

53 

4 

ambient 

59.0 

237 

84 


DSTO-TR  1590 


(is^)  ssaa^s  Suipuag 


Figure  Dl:  Plot  of  stabilator  spar  bending  stress  for  the  manoeuvre  group 
of  autos.  These  stresses  were  sorted  by  the  average  steady  stress  for  each 
run.  These  strain  measurement  were  taken  at  the  zeroth  gauge ,  which  is 
located  9.2"  from  the  stabilator  root. 
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Figure  D3:  Plot  of  stabilator  spar  bending  stress  for  the  manoeuvre  group 
of  hover.  These  stresses  were  sorted  by  the  average  steady  stress  for  each 
run .  These  strain  measurement  were  taken  at  the  zeroth  gauge,  which  is 
located  9.2"  from  the  stabilator  root 
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Figure  D21:  Plot  of  stabilator  spar  bending  stress  for  the  manoeuvres 
group  of  sideslip.  These  stresses  were  partitioned  into  manoeuvres , 
and  then  sorted  by  the  average  steady  stress  for  each  run  within  that 
manoeuvre.  These  strain  measurement  were  taken  at  the  zeroth  gauge , 
which  is  located  9.2 "  from,  the  stabilator  root. 


SYMMETRIC  PULLOUT 


DSTO  TR  1590 


Appendix  E  Stabilator  Frequency  Response: 
Rolling  Pullout  and  Rough  Approach 

In  Section  6.1  the  frequency  response  of  the  stabilator  to  level  flight  was  investigated. 
In  order  to  check  that  those  results  represent  a  "typical”  frequency  response,  we  investigate 
the  frequency  response  of  two  further  manoeuvres  in  this  appendix. 

The  two  manoeuvres  investigated  were: 

•  a  rolling  pullout  manoeuvre  from  flight  45  run  44  (R0LLuP0uLT,uVH,u2 . 2G)  and 

•  a  rough  approach  from  flight  69  run  42  (APPROACH  ,uR0UGH). 

These  manoeuvres  were  chosen  as  the  antithesis  of  level  flight  for  the  following  reasons. 
The  left  rolling  pullout  was  chosen  because  it  was  thought  that  the  left  turn,  combined 
with  the  pullout  motion,  would  force  the  unsteady  downwash  (from  the  main  rotor)  onto 
the  stabilator.  In  contrast,  the  rough  approach  would  provide  frequency  information  for 
low  speed  flight  combined  with  continuously  fluctuating  stabilator  movements. 

Figures  El  and  E2  show  the  frequency  response  of  the  stabilator  respectively  for  the 
rolling  pullout  and  rough  approach  manoeuvres. 

As  in  Section  6.1.  the  thick  black  lines  represent  a  twentieth  order  Chebyshev  poly¬ 
nomial  through  the  frequency  response  shown  in  the  plot.  The  red  lines  denote  selected 
multiples  of  the  main  and  tail  rotor  frequencies.  The  box  that  labels  these  red  lines  denotes 
which  multiple  of  the  rotor  frequency  each  line  represents. 

For  a  more  detailed  analysis  of  these  plots  see  the  comments  made  for  the  frequency 
response  of  the  level  flight  manoeuvre  in  Section  6.1. 

The  Mathematica  code  used  to  generate  these  frequency  response  plots  is  shown  in 
Figures  E3  E6 
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Figure  E2:  Frequency  response  of  stnbilator  loading  to  rough  approach. 
The  red  vertical  lines  labelled  with  a  box  demote  multiples  of  the  main 
or  tail  rotor  frequency,  which  are  ujm  =  f  .3  Hz  and  cvt  —  19.8  Hz , 
respectively. 
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stabFreq.nb 


Stabilator  Loading  Frequency 

(using  data  from  Black  Hawk  Flight  Loads  Survey) 

Work  carried  out  for  Stabilator  lug  cracking  on  Seahawk  (Jul  -  Aug  2003)  ] 


In  [l]  :=  «  myColors  .m; 

<<  Graphics 'Graphics'; 

in [3] :  =  Off [General: :spelll] ; 

Off [General: : spell] ; 

In  [5]  :=  SetDirectory[ 

ToFileName [Extract [ "FileName "  /.  No tebooklnformat ion [EvaluationNot ebook [ ] ] , 
{1}  ,  FrontEnd' FileName]  ]  ] 

Out(5)=  C:\MyFolder\Work  for  Others\Boykett  Stabilator\StatData 


■  Plotting  Functions  J 

In[€]:=  maxAmp [v_,  ul_,  w2_]  :=  (*  max  amplitude  between  two  frequencies  *) 

First [Reverse! 

Sort [DeleteCases[Abs  [v]  /.  {w_,  amp_)  -»  If  [ul  <  u  <  oi2,  {amp,  to},  Null]  ,  Null]  J  ]  ] 

In  [7]  :=  lgPlt[v_,  opts _ ]  :=  LogListPlot  [v,  PlotRange  -*  All,  ImageSize  -*  500 ,  1 

Frame  -+  True,  Axes  ^  False,  FrameLabel-*  {"Frequency",  "StressAmpli tude" }  ,  opts]  J 

InfSJ:=  logTks[xmin_,  xmax_,  lst_:{l}]  :=  Flatten [Table [ 

{Log [10 . ,  j  10 .i] ,  If [MemberQ[lst,  j ]  ,  ToString[ j ]  <>  "E"  <>  ToString[i] ,  ""]}, 

{i,  Floor [Log[10,  xmin] ] ,  Ceiling[Log[10,  xmax] ] } ,  {j,  9}],  1] 
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Figure  E3:  Mathematica  code  for  frequency  plots,  Page  1. 
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stabFreq.nb 


2 


In  [9]  ;  = 


pltDat[pts_,  n_,  polyOrd_,  tail_]  :  = 

Module[{v,  fn,  pi,  p2,  xO,  xl,  chebPoly,  xa,  xb,  ya,  yb,  tks,  dsh, 


tksX,  tksY,  tailRotorFreq =  19.8,  rotorTxt,  yO ,  colRotor  =  Red} , 
rotorTxt [i_]  :=  {colRotor, 

Text [ToString [i] ,  {tailRotorFreq*  i,  yO},  {0,  -1},  {0,  1}]}; 


chebPoly[k_,  x_]  ;  =  chebPoly[k,  x]  =  ChebyshevT[k, 


2  {x-xO) 

{xi -xor 


-1]; 


dsh = Dashing [{0.001,  0.01}]; 

yO  =  Log [10,  Min  [Abs  [Take  [Transpose [pts]  pi ,  n]  ]  ]  ] ; 

V  =  Map[{#|ri]I,  lf[#I2]]  s  0,  10'S  Log [10 ,  #[[2]]]]}  &,  Abs [Take [pts ,  n]  ]  ]  ? 
{x0,  xl}  =  {v[[l,  II,  v[[-l,  II}; 

pi  =  lgplt  [Map[  {#111 ,  10*121}  &,  v]  ,  DisplayFunction Identity, 

PlotStyle  -»  {Cornf lowerBlue}  ,  GridLines  -*  Automatic]  ? 

{{xa,  xb},  {ya,  yb}}  =  PlotRange  /.  FullOptions[pl] ; 
tks  =  logTks[10y\  10yb]  ; 

tksX = Prepend [ (GridLines /.  FullOptions [pi] )  III ,  {0,  {GrayLevel [0]  } } ]  /. 
{x_,  Y_>-  (x,  {dsh}}; 

tksX =  Join[ tksX,  Map[{#,  {colRotor} }  &,  tailRotorFreq*  tail]  ] ; 
tksY  =  DeleteCases [tks,  {x_,  ""}]  /.  {x_,  y_}  -+  {x,  {dsh}}; 
fn =  Simplify [Fit [v.  Table [chebPoly [i,  x] ,  {i,  0,  polyOrd}],  x]  ] ; 
p2  =  Plot[fn,  {x,  0,  vj-l,  II}, 

PlotStyle -»  {Thickness [0 . 005]  }  ,  Di spl ay Funct ion -*  Identity] ; 

Show [{pi,  p2,  Graphics [Map [rotorTxt,  tail]]}, 

DisplayFunction  -*  $DisplayFunction, 

FrameTicks -+ {Automatic,  tks.  Automatic,  tks/.  {{x_,  y_}  -*  {x,  ""}}}, 
GridLines-*  {tksX,  tksY},  Text  Style -*  {FontFamily -»  "Courier"}] 

] 


■  Importing  Data  (Tape  24  Run  82) 


In[10]:=  Dimensions [dat2 4  = ReadList [ "stbnbmlrT024R082 .datn] ]  1 

Out  [1 0]  =  {12907}  ^ 


In[ll]  ;=  f  d24  =  Module  [  {v,  n,  At  =  1  /  832  .}, 

v  =  2  Fourier  [dat24,  FourierParameters -»  {-1,  1}]; 
v  =  ReplacePart[v,  vpI/2,  1]; 
n  =  Round  [Length  [v]  /2.]  -1; 

Transpose [ {Range [0,  — - — ,  — - — ],  Take[v,  n+1]}] 
]> 


In  [12]  :  =  picFullFreq24  =  pltDat[fd24,  Length[fd24]  ,  20,  {2,  3,  4,  6,  8,  9,  10,  11,  12}]  ^ 

In  [13]  :  =  maxAmp  [fd24,  360,  370]  ^ 

Out  [13]  =  {24.7477,  367.457} 
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Figure  E4:  Mathematica  code  for  frequency  plots,  Page  2. 
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In  [14]  picInitPreq24  = 

pltDat [£d24,  1450,  20,  Join[0.2172  {1,  2,  3,  4,  6,  8,  10),  {2,  4}]] 


L  S'*  v*.  » 


40 

Frequency 


Out  [14]“  -Graphics- 


in [15]  ;=  maxAmp [£d2 4,  60,  70] 

Out  [15]  “  {126.822,  62.2743} 

In  [16]  :*=  Export ["freqFull.eps",  picFullFreq24,  "EPS"]; 

Export ["freqlnit.eps",  pic!nitFreq24,  "EPS"] ; 


Importing  Data  (Tape  45  Run  44) 


Dimensions [dat4 5  =  ReadList [ "stbnbmlrT045R044 .dat"] 


f d45  =  Module [{v,  n.  At  *  1  /  832  .), 

v  s  2  Fourier  [dat45,  Fourier  Parameters  -*•  {- 1,  1}]; 
v  tt  ReplacePart[v,  v [[1]]  /  2,  1] ; 
n  =  Round  [Length  [v]  /  2.)  -  1 ; 

Transpose  [{Range  [0,  — - ,  — - —  ],  Take[v,  n+1]}] 

L  L  1  2  At  2  n  At 

]> 

picFullFreq45  =  pltDat  [fd45,  Length  [  fd45]  ,  20,  {2,  3,  4,  6,  8,  9,  10,  11,  12}] 
maxAmp  [fd4 5,  3  60,  370] 

{26.3868,  367.554} 

pic!nitFreq45  =  pltDat  [£d45,  1450,  20,  Join[0.2172  {1,  2,  3,  4,  6,  8),  {2,  4}]] 
maxAmp  [  fd4  5 ,  60,  70] 

{78.355,  62.3204} 
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Figure  E5:  Mathematica  code  for  frequency  plots,  Page  3. 
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stabFreq.nh 


4 


In  [24]  :  =  Export [ " f RollFul . eps " ,  picFullFreq45,  "EPS"]; 

Export ["fRollIni .eps" ,  pic!nitFreq45,  "EPS"]; 


^1 


■  Importing  Data  (Tape  69  Run  42) 

In  [26]  : -  Dimensions [dat69  =  ReadList [ "stbnbmlrT069R042 .dat"  ]  ] 

Out  [26]  =  {68320} 

In  [27]  :  =  fd6  9  =  Module  [{v,  n,  At  =  1  /  832.}, 

v = 2  Fourier [dat69,  FourierParameters -♦  {-1,  1}]; 
v  =  ReplacePart  [v,  v[[l]]/2,  1]; 
n  =  Round  [Length  [v]  /  2  .  ]  -  1; 

Transpose [ {Range [0,  - „  - ],  Take  [v,  n.l]}] 

h 


T 

1 


In  [28] :  -  picFullFreq69  =  pltDat[fd69,  Length  [fd69]  ,  20,  {1,  2,  4,  8}] 

In[29}:=  maxAmp [ f d6 9 ,  360,  370] 

Out  [29]=  {26.5608,  367.603} 

In  [30]:=  picInitFreq69  = 

pltDat  [fd69 ,  7400,  20,  Join[0 .2172  {1,  2,  4,  6,  8,  12,  16,  20},  (1,  2,  4}]] 
In  [31]  :  =  maxAmp  [fd6  9,  60,  65] 

Out[3I]=  {14.5888,  62.2313} 

In  [32]  :  =  Export [ "fAprFul . eps" ,  picFullFreq69,  "EPS"]; 

Export ["fAprIni.eps",  pic!nitFreq69,  "EPS"]; 


U 

]J 

3 

]. 
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Figure  E6:  Mathematica  code  for  frequency  plots,  Page  4 
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Appendix  F  Bending  Stress  on  Wing  Panel  Lug 

In  this  section,  given  a  value  of  vertical  load  on  the  stabilator,  we  calculate  the  bending 
stress  on  the  top  lug  of  the  forward  spar.  We  want  to  compare  the  bending  stress  from 
two  separate  Sikorsky  reports  [1,  9]  with  our  predictions  based  on  measurements  from  the 
Flight  Loads  Survey. 

In  Sikorsky’s  fatigue  report  [1],  a  diagram  (on  page  H-10)  shows  the  strain  gauge 
locations  during  the  fatigue  test.  In  particular,  a  gauge  labelled  Number  2  (with  mnemonic 
STABL-2)  on  the  top  lug  of  the  forward  spar  appears  to  be  the  strain  gauge  used  to  check 
that  the  correct  test  loading  was  applied. 

In  Sikorsky’s  flight  loads  report  [9],  instead  of  the  bending  stress  on  the  forward  spar’s 
lug,  the  vertical  loading  on  the  stabilator’s  wing  panel  is  given.  The  geometry  calculations 
in  this  section  are  used  to  determine  an  equivalent  bending  stress  on  the  top  lug  of  the 
forward  spar,  which  is  where  the  bending  stress  from  Sikorsky’s  fatigue  report  [1]  was 
measured.  We  will  first  convert  this  vertical  load  into  a  bending  moment,  and  then  convert 
this  bending  moment  into  a  bending  stress. 

To  convert  a  bending  moment  to  a  bending  stress,  we  must  know  the  second  moment 
of  area  (see  Equation  (3.4)  on  page  16)  on  the  four  lugs  supporting  each  stabilator  wing 
panel.  Both  the  forward  and  aft  spar  have  two  lugs  located  on  the  top  and  bottom  of  the 
spars.  These  lugs  attach  the  stabilator’s  wing  panels  to  the  centre-box,  and  they  form 
the  male  part  of  the  clevis  lug  connection.  Figures  A4  and  A12  (on  pages  55  and  63) 
respectively  show  the  lugs  attached  to  the  forward  and  aft  spars.  From  these  figures  we 
can  determine  the  cross-sectional  geometry  of  the  lugs  for  use  in  second  moment  of  area 
calculations,  see  Figure  FI. 

From  the  strain  gauge  location  diagram  in  Sikorsky’s  fatigue  report  [1],  strain  gauge 
Number  2  appears  to  be  mounted  on  the  lug's  upper  surface,  away  from  the  lug’s  hole 
and  towards  the  spar’s  root.  This  mounting  (away  from  the  hole)  would  reduce  the  stress 
concentration  effect  produced  by  the  lug’s  hole.  We  see  from  Figure  A4  (and  more  clearly 
from  Figure  FI)  that  between  the  lug’s  hole  and  the  spar’s  root,  the  lug  has  a  varying 
cross-sectional  area,  which  is  a  function  of  the  distance  £.  This  distance  is  measured  from 
the  spars  root  (at  BL  9)  towards  the  lug’s  hole,  see  Figure  FI.  To  avoid  the  bending  gauge 
being  located  above  the  forward  spar  or  above  the  hole  of  the  forward  spar’s  lug,  we  bound 
the  distance  £  within  the  range 

0  <  £  <  0.884. 

For  ease  of  calculation,  the  cross-sectional  area  of  the  forward  spar’s  lug  is  partitioned 
into  two  regions.  These  two  regions  are  the  constant  and  fillet  regions,  which  are  denoted 
respectively  by  the  subscripts  1  and  2  on  the  centroid  symbols.  (The  lug’s  fillet  region 
is  shown  cross-hatched  in  Figure  FI.)  The  height  of  this  fillet  region  is  a  function  of  the 
distance  £,  and  may  be  calculated  from  the  geometry  of  a  square  inscribed  in  a  circle.  The 
radius  of  the  fillet  is  1.00",  and  so  the  fillet’s  height  is  given  by  the  relation 

»/  =  1  -  \/£(2  -  £)• 

The  distance  from  the  fillet’s  centroid  to  the  neutral  axis  is  given  by 

Af  =  1.59-77/2. 
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The  angle  the  aft  spar’s  lug  makes  with  the  horizontal  (denoted  by  d)  was  not  shown 
in  any  available  drawing,  and  so  it  was  calculated  from  the  given  geometry.  Figure  FI 
shows  the  quadrilateral  OAJBC^  where  the  point  T)  lies  on  the  line  OC .  The  lines  *AD  and 
BC  are  parallel,  and  are  both  perpendicular  to  the  line  OC.  Using  simple  geometry,  an 
equation  in  terms  of  d  is  obtained,  1.866  sin  d  =  0.78  +  0.145  cos  ,  whose  only  solution 
in  the  range  0  <  d  <  90°  is  ~  28.7°.  The  centroid  of  the  Region  3  is  a  distance 

Aa  =  2.500  -  (1.196  -  0  tan 
from  the  spar’s  centreline  of  symmetry. 

A  rectangle  whose  centroid  is  a  distance  y  from  the  neutral  axis  has  a  second  moment 
of  area  given  by  (see  Roark  [28])  the  equation 

1  =  bd3/l2  +  bdy2, 

where  b  and  d  are  the  rectangle’s  width  and  height,  respectively.  The  second  moment  of 
area  for  Regions  1  and  2  of  the  forward  spar’s  lug  are  then 

h  =  5.27  in4 


and 

I2  =  n  (0.037??2  +  0.444A2 ) 

=  0.565  +  £(0.262£  -  0.524)  +  0.148(£  +  0.430)(£  -  2.43)v/^(2-0  in4. 

The  second  moment  of  area  for  Region  3  of  the  aft  spar’s  lug  is 

h  =  0.181  +  0.685 A2 

=  2.51+^(0.205^  +  1.38)  in4. 

Summing  I\ .  /2,  and  1%  and  multiplying  by  two  (to  account  for  symmetry)  gives  the 
total  second  moment  of  area  of  the  stabilator  wing  panel  lugs  a  distance  £  inboard  from 
the  spar  roots  (at  BL  9.00) 

Tvp  =  2  (/i  +  I2  +  h) 

=  16.7  +  £(0.935£  +  1.72)  +  0.296(£  +  0.430)^  -  2.43)  ^(2  -  C)  in4- 

This  curve  is  close  to  linear  with  a  quick  turn  at  £  =  0.  The  minimum  and  maximum 
second  moments  of  area  are 

min  (7wp)  =  16.7  in4  (which  occurs  at  £  =  0.0184)  (FI) 

and 

max  (7wp)  =  18.4  in4  (which  occurs  at  =  0.884).  (F2) 
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From  Sikorsky’s  loads  report  [9,  p.  47],  the  maximum  vertical  load  on  the  right  stabi- 
lator  from  the  flight  loads  measurements  was  Vnmx  =  1848  lb.  Assume16  that  this  vertical 
load  acts  through  the  horizontal  centroid  of  an  elliptical  lifting  distribution.  If  the  spar’s 
length  is  L  =  74  in,  then  the  bending  moment17  at  the  spar’s  root  is  given  by 

max  (Mr)  =  4Vrmax£/(37r) 

=  58xl03  in  lb.  (F3) 

For  the  forward  spar,  the  height  from  the  centreline  of  symmetry  to  the  top  of  lug  is 
hf  =  3.41",  see  Figure  FI.  While  from  Equation  (3.4),  the  bending  stress  at  the  top  of  the 
forward  spar’s  lug  is  given  by  aY  =  Mrftf/Jwp.  Using  the  maximum  bending  moment  given 
by  Equation  (F3)  and  the  range  of  second  moment  of  area  given  by  Equations  (FI)  and 
(F2),  yields  the  range  of  maximum  bending  stress  at  the  top  of  the  forward  spar’s  lug 

11  ksi  <  max(<rr)  <  12  ksi.  (derived  from  vertical  load  on  stabilator).  (F4) 

A  range  of  maximum  stresses  was  obtained  because  the  second  moment  of  area  varied 
depending  on  where  the  bending  strain  gauge  was  located. 

Sikorsky  [1,  p.  H-16]  appear  to  have  carried  out  the  fatigue  testing  of  the  stabilator ’s 
wing  panel  at  a  steady  stress  of  6.00  ksi  and  vibratory  stresses  of  3.75-7.50  ksi.  These 
steady  and  vibratory  stresses  sum  to  a  maximum  bending  stress  at  top  of  the  forward 
spar’s  lug  of 


13.5  ksi  (obtained  from  fatigue  testing  on  stabilator).  (F5) 

We  now  see  that  the  maximum  bending  stresses  given  by  Equations  (F4)  and  (F5)  (that 
is,  the  two  different  Sikorsky  reports)  are  comparable. 


16 The  assumptions  made  in  this  appendix  are  not  numbered  because  they  are  not  part  of  the  main 
analysis. 

1 7 The  horizontal  centroid  (of  the  quarter-ellipse)  occurs  at  a  distance  4L/(37t)  from  the  root,  see  page  16. 
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Appendix  G  Fractographic  Analysis 

In  this  appendix,  we  examine  the  cracked  lug  fractographically  in  order  to  obtain  a 
simple  lower  bound  on  the  number  of  cycles  to  failure. 

Five  regions  from  the  cracked  lug’s  fracture  surface  were  chosen  to  estimate  the  overall 
crack  growth  rate.  Figure  G1  shows  the  cracked  lug,  as  well  as  the  five  regions  inspected 
more  closely.  For  the  location  of  this  lug  tip  within  the  stabilator  see  Figure  1.3  (shown 
on  page  4). 


Figure  Gl:  Upper  section  of  cracked  lug  tip  showing  positions  of  the  five 
selected  regions  for  striation  sampling.  The  origin  and  positive  direction 
of  the  x  and  y  coordinates  are  shown  by  the  axes.  (Unlabelled  photograph 
courtesy  of  Rohan  Byrnes.) 


In  each  of  these  five  inspection  regions,  the  coordinates  of  the  visible  striations  were 
measured  in  three  coordinate  x,  y ,  and  z.  The  x  and  y  coordinates  lie  on  the  crack 
propagation  plane  (see  Figure  Gl).  The  z  coordinate  was  measured  perpendicular  to  the 
crack  propagation  plane,  and  was  determined  using  the  focal  length  of  the  microscope. 
The  measured18  coordinates  (which  are  in  millimetres)  are  shown  in  Table  Gl. 

These  regions  were  chosen  in  an  approximately  straight  line  from  the  crack  initiation 
site  (the  x-y  origin  in  Figure  Gl)  to  the  crack  edge  (which  is  near  Region  E).  Crack 
striations  near  the  origin  were  difficult  to  resolve,  so  the  first  region  sampled  (Region  A) 
was  approximately  10  mm  from  the  origin.  In  order  to  determine  the  coordinates  of  the 
striations,  the  sampled  regions  were  magnified  1000  times  (xlOO  for  the  lens  and  xlO 

18 Those  striation  coordinates  were  measured  by  Rohan  Byrnes  at  DSTO. 
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Table  Gl:  Coordinates  of  sampled  striations  from  five  selected  regions . 
(All  coordinates  are  in  millimetres. ) 
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for  the  eye  piece).  The  resolution  of  coordinate  measurements  and  of  the  motion  in  the 
microscope’s  stage  was  0.25  /mi. 

One  way  to  estimate  the  crack  growth  rate  is  to  determine  the  distance  between  visible 
striations.  In  other  words,  if  a  denotes  the  crack  length  from  the  crack  initiation  point, 
then  we  can  estimate  the  crack  growth  rate  (per  loading  cycle)  as 


da  ^  A  a. 

In  k  an' 


(Gl) 


where 


(G2) 


is  the  perpendicular  distance  between  visible  striations  and  AN  is  the  number  of  loading 
cycles  between  these  striations.  We  need  to  make  an  assumption  about  the  number  of 
cycles  between  visible  crack  striations  in  order  to  estimate  the  number  of  cycles  it  took  to 
completely  crack  the  lug: 


<^>  Assumption  17  There  is  only  one  loading  cycle  between  visible  crack  striations,  that 

is,  A N=1  in  Equation  (Gl). 


In  other  words,  at  a  magnification  of  xlOOO,  only  loading  cycles  that  cause  these  visible 
crack  striations  are  considered.  This  assumption  means  that  any  estimate  we  obtain  for 
the  number  of  cycles  to  failure  has  to  necessarily  be  a  lower  bound  on  the  true  value.  This 
statement  follows  from  the  fact  that  there  are  more  striations  than  could  be  seen  at  a 
magnification  of  xlOOO.  ■  , 

The  crack’s  total  length  (which  is  approximately  16  mm)  can  only  be  measured  in 
the  plane  of  the  crack,  that  is,  in  the  x-y  plane.  (It  would  be  impossible  to  measure  the 
millions  of  out-of-plane  surface  variations  along  the  crack’s  length.)  Hence  when  measuring 
crack  growth  we  ignore  the  variation  in  the  z  coordinate  both  (i)  between  striations  (in 
Equation  (G2))  and  (ii)  of  the  total  crack  length. 

Using  Equation  (Gl),  the  crack  growth  rate  can  be  estimated  using  the  measured 
coordinates  of  the  visible  crack  striations.  If  (aq,  2/o  -h)  represents  the  coordinates  of  the 
ith  measured  striation,  then 


since  AN -l  from  Assumption  17.  Estimates  of  the  crack  growth  rate,  calculated  using 
Equation  (G3),  are  plotted  in  Figure  G2.  Note  that  the  crack  growth  A  a  and  crack  growth 
rate  Aa/AN  are  equivalent  for  these  striation  measurements  because  we  have  assumed 
that  A7V=1. 

In  Figure  G2,  the  box-plots  show  the  5th,  25th,  50th,  75th,  and  95th  percentiles  of 
the  relevant  data.  The  5th  and  95th  percentiles  are  denoted  by  the  lower  and  upper 
parts,  respectively,  of  the  box-plot’s  whiskers.  The  25th  and  75th  quartiles  are  denoted 
by  the  lower  and  upper  sections,  respectively,  of  the  box-plot’s  rectangle.  The  median  (or 
50th  percentile)  is  denoted  by  the  horizontal  bar  through  the  mid-region  of  the  box-plot’s 
rectangle,  which  extends  through  the  domain  of  the  relevant  data.  The  solid  black  line 
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Figure  G2:  Estimate  of  crack  growth  rate  from  five  selected  regions .  Box- 
plots  show  5th ,  25th,  50th,  15th,  and  95th  percentiles .  Horizontal  dashed 
lines  and  dashed  box-plots  denote  median  and  percentiles,  respectively, 
of  individual  regions.  Solid  black  line  and  solid  box-plot  denote  median 
and  percentiles,  respectively,  of  five  regions  aggregated. 


through  all  five  regions  and  the  solid  box-plot  on  the  left  denote  the  median  and  percentiles 
of  the  aggregated  data  from  all  five  regions. 

We  can  estimate  the  number  of  cycles  to  failure  by  dividing  the  total  crack  length  by 
the  crack  growth  rate,  but  first  we  need  to  make  an  assumption  about  the  way  in  which 
the  crack  grows: 

Assumption  18  The  crack  grows  linearly  (that  is,  at  a  constant  rate)  throughout  its 
length. 

In  reality,  the  crack  growth  rate  probably  follows  some  sort  of  power  law  as  a  function 
of  crack  length.  For  example,  a  simple  law  for  crack  growth  rate  is  given  by 

%  =  C(AK)m,  (G4) 

where  C  is  a  coefficient,  A K  is  the  range  of  the  stress  intensity  factor,  and  m  is  the 
characteristic  slope  for  the  curve  of  the  crack  growth  rate.  (See  Megson  [20]  or  the  fracture 
mechanics  ESDU  data  sheet  [21]  for  further  details.)  Both  C  and  m  are  dependent  on  the 
material  properties  and  the  mean  stress  under  consideration. 

Let  Sm\n  and  5max  denote  respectively  the  minimum  and  maximum  values  for  stress 
within  a  cycle  of  a  constant  amplitude  loading.  In  its  simplest  form,  the  range  of  the  stress 
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intensity  factor  is  given  by 

A K  —  (S'max  —  S'min)  \J  TUI  Q (G5) 

where  a  is  a  non-dimensional  length  coefficient,  which  usually  expressed  as  the  ratio  of 
crack  length  to  any  convenient  local  dimension  in  the  plane  of  the  component. 

Combining  Equations  (G4)  and  (G5),  we  see  that  the  crack  growth  rate  (given  by  this 
simple  law)  is  a  power  function  of  the  crack  length.  The  data  shown  in  Figure  G2  suggests 
that  the  crack  growth  rate  is  relatively  constant,  and  hence  at  least  for  the  later  part  of 
the  crack  the  characteristic  slope  m  must  be  close  to  zero. 

From  Figure  G2  we  can  determine  the  median  crack  growth  rate  of  the  five  sampled 
regions,  as  well  as  the  five  regions  aggregated.  These  median  values  for  the  crack  growth 
rate  are  shown  in  Table  G2.  As  can  be  seen,  the  crack  growth  rate  appears  to  remain 
relatively  constant  (at  approximately  0.7  //m/cycle)  over  the  five  sampled  regions. 


Table  G2:  Median  of  crack  growth  rate  for  five  sampled  regions  and 


regions  aggregate a. 

Region 

ABODE 

Aggregate 

Crack  growth  rate  (//m/cycle) 

0.8  0.7  0.7  0.6  0.8 

0.7 

Using  Assumption  18,  we  can  divide  the  total  crack  length  (which  is  16  mm)  by  the 
median  crack  growth  rate  (which  is  0.7  /tm/cyclc)  to  obtain  a  lower  bound  on  the  number 
of  cycles  to  failure:  20  x  103  cycles. 

Dividing  the  number  of  flight  hours  before  the  crack  was  detected  (which  is  approx¬ 
imately  400  hours)  by  the  number  of  cycles  (which  is  20  x  103  striations),  we  obtain  the 
cycle  period: 

1  minute  (between  cycles  causing  visible  striations). 

In  obtaining  this  period  we  have  made  an  implicit  assumption  about  when  the  crack  began: 


Assumption  19  The  crack  began  when  the  component  was  installed  (which  was 
402  flight  hours  before  detection). 


If  the  crack  began  significantly  after  installation,  then  the  time  between  cycles  given  above 
would  be  significantly  shorter. 

This  1  minute  period  between  cycles  suggests  that  at  least  the  visible  striations  were 
caused  by  loading  that  occurred  approximately  once  per  manoeuvre. 
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